Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB143	2,094,808	A→C	41.3%	L186R (CTG→CGG)	CJ019_01995 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB143	2,094,808	0	A	C	41.3%	1.2 / 3.6	33	L186R (CTG→CGG)	CJ019_01995	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (9/8);  new base C (7/7);  total (16/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCTGCCGGGGCCATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATAC‑TGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAG  >  ACB143/2094675‑2094936                                                                                                                                       |                                                                                                                                 ttctGCCGGGGCCATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATAC‑TGCACGCCGCGCTCGAAGGCCCgg                                                                                                                                >  4:183872/1‑136 (MQ=255)                                 cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTTCGCAAATAC‑TGCACGCCGCGCTCGAAGGCCCCGCGCGCCAGGCCCGGCAGGCGGTTTTTGAgggg                                                                                                >  7:55640/1‑135 (MQ=255)                                 cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATAC‑TGCACGCCGCGCTCGAAGGCCACGCGCCCCAGGCCGGGCTGGCCATTTTGCAGGGc                                                                                                >  6:281877/1‑136 (MQ=255)                                 cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATAC‑TGCACGCCGCGCTCGAAGGCCACGCCCGCCCGGCCGGGCGGGCGGGTTTTCAgggg                                                                                                >  3:83533/1‑135 (MQ=255)                                 cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCAGCTCTTCATACACCTTCGGCGCCGGCAGCTTGCGCAACTAC‑TGCACTCCGCGCGCGAAGGCCCAGCGCGCCGGGACCGGCTGGCCGTTTTgtcatgc                                                                                                >  8:270247/1‑130 (MQ=255)                                      gCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGTCTCCGGCTGCTTGCGCAAATAC‑TGCACGCCGCGCTCGAAGGCCCTGCGCGCGGGGCCCGGGAGGGCGTTTTGACGGGGCTggt                                                                                           >  2:186010/1‑134 (MQ=255)                                      gCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATAC‑TGCACGCCGCGCTCGAAGGCCCCGCGCGCCGGGGCGGGCGGGGCGGTTTGGAGGGcggggg                                                                                           >  3:132095/1‑131 (MQ=255)                                              cAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATAC‑TGCACGCCGCGCTCGAAGGCCCGGCGCGCCGGGCCGGGGGGGGCATTTTGGAGGGCCTGGTGGGCGAgg                                                                                   >  6:22150/1‑136 (MQ=255)                                                   gTTTCACATGATAGGTAAGGTGCTTTTCATACACGTTCGGCTCCGGGTGGTTGCTCAAATAC‑TCCACGCCGCGCTGGAAGGCCAGGGGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGgttgtt                                                                              <  2:179444/136‑1 (MQ=255)                                                   gTTTCACAAGATACTTCAGCTGCTCTACATACACCTTCGGCGCCGGCTGCTTGCCCAACTAC‑TCCACGCCGCGCGCGAAGGCCTGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGgttgtt                                                                              <  2:175446/136‑1 (MQ=255)                                                   gTTGCTCAAGATACGTCTGCTGCTCTTCATACACCTTCGGCGCCGGCTGCTTGCGCAAATAC‑TGCACGCCGCGCGCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGgttgtt                                                                              <  2:41191/136‑1 (MQ=255)                                                    tttCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATAC‑TGCACGCCGCGCTCGAAGGCCACGCGCGCCTGGCCCGGCTGGCCGTTTTGCAGGGCCTGGTGGCCGAGggtgtt                                                                              >  3:117463/1‑135 (MQ=255)                                                     ttgtcAATATCCGTCAGCTGCTCTTCAAACAGATTGGGGTCGGGGTGCTTTC‑AAAACACATGCCCCCGGCGCTAGAAGGCCAGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGttgttg                                                                             <  6:41456/132‑1 (MQ=255)                                                      tctcAAGATACGGCACCTGCTATCAAAACACCTTCGGCTCCGGCTGGTTGCCAAAATAC‑TCCCCGCCGCGCGCGACGGCCCGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTga                                                                            <  8:88898/135‑1 (MQ=255)                                                         cAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATAC‑TGCACGCCGCGCTCGAAGGCCAGGCCCGCCTGGCCCGGCGGGCCATTTTGAAGGGCCAGTTGGCCGAGGTTGTTgaaga                                                                         >  8:96483/1‑135 (MQ=255)                                                                    aGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATAC‑TGCACGCCGCGCTCGACGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATgc                                                              <  8:269084/135‑1 (MQ=255)                                                                     gctgctCTTCATACACCTTCGGCTCCGGCTGCTTGCGAAAATAC‑TGCCCGCCGCGCTCGAAGGCCCGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATgcag                                                            <  8:156388/136‑1 (MQ=255)                                                                        gatctTCATACACGTTCCGCTCGGGCTGCTTGCGAAAATAC‑GGCCCGCGGCGCTCGAAGGCCCGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATgcagcag                                                         <  1:93640/134‑1 (MQ=255)                                                                         ctctTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATAC‑TGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATgcagcag                                                         <  3:149941/135‑1 (MQ=255)                                                                         ctctTCAAACACCTTCGGCTCCGGCTGCTGCCGAAAATAC‑TGCACGCCGCGCCCGAAGGCCTGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATgcagcag                                                         <  6:218162/135‑1 (MQ=255)                                                                            ttCATACACCTTCGGCTCCGGCTGCTTGCGCAAATAC‑TGCACGCCGCGCTCGAAGGCCACGCGCGCCTGGCCGGGCTGGCCTTTTTTCAGGGCCTGTTGGGCGAGGTTGTTGAAGAACACTATTTGCAGCAGCAcc                                                     >  6:242942/1‑136 (MQ=255)                                                                                      ttCGGCTCCGGCTGCTTGCGCAAATAC‑TGCACGCCGCGCTCGAAGGCCAGGCGCGCCCGGCCGGGCGGGCCGGTTTTGAGGGCCGGGTGGGCGCGGTTGGTTAAGAAATCTATATGAAGGAGCCCCCCGATATGGc                                           >  6:29980/1‑136 (MQ=255)                                                                                        cGGCGCGGGCTGGTTGCCAAAATAC‑CGCCCGCCGCGCTCGCCGGCCCGGCGGGCGGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGc                                           <  4:149403/134‑1 (MQ=255)                                                                                          gCTCCGGCTGCTTGCGCAAATAC‑TGCACGCCGCGCTCGAAGGCCCGGCGCGCCGGGGCGGGGTGGGCGTTTTTCAGGGCCGGTTGGGCGAGGTTGTTGTAGAACTCGATATGCAGCAGCACCAGGATATGGCGGAt                                       >  2:169278/1‑136 (MQ=255)                                                                                            cccGGCGGCTTGCCAAAATAC‑TCCACGCCGCGCTCGAAGCCCCGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGAttt                                     <  3:223859/135‑1 (MQ=255)                                                                                                       gcaaaaaTAC‑CCCCCGCCGCGCGCCACGGCCAGCCGCGCGGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCaa                          <  2:156009/132‑1 (MQ=255)                                                                                                                     ccgccgcGCTCGACGGCCAGGGGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGa             <  3:52750/135‑1 (MQ=255)                                                                                                                      cgccgcGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCGTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGa             >  7:212710/1‑135 (MQ=255)                                                                                                                       gccgcGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGGTGGTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATTCTTTGCCTCGACa           >  1:171542/1‑136 (MQ=255)                                                                                                                           cgcTCGAAGGCCCGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCgg       <  8:157014/136‑1 (MQ=255)                                                                                                                            gcTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCGTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATTGCGGATTTCCTTTCCCCAATGCTTGGCCTCGACAGCggg      >  5:166443/1‑136 (MQ=255)                                                                                                                             cgcGACGGCCCGGCGCGCGGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGc     <  3:9120/134‑1 (MQ=255)                                                                                                                                gAAGGCCAGGCGGGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCaag  <  6:205886/136‑1 (MQ=255)                                                                                                                                       |                                                                                                                                 TTCTGCCGGGGCCATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATAC‑TGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAG  >  ACB143/2094675‑2094936

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑