Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	637,033	G→T	19.3%	V97V (GTC→GTA)	uspG ←	universal stress protein UP12

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	637,033	0	G	T	19.3%	77.9 / 4.9	31	V97V (GTC→GTA)	uspG	universal stress protein UP12
Reads supporting (aligned to +/- strand):  ref base G (11/14);  new base T (3/3);  total (14/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 = 2.85e-01

GGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGAACCTAACAGATGGGTCGAAATCGATGGGTTGCGAGAACCAATAACTACAACATCAGCCCCCAGTTCTTCTGCCAACTCATTGACTTCATCCCGCACGCTACCAAAACGGACATGTTGTTTAATGCGGGAAGGATCGATGGTGAAGTGGCTGACCATCGTTTGCAGACGTTCTTGTGCTTCATGTTGCAGATGCTCTTCAAAACGACGCACATCAGCGGCAAAACGGTGCAGGCTCAGGCTGGCTGACCCGGGTAGTACGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCATTTCAAA  >  CP009273/636813‑637290                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                  ggTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACNACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGAACCTAACAGATGGGTCGAAATCGATGGGTTGCGATAACCAATAACTACAACATCAGCCCCCAGTTCTTCTGCCAACTCATTGACTTCATCCCGCACGCTACCAAAACGTACATGTTGTTTAATGCGGGAAGGATCg                                                                                                                                                                                                             <  1:526692/275‑1 (MQ=255)                                          cggggATGGCTCCGGGTAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGAACCTAACAGATGGGTCGAAATCGATGGGTTGCGAGAACCAATAACTACAACATCAGCCCCCAGTTCTTCTGCCAACTCATTGACTTCATCCCGCACGCTACCAAAACGGACATGTTGTTTAATGCGGGAAGGATCGATGGTGAAGTGGCTGACAATCGTTTGCAGACGTTCTtgtg                                                                                                                                                                     <  1:319667/272‑1 (MQ=255)                                                        ggCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGAACCTAACAGATGGGTCGAAATCGATGGGTTGCGAGAACCAATAACTACAACATCAGCCCCCAGTTCTTCTGCCAACTCATTGACTTCATcc                                                                                                                                                                                                                                                          >  2:326076/1‑175 (MQ=255)                                                            gggCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGAACCTAACAGATGGGTCGAAATCGATGGGTTGCGAGAACCAATAACTACAACATCAGCCCCCAGTTCTTCTGCCAACTCATTGACTTCATCCCGCa                                                                                                                                                                                                                                                      <  1:326081/175‑1 (MQ=255)                                                                                              aagagagacagaGATTGGCGTTGCGGATCCCGCTCGAGGCGTTAGATCCTAACAGATGATTCGTAATCGATGGGTTGCGAGAACCAATAACTACAACATCAGCCCCCAGTTCTTCTGCCAACTCATTGACTTCATCCCGCACGCTACCAAAACGGACATGTTGTTTAATGCGGGAATGATCGATGGTGAAGTGGCTGACCATCGTTTGCAGACGTTCTTGTGCTTCATGTTGCAGATGCACTTCAAAACGACGCACATCAGCGGCAAAACGGTg                                                                                                                 <  2:568767/263‑1 (MQ=255)                                                                                                      ggCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGAACCTAACAGATGGGTCGAAATCGATGGGTTGCGAGAACCAATAACTACAACATCAGCCCCCAGTTCTTCTGCCAACTCATTGACTTCATCCCGCACGCTACCAAAACGGACATGTTGTTTAATGCGGGAAGGATCGATGGTGAAGTGGCTGACCATCGTTTGCAGACGTTCTTGTGCTTCATGTTGCa                                                                                                                                                         >  2:490973/1‑226 (MQ=255)                                                                                                          gATTGGCGTGGCGGATCACGCTCGAGGCGTTAGAACCTAACAGATGGGTCGAAATCGATGGGTTGCGAGAACCAATAACTACAACATCAGCCCCCAGTTCTTCTGCCAACTCATTGACTTCATCCCGCACGCTACCAAAACGGACATGTTGTTTAATGCGGGAAGGATCGATGGTGAAGTGGCTGACCATCGTTTGCAGACGTTCTTGTGCTTCATGTTGCAGATg                                                                                                                                                     <  1:490981/226‑1 (MQ=255)                                                                                                                    gcgGATCACGCTCGAGGCGTTAGAACCTAACAGATGGGTCGAAATCGATGGGTTGCGAGAACCAATAACTACAACATCAGCCCCCAGTTCTTCTGCCAACTCATTGACTTCATCCCGCACGCTACCAAAACGGACATGTTGTTTAATGCGTGNAGGATCGATGGTGAAGTGGCTGACCATCGTTTGCAGACGTTCTTGTGCTTCATGTTGCAGATGCTCTTCAAAACGACGAACATCAGCGGCAAAAc                                                                                                                     >  1:568778/1‑248 (MQ=255)                                                                                                                                                aaCAGATGGGTCGAAATCGATGGGTTGCGAGAACCAATAACTACAACATCAGCCCCCAGTTCTTCTGCCAACTCATTGACTTCATCCCGCACGCTACCAAAACGGACATGTTGTTTAATGc                                                                                                                                                                                                                        >  2:74805/1‑121 (MQ=255)                                                                                                                                                    gATGGGTCGAAATCGATTGGTTGCGAGAACCAATAACTACAACATCAGCCCCCAGTTCTTCTGCCAACTCATTGACTTCATCCCGCACGCTACCAAAACGGACATGTTGTTTAATGCGGGa                                                                                                                                                                                                                    <  1:74805/121‑1 (MQ=255)                                                                                                                                                        ggTCGAAATCAATGGGTTGCGAGAACCAATAACTACAACATCAGCCCCCAGTTCTTCTGCCAACTCATTGACTTCATCCCGCACGCTACCAAAACGGACATGTTGTTTAATGCGGGAAGGATCGATGGTGAAGTgg                                                                                                                                                                                                 >  1:647821/1‑136 (MQ=255)                                                                                                                                                            gAAATCGATGGGTTGCGAGAACCAATAACTACAACATCAGCCCCCAGTTCTTCTGCCAACTCATTGACTTCATCCCGCACGCTACCAAAACGGACATGTTGTTTAATGCGGGAAGGATCGATGGTGAAGTGGCTGa                                                                                                                                                                                             <  2:647810/136‑1 (MQ=255)                                                                                                                                                                   aTGGGTTGCGAGAACCAATAACTACAACATCAGCCCCCAGTTCTTCTGCCAACTCATTTACTTCATCCCGCACGCTACCAAAACGGACATGTTGTTTAATGCGGGAAGGATCGATGGTGAAGTTGCTGACCATCGTTTGCAGACGTTCTTGTGCTTCATGTTGCAGATGCTCTTCAAAACGACGCACATCAGCGGCa                                                                                                                         >  1:571366/1‑197 (MQ=255)                                                                                                                                                                       gTTGCGAGAACCAATAACTACAACATCAGCCCCCAGTTCTTCTGCCAACTCATTTACTTCATCCCGCACGCTACCAAAACGGACATGTTGTTTAATGCGGGAAGGATCGATGGTGAAGTTGCTGACCATCGTTTGCAGACGTTCTTGTGCTTCATGTTGCAGATGCTCTTCAAAACGACGCACATCAGCGGCAAAAc                                                                                                                     <  2:571355/197‑1 (MQ=255)                                                                                                                                                                                aCCAATAACTACAACATCAGCCCCCAGTTCTTCTGCCAACTCATTGACTTCATCCCGCACGCTACCAAAACGGACATGTTGTTTAATGCGGGaa                                                                                                                                                                                                                   >  1:408027/1‑94 (MQ=255)                                                                                                                                                                                 ccAATAACTACAACATCAGCCCCCAGTTCTTCTGCCAACTCATTGACTTCATCCCGCACGCTa                                                                                                                                                                                                                                                 >  1:333833/1‑63 (MQ=255)                                                                                                                                                                                    aTAACTACAACATCAGCCCCCAGTTCTTCTGCCAACTCATTGACTTCATCCCGCACGCTACCAAAACGGACATGTTGTTTAATGCGGGAAGGAt                                                                                                                                                                                                               <  2:408022/94‑1 (MQ=255)                                                                                                                                                                                     tAACTACAACATCAGCCCCCAGTTCTTCTGCCAACTCATTGACTTCATCCCGCACGCTACCaa                                                                                                                                                                                                                                             <  2:333828/63‑1 (MQ=255)                                                                                                                                                                                     tAACTACAACATCAGCCCCCAGTTCTTCTGCCAACTCATTGACTTCATCCCGCACGCTACCAAAACGGa                                                                                                                                                                                                                                       >  2:597454/1‑69 (MQ=255)                                                                                                                                                                                         tACAACATCAGCCCCCAGTTCTTCTGCCAACTCATTGACTTCATCCCGCACGCTACCAAAACGGACAtg                                                                                                                                                                                                                                   <  1:597465/69‑1 (MQ=255)                                                                                                                                                                                                tCAGCCCCCAGTTCTTCTGCCAACTCATTTACTTCATCCCGCACGCTACCAAAACGGACATGTTGTTTAATGCGGGAAGGATCGATTGTGAAGTGGCTGACCATCGTTTGCAGACGTTCTTGTGCTTCATGTTGCAGATGCTCTTCAAAac                                                                                                                                          >  1:621257/1‑151 (MQ=255)                                                                                                                                                                                                   gCCCCCAGTTCTTCTGCCAACTCATTTACTTCATCCCGCACGCTACCAAAACGGACATGTTGTTTAATGCGGGAAGGATCGATGGTGAAGTGGCTGACCATCGTTTGCAGACGTTCTTGTGCTTCATGTTGCAGATGCTCTTCAAAac                                                                                                                                          >  2:364898/1‑148 (MQ=255)                                                                                                                                                                                                   gCCCCCAGTTCTTCTGCCAACTCATTGACTTCATCCCGCACGCTCCCAAAACGGACATGTTGTTTAATGCggcca                                                                                                                                                                                                                   >  2:476829/1‑72 (MQ=255)                                                                                                                                                                                                   gCCCCCAGTTCTTCTGCCAACTCATTGACTTCATCCCGCACGCTACCAAAACGGACAtgt                                                                                                                                                                                                                                  >  2:47320/1‑60 (MQ=255)                                                                                                                                                                                                    cccccAGTTCTTCTGCCAACTCATTTACTTCATCCCGCACGCTACCAAAACGGACATGTTGTTTAATGCGGGAAGGATCGATTGTGAAGTGGCTGACCATCGTTTGCAGACGTTCTTGTGCTTCATGTTGCAGATGCTCTTCAAAacgacg                                                                                                                                      <  2:621246/151‑1 (MQ=255)                                                                                                                                                                                                       ccAGTTCTTCTGCCAACTCATTTACTTCATCCCGCACGCTACCAAAACGGACATGTTGTTTAATGCGGGAAGGATCGATGGTGAAGTGGCTGACCATCGTTTGCAGACGTTCTTGTGCTTCATGTTGCAGATGCTCTTCAAAacgacg                                                                                                                                      <  1:364903/148‑1 (MQ=255)                                                                                                                                                                                                       ccAGTTCTTCTGCCAACTCATTGACTTCATCCCGCACGCTACCAAATCGGCCATGTTGTTTAATGCGGGAAGGAt                                                                                                                                                                                                               <  1:476837/75‑1 (MQ=255)                                                                                                                                                                                                       ccAGTTCTTCTGCCAACTCATTGACTTCATCCCGCACGCTACCAAAACGGACATgttgtt                                                                                                                                                                                                                              <  1:47320/60‑1 (MQ=255)                                                                                                                                                                                                             taagagacatcTCATTGACTTCATCCCGAACGCTACCAAAACGGACATGTTGTTTAATGCGGGAAGGATCGATGGTGATGTGGCTGACCATCGTTTGCAGACGTTCTTGTGCTTCATATTGCCGATGCGCTTCAAAACGACGCACATCAGCGGCAAAACGGTGCAGGCTCAGGCTGGCTGACCCGGGTAGTACGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCATTTCaaa  <  1:641018/264‑1 (MQ=255)                                                                                                                                                                                                              ttctGCCAACTCATTGACTTCATCCCGCACGCTACCAAAAc                                                                                                                                                                                                                                          >  2:372853/1‑41 (MQ=255)                                                                                                                                                                                                                  gCCAACTCATTGACTTCATCCCGCACGCTACCAAAACGGAc                                                                                                                                                                                                                                      <  1:372858/41‑1 (MQ=255)                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                  GGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGAACCTAACAGATGGGTCGAAATCGATGGGTTGCGAGAACCAATAACTACAACATCAGCCCCCAGTTCTTCTGCCAACTCATTGACTTCATCCCGCACGCTACCAAAACGGACATGTTGTTTAATGCGGGAAGGATCGATGGTGAAGTGGCTGACCATCGTTTGCAGACGTTCTTGTGCTTCATGTTGCAGATGCTCTTCAAAACGACGCACATCAGCGGCAAAACGGTGCAGGCTCAGGCTGGCTGACCCGGGTAGTACGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCATTTCAAA  >  CP009273/636813‑637290

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