Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	2,719,426	T→C	27.4%	intergenic (‑321/+2)	kgtP ← / ← rrfG	alpha‑ketoglutarate transporter/5S ribosomal RNA of rrnG operon

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	2,719,426	0	T	C	27.4%	39.9 / 10.6	28	intergenic (‑321/+2)	kgtP/rrfG	alpha‑ketoglutarate transporter/5S ribosomal RNA of rrnG operon
Reads supporting (aligned to +/- strand):  ref base T (10/10);  new base C (4/4);  total (14/14)
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

CGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCA  >  CP009273/2719281‑2719562                                                                                                                                                  |                                                                                                                                         cGTTTGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGAt                                                                                                                                       <  1:188602/149‑1 (MQ=255)   tATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAGTGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTCGATGc                                                                                                                                     >  1:109272/1‑149 (MQ=255)       gAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCCAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTgg                                                                                                                                 >  1:27615/1‑149 (MQ=255)                   ttGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTCGATGCCTGGCAGTTc                                                                                                                           >  2:136191/1‑143 (MQ=255)                        gAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGc                                                                                                                >  1:126375/1‑149 (MQ=255)                         aaGCCGTATCCGTTGCTGAATCTGGCATATGTGGGGGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCa                                                                                                               <  1:20666/149‑1 (MQ=255)                                  ccGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCGAAAAGGCCCTCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGAc                                                                                                      <  1:129908/149‑1 (MQ=255)                                         tGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCACATGGGTAGAcc                                                                                                     <  2:92602/143‑1 (MQ=255)                                                 gCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCGCACTACCATCg                                                                                       >  1:102321/1‑149 (MQ=255)                                                       gtgGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTCGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTa                                                                                 <  1:155515/149‑1 (MQ=255)                                                           gAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGgcgc                                                                                   >  2:59932/1‑143 (MQ=255)                                                                            cGTCGCATCAGGCATCTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTCGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTc                                                            <  1:136459/149‑1 (MQ=255)                                                                             gTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCg                                                           <  1:38599/149‑1 (MQ=255)                                                                             gTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTGCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCg                                                           <  1:38628/149‑1 (MQ=255)                                                                              tCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGGGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCgg                                                          >  1:90108/1‑149 (MQ=255)                                                                                    cAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTCGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCgg                                                          <  2:96556/143‑1 (MQ=255)                                                                                         aTTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAg                                               >  1:2417/1‑149 (MQ=255)                                                                                                 catcaaattatgcgaaaGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGGGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGTCATGGGGTCAGGt                                             >  2:154869/15‑143 (MQ=17)                                                                                                   gcgcCAAGGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCAGAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTgg                                           <  2:83335/143‑1 (MQ=255)                                                                                                             aaaaaGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTGACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCCCcgcg                                 >  2:20854/1‑143 (MQ=255)                                                                                                                   gCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCa                     >  1:305573/1‑149 (MQ=255)                                                                                                                     cATCCGTCAGGATGGCCTTTCGCATAATACGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAgg                   >  1:216398/1‑149 (MQ=255)                                                                                                                        ccGTCAGGATAGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGgccgcc                      <  2:234755/143‑1 (MQ=255)                                                                                                                               gATGGCCTTTCGCATAATTCGATGCCTGGCAGTCCCCGACTCTCGCATGGGGAGACCCCACACTCCCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGt         <  1:271928/149‑1 (MQ=255)                                                                                                                                aTGGCCTTTCGCATAATTCGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAAt              >  2:57071/1‑143 (MQ=255)                                                                                                                                    ccTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTAt    <  1:237323/149‑1 (MQ=37)                                                                                                                                      tttCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCa  <  1:67343/149‑1 (MQ=35)                                                                                                                                            aTAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACGACCGCGCTACGGCCGGCGGGCAGATTCTGTTTTGTCa  >  2:236816/1‑143 (MQ=11)                                                                                                                                                  |                                                                                                                                         CGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCA  >  CP009273/2719281‑2719562

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑