Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,575,956	0	A	C	53.3%	1.2 / 13.0	17	V369G (GTG→GGG)	gntU	gluconate transporter, low affinity GNT 1 system
Reads supporting (aligned to +/- strand):  ref base A (5/2);  new base C (0/8);  total (6/11)
Fisher's exact test for biased strand distribution p-value = 6.99e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.29e-03

TCGCGCCGGTAAATTTACCGAACAACCAGAAACCGGCGTCGTTAACGTGGCTGACAACAATCGAACCACCAGCGATACAAATC‑GACAGCGCCGCCATTTGCGCACCGGAGTAGTTCAGTTGTTCAATAACCGGCATCACCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTGATAGCAATCGGCAGGCCCATGCCGGTTAACGCTTCGCCCAGTGCCGGACCTACGCCAGAGTCAACCAG  >  NC_000913/3575819‑3576104                                                                                                                                           |                                                                                                                                                     tCGCGCCGGTAAATTTACCGAACAACCAGAAACCGGCGTCGTTAACGTGGCTGACAACAATCGAACCACCAGCGATACAAATC‑GACAGCGCCGCCATTTGCGCACCGGAGTAGTTCAGTTGTTCAATAACCGGCATCACCAGTCCTAccgc                                                                                                                                         <  1:447599/151‑1 (MQ=255)                   caaacaacCAAAACCCGGCGTCTTTACCGTGGCTAACAACAACCAACCCCCCAGCAATACAAATC‑GAAGGCCCCCCCATTTGCGCCCCGGAGTATTTCGGTTTTTCAATAACCGGCATCCCCAGTCCTACCGCCGTTAAACAGGCTACGGt                                                                                                                       <  1:433045/149‑1 (MQ=255)                         acCAGAAACCGGCGTCGTTAACGTGGCTGACAACAATCGAACCACCAGCGATACAAATC‑GACAGCGCCGCCATTTGCGCACCGGATTAGTTCAGTTGTTCAATAACCGGCATCACCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGa                                                                                                                 >  2:402331/1‑151 (MQ=255)                                  cGGCGTCTTAAGCGGGGCTGACAACAACCGACCCACCGGC‑AAAAAAATCAGACCGCGCCCCCATTGGCGCCCCGGAGTAGTTCAGTTTTTCAATAACCGGCATCCCCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAat                                                                                                        <  1:240159/151‑1 (MQ=255)                                     cggggtTAAGGGGGCGGACAACAATCGCACCACCAGCAAAAAAAATC‑GACCGCGCCGCCTTTTGCCCACGGGAATGTTTAAGTGGTTAAATAGCCGCCATCGCCCATCCTGCCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAatga                                                                                                      <  2:134049/146‑1 (MQ=255)                                           tAACGGGGCTGACAACACCCAACCCCCCAGCGATACAACCC‑GACAGCGCCGCCATTTGCGCACCGGAGTAGTTCAGTTTTTCAATACCCGGCATCACCAGTCCTCCCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCAc                                                                                               <  1:270630/151‑1 (MQ=255)                                                  gCTGACAACAATCGAACCACCAGCGATACAAATC‑GACAGCGCCGCCATTTGCGCACCGGAGTAGTTCAGTTGTTCAATAACCGGCACCACCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAg                                                                                          >  2:497220/1‑149 (MQ=255)                                                          caaTCGAACCACCAGCGATACAAATC‑GACAGCGCCGCCATTTGCGCACCGGAGTAGTTCAGTTGTTCAATAACCGGCATCACCAGTCCTACCGCCGTTAAACAGGCTACGGTTGCAGAACCCTCAATGATGCGCACTGCAGCTGCCAGCa                                                                                 >  2:565909/1‑150 (MQ=255)                                                                  ccccaagcaaaaaaaaCC‑GACGGCCCCCCCTTTGGCCCCCCGGAGTAGTTCAGTTGTTCAAAACCCGGCACCCCCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAgg                                                                        <  2:401961/138‑1 (MQ=255)                                                                    cccaGGGATAAAAATC‑GACAGCGCCGCCTTTTGCCCACCGGAGTAGTTCAGTTGTTCAATACCCGGCATCCCCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTg                                                                      <  2:278419/150‑1 (MQ=255)                                                                    cccaGCAATACAAATC‑GACGGCGCCGCCATTGGCGCACGGGAGGAGTGCAGTTGTACAATACCGGGCACCCCCAGCCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTg                                                                      <  1:140962/150‑1 (MQ=255)                                                                    acaaccGAAACAAACC‑GACGGCACCGCCATTGGCGCCCGGGGGTATTTAATTTGTTAAAAAACCGGCGTCCCCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTg                                                                      <  2:362636/146‑1 (MQ=255)                                                                                                                  gtTCAGTTGTTCAATAACCGGCATCTCCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTGATAGCAATCGGCAGGCCCATGCCGGTTAACGCTTCGCCCAGTGc                          >  1:207535/1‑150 (MQ=255)                                                                                                                       tttgttCAATAACCGCCATCCCCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTGATAGCAATCGGCAGGCCCATGCCGGTTAACGCTTCGCCCAGTGCCGGAcc                    <  2:482755/150‑1 (MQ=255)                                                                                                                        ttgtaCAAAACCCGCCATCCCCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTGATAGCAATCGGCAGGCCCATGCCGGTTAACGCTTCGCCCAGTGCCGGAcc                    <  1:450419/150‑1 (MQ=255)                                                                                                                                    cGGCATCACCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTGATAGCAATCGGCAGGCCCATGCCGGTTAACGCTTCGCCCAGTGCCGGACCTACGCCAGAGTCa       >  1:364439/1‑151 (MQ=255)                                                                                                                                          cACCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTGATAGCAATCGGCAGGCCCATGCCGGTTAACGCTTCGCCCAGTGCCGGACCTACGCCAGAGTCAACCAg  >  1:540476/1‑150 (MQ=255)                                                                                                                                           |                                                                                                                                                     TCGCGCCGGTAAATTTACCGAACAACCAGAAACCGGCGTCGTTAACGTGGCTGACAACAATCGAACCACCAGCGATACAAATC‑GACAGCGCCGCCATTTGCGCACCGGAGTAGTTCAGTTGTTCAATAACCGGCATCACCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGCAGGTGATAGCAATCGGCAGGCCCATGCCGGTTAACGCTTCGCCCAGTGCCGGACCTACGCCAGAGTCAACCAG  >  NC_000913/3575819‑3576104

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