Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,069,743	C→T	A33T (GCA→ACA)	rutG ←	pyrimidine:H(+) symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,069,743	0	C	T	97.1%	90.7 / NA	34	A33T (GCA→ACA)	rutG	pyrimidine:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (18/15);  total (19/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

GCCCAGCGCAATGCTGATATTCGGGTTAATGCCCTGACCGTTAAAGCCCGTCGCGGCGATCACCACGCCGACAAAAGCCGCGCTGGAGCCGAGATAACTGGGGACGCGCCCGCC‑GGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCA  >  NC_000913/1069498‑1069984                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  ccccAGCGCAATGCTGATATCCGGTTAAATGCCCGGCCCGTTAAACCCCGCCGGGGCGACCACCCCTCCGAAAAAAGCCGCCCTGGGGCGGAGAAAACTGGGGACGCGCCCGCC‑GGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCt                                                                                                                                                                                                                                              <  1:47781/250‑1 (MQ=255)              cTGATATGCGGGTTCATGCCCTGACCGTTCAAGCCCGTCGCGGCGATCACCACGCCGCCAAAAGCCGCGCTGGAGCCGAGATAACTGGGGACGCGCCCGCC‑GGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAACATGGAAAGATTGGGATCCACCCCCATCATAATCGGCATTAACACCGTCGCTCCGAACATCGCCCCCCCGAGTTGAACGCCCATTACTGCCTTCTGCGCAAACGGCa                                                                                                                                                                                                                                  >  2:125/1‑250 (MQ=255)                                                  gtcgccgagctaaagacgacGACAACAGCAGCGCTGGAGTCGAGATAACTGGGGCCGCGCCAGAC‑GGTAATCAGGAAAGCCAGTACCGTGTCGATACCAGACATTACGATGGAAAGATTGGGATCCAGCCCCATCAGAAACGGCCTTCACACCGTCGCGCCCAACATCGCCACCGCGTGTTGAGCGCCCCATACTGTCGTCTGAGCCACCGGCCATCGATCAGCGGGCACGACCCAGCCGCACTCTGtct                                                                                                                                                                                              >  2:163602/20‑248 (MQ=255)                                                     ggggcgACCACCCCCCCGAAAAAGGCCGCGCTGGGGCGGAAATAACGGGGGCCCCGCCCGCC‑GGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACCCCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGt                                                                                                                                                                                          <  1:333472/249‑1 (MQ=255)                                                                             gccgcgtgggagcggaaaaaagtGGGGCCGCGCCCGCC‑GGTAATAAAGAAAAACAGTACCGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGaaa                                                                                                                                                                  <  1:290099/229‑1 (MQ=255)                                                                                  gcTGGAGCCGAGATAACTGGGGACGCGCCCGCC‑GGTAATAAAGAAACACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCg                                                                                                                                                              >  1:177008/1‑250 (MQ=255)                                                                                    tGGAGCCGAGATAACTGGGGACGCGCCCGCC‑GGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAAc                                                                                                                                                           >  1:200232/1‑251 (MQ=255)                                                                                           gaaaaaaCTGGGGGCGCCCCCACCCGGTAATAAAGAAAAACAGTACCGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCTGCCCCATCAGAATCGGCATTAACACCTTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCAATTACTTTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGcc                                                                                                                                                     <  2:340850/246‑1 (MQ=255)                                                                                                 acggggaccgcccccgcc‑GGTAATAAAAAAAAACAGTAACTTGCCGACCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAACCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTCCTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCt                                                                                                                                              <  1:53125/239‑1 (MQ=255)                                                                                                    ggggACGCGCCCGCC‑GGTAATAAAGAAAAACAGTAACGTGCCGCTCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGtc                                                                                                                                           >  1:313625/1‑251 (MQ=255)                                                                                                      gaccgcgCCCGCC‑GGGAAAAAAGAAAAACAGTAACGTGCCGACCCCCGACATAAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGGGTTGAACGCCAATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGTCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGtctc                                                                                                                                         <  1:225688/248‑1 (MQ=255)                                                                                                             cccgcc‑GGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAACTGGCAAGATCGGGATCAAGCCACCTCAGAATCGGCATTAACACCGTCGCGCCAAAAATCACCACCACACGTTGAAAGCCCATTAATGTCGTCTGAGCAAACGGCACTCGTTCATCAGGCAAAACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCCTCAGCTGTCTCCTTAagg                                                                                                                                  >  1:56811/1‑251 (MQ=255)                                                                                                                        ataaAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACa                                                                                                                         >  1:272102/1‑250 (MQ=255)                                                                                                                              aaaaaCAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTACGGAGGTTAACAAGCAggg                                                                                                                  >  1:247718/1‑251 (MQ=255)                                                                                                                                       aacgtctcGATCCCCGACATGAAAAGTGAAAGATTGGGATCCAGCCCCATCAGAAGCGGTATTAACACCGTCTCTCCAAACATCGCCACCGCGTGTTGATCGCCCATTACTTTCGTCTGATCAAACTGCAATCGGTCATCCGTCGCGACCACGCCGATCTCTGTTGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGACCATTGCCAACAGCTGTAGCCTTAAGGAAGTTAACAAGCAGGGCGCATCAgc                                                                                                         <  2:254290/244‑1 (MQ=255)                                                                                                                                           tGCCGCTCCCCGCCATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACCCCCATTACTGTCGTCTGCGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGCTTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGt                                                                                                     >  2:390842/1‑251 (MQ=255)                                                                                                                                                         aTTAACATGGAACGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGGGGAAAACCGAACATTGCCATCAGCTGTCTACTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCg                                                                                       >  2:181050/1‑251 (MQ=255)                                                                                                                                                          ttAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGCACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGACGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCg                                                                                       >  1:382080/1‑250 (MQ=255)                                                                                                                                                                    aaGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACAcc                                                                              >  1:331246/1‑249 (MQ=255)                                                                                                                                                                       aTTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGccc                                                                         >  1:27988/1‑251 (MQ=255)                                                                                                                                                                                 agcccactCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTACCTGCCAGTGAGGAAAACCGCACATTGCCATAAGCTGTATCCTTAAGGAGGGTAAAAAGCAGGGCGCATCAGCGCGTGATAACTGAGATCGAACCACACCAGCCCGCAGggtgcg                                                               >  2:36887/8‑249 (MQ=255)                                                                                                                                                                                  cccccACCAGAATCGGCATTACAACCGTCGCGCCAAACACCGCCACGGGGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGggtgtgg                                                              <  2:152597/250‑1 (MQ=255)                                                                                                                                                                                   ccccATCCGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGggtgtggg                                                             >  1:288920/1‑250 (MQ=255)                                                                                                                                                                                              ccGGCATTAACCCCGTCGCGCCAAACATCGCCACCGCTTGTGGAACGCCCATTACTGTCGTCGGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGaa                                                  <  2:200232/250‑1 (MQ=255)                                                                                                                                                                                                 gCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCg                                               >  1:202487/1‑251 (MQ=255)                                                                                                                                                                                                  cATTAACCCCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCCTCAGCTGTCTCCTTAAGGAGGTTAACAAGCCGGGCGCATCAGCGCGTGATAACTGCGCTCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGc                                              >  1:51116/1‑251 (MQ=255)                                                                                                                                                                                                     aaaccccttCGCGCCAAACATCGCCCCCGCGTGTTGAACGCCCATTCCTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTc                                           <  1:9805/243‑1 (MQ=255)                                                                                                                                                                                                            gTCGCGCGAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGACTGAGCAAACGGCAATCGTTCATCCGGCGCGACCAAGACGCTCTCTGAAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGGCGATGAATCGCCTCGCTGgcg                                    >  2:387771/1‑251 (MQ=255)                                                                                                                                                                                                                           cccACCTCGTGTTGAACGCCAATTACGGTCGTCTGAGCAAACGGCAATCGTTCACCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCg                     <  1:346935/250‑1 (MQ=255)                                                                                                                                                                                                                             cACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAACCGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGACGGCGCAAAACAGAATGTCGTg                   <  1:277677/251‑1 (MQ=255)                                                                                                                                                                                                                                           cGCCCATTACGGTCGTCTGAGCAAACGGCAATCGTCCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTca     <  2:309074/251‑1 (MQ=255)                                                                                                                                                                                                                                              ccATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGa                                                         >  2:271572/1‑196 (MQ=255)                                                                                                                                                                                                                                              ccATTACTGTCGTCTGAGCAAACGGCAATCGTTCAGCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCGCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGa                                                         <  1:271572/196‑1 (MQ=255)                                                                                                                                                                                                                                              ccATTACTGTCGTCTGACCAAACGGCAATCGTTCATCCGGCGCTCCCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCTAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTcacca  <  1:387771/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  GCCCAGCGCAATGCTGATATTCGGGTTAATGCCCTGACCGTTAAAGCCCGTCGCGGCGATCACCACGCCGACAAAAGCCGCGCTGGAGCCGAGATAACTGGGGACGCGCCCGCC‑GGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCA  >  NC_000913/1069498‑1069984

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