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	96.5%	84.2 / NA	29	A33T (GCA→ACA)	rutG	pyrimidine:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (19/9);  minor base G (0/1);  total (19/10)
Fisher's exact test for biased strand distribution p-value = 3.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGCAATGCTGATATTCGGGTTAATGCCCTGACCGTTAAAGCCCGTCGCGGCGATCACCACGCCGACAAAAGCCGCGCTGGAGCCGAGATAACTGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCT  >  NC_000913/1069502‑1069979                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                             aGCGCAATGCTGATATTCGGGTTAATGCCCTGACCGTTAAAGCCCGTCGCGGCGATCACCACGCCGACAAAAGCCGCGCTGGAGCCGAGATAACTGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAACATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCCCCGCGTGTTGAACGCCCCTTACTGTCGGCTGAGc                                                                                                                                                                                                                                     >  1:378114/1‑251 (MQ=255)    gcAATGATTATATTCGGTTAAATCCCCTGCCCGTAAAACCCCGCCGCGGCAACCCCCCCGCCGACAAAAGCCGCGCTGGAGCCGAGATAACTGGGGCCCCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACAGTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCaaa                                                                                                                                                                                                                                  <  1:440823/251‑1 (MQ=255)    caattccTGATATTCGGTTTAACCCCCTGCCCGTTAAACCCCGCCGCGGCGACCCCCACGCCGACAAACGCCGCGCTGGAGCCAATAAACCGGGGGCCGCCCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCaaa                                                                                                                                                                                                                                  <  1:102333/245‑1 (MQ=255)              tatTCGGGTTAATGCCCTGACCGTTAAAGCCCGTCGCGGCGATCACCACGCCGACAAAAGCCGCGCTGGAGCCGAGATAACTGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGt                                                                                                                                                                                                                        >  2:320255/1‑251 (MQ=255)                      taaatgcctggACCTTTAAACCCCGCCGCGGCGATCCCCACGCAGAAAAAAGCCGCGCGGGGGCCAAGAAAACGGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCg                                                                                                                                                                                                                 <  1:114693/240‑1 (MQ=255)                               ggcccGTAAAAGCCCGTCGCGGCGACCCCCACGCCAAAAAAACCCGCGCTGGAGCCGAAAAAACGGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCAc                                                                                                                                                                                                       <  1:138350/248‑1 (MQ=255)                                    gTTAAAGCCCGTCGCGGCGATCACCACGCCGACAAAAGCCGCGCTGGAGCCGAGATAACTGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGAAAACGGCAATCGTTCATCCGGCGCCACCAcgccgc                                                                                                                                                                                                  >  2:199515/1‑251 (MQ=255)                                       aaaGCCCGACGCGGCGATCACCAAGACGACACAAGCCGCGCTGGAACCGAGATAAATGGGGACGCGCCCACCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATAAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGCGCAAACGGCAATCGTTCATCCGGCGCGCACACCCCGctct                                                                                                                                                                                               >  2:440034/1‑251 (MQ=255)                                                                    cAAAAGCCGCGCTGGAGCCGAGATAACTGGGGAAGCGACCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTg                                                                                                                                                                   >  1:248017/1‑250 (MQ=255)                                                                                     cgaaaacaccgggcgcctcccctgccgGTAAAAAAGAAAAACAGTAACGAGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTTTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGCCCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGc                                                                                                                                                 <  2:11462/228‑1 (MQ=255)                                                                                      cGAGATAACTGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGcc                                                                                                                                                >  1:299638/1‑251 (MQ=255)                                                                                           tAACTGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGAGCCCCGACATTAAAATGGAAAGGTTGGGATCCAGACCCATCAGAATCGGCATTAACACCGTCTCGCCAAACATCGAAACCGCGCGTTGAACGCCTCTTACTGTCGTCTGAGCAAACGGCACTCGTTCATCCGGCGCGACCACACCGCACTCTGAATAGGCCGATTTTCACTGCCAGTGAGGAAAAACGAAAATTACCAta                                                                                                                                             >  2:15223/1‑248 (MQ=255)                                                                                                           cgccgGTAATAAAGAAAAACAGTAACGTGCCGATCCACGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAagg                                                                                                                             >  1:44241/1‑249 (MQ=255)                                                                                                             ccgGTACTAAAGCAACACAGTAACGTGCAGATCCCCCACATTAAAATGGACAGATTAGGATCCAGCCCCATCAGAATCGGCATTCACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGCTTTTAACTGCCAGTGAGGAAAACCGCACATTGCCATCAGCTGTCTCCTTAAGGAGGt                                                                                                                         >  2:44866/1‑251 (MQ=255)                                                                                                                gTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTa                                                                                                                       >  1:70465/1‑250 (MQ=255)                                                                                                                              cAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGACCGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGACGATTTTAACTGCCAGTGAGGAAAACCGAACATAGCCATCAGCTGTCTCCTTAAGGAGGTTAACAATCAGGACGCAt                                                                                                        >  2:428757/1‑251 (MQ=255)                                                                                                                                gaaaCGTGCGGACCCCGGACATGAAAATGGAAAGATGGGGATCCACCCCCATCAGAATCGGCATTAACCCCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACGGGCGTCTGAGCAAACGGCAATCTTTTACCCGTCGCGACCACGCCGCTCTCGGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCa                                                                                                      <  1:44866/249‑1 (MQ=255)                                                                                                                                            tCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATAAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAAc                                                                                          >  1:239959/1‑251 (MQ=255)                                                                                                                                             ccccGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACt                                                                                         <  1:77189/251‑1 (MQ=255)                                                                                                                                                  aCATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTTAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGAt                                                                                    >  2:262740/1‑251 (MQ=255)                                                                                                                                                   caatAAAATGGAAAGCTTTGTACCCAGCCCCATCAGACCCTGCATTACCCCCGTCGCGCCAATCACCGCCCCCGCGTGTTGACCGCCCATTACTGTCGTCTGAGCAAACGGCAATCTTTCTCCCGGCGCGCCCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTTAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGTATCAGCTCTTGATAACTGCGATc                                                                                   <  2:44241/248‑1 (MQ=255)                                                                                                                                                     ttAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTAGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGACGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTACAAAGCAGGGAGAATCAGCGCGAGATAACTGCGAt                                                                                    >  2:222713/1‑248 (MQ=255)                                                                                                                                                            ggAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCacac                                                                          <  1:93537/251‑1 (MQ=255)                                                                                                                                                                            aGCCCCATCCGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAAAGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAACACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACACGCAGGGCGCATCCGCGCGTGATACCTGCGCTCGAACCACCCCAGCCCGTCGggcgtg                                                          >  2:173160/1‑251 (MQ=255)                                                                                                                                                                                              aTTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCt                                        >  1:89003/1‑251 (MQ=255)                                                                                                                                                                                              aTTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCt                                        >  1:342213/1‑251 (MQ=255)                                                                                                                                                                                                                            gcgTGTTGAACGCCCATTACTGTCGTCTGAGAAAACGGAAATCGTCCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTg          <  2:70465/251‑1 (MQ=255)                                                                                                                                                                                                                                 ttGAACGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGCACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAACACAGAATGTCGTGGGTGCCGAc     >  1:166524/1‑251 (MQ=255)                                                                                                                                                                                                                                    aaCGCCCATTACTGTCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGCGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGCGTGGGGTGCCGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGcc  >  1:290620/1‑250 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                             AGCGCAATGCTGATATTCGGGTTAATGCCCTGACCGTTAAAGCCCGTCGCGGCGATCACCACGCCGACAAAAGCCGCGCTGGAGCCGAGATAACTGGGGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCT  >  NC_000913/1069502‑1069979

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: ‑