Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,131,440	0	T	C	60.0%	9.0 / 21.5	20	G323G (GGA→GGG)	yghQ	putative inner membrane polysaccharide flippase
Reads supporting (aligned to +/- strand):  ref base T (7/1);  new base C (0/12);  total (7/13)
Fisher's exact test for biased strand distribution p-value = 1.03e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.77e-05

AAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCC‑TAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGAT  >  NC_000913/3131292‑3131580                                                                                                                                                      |                                                                                                                                             aaatgaCTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCC‑TAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTc                                                                                                                                             <  2:625999/150‑1 (MQ=255)                     aaaaaattCCCCCCTTTTTCCCCCCCTCCCAAAAACCAGTTCATCCCCTTGGAAATCCTAAACCCAAACCCCCAGAGAAAGGGGGGTTTTGCAGACGACCACCCCCGCGAGCCCCACAAAAATACCGTTCCCGCCCGCCAGTAAACCGGAc                                                                                                                         <  1:503515/145‑1 (MQ=255)                                     ttcccccccTCCCC‑AAAATCAGGTCATCCGCTTCGAGATACTTACCGCCAAACCCCAGCAAAAGGAGCGGCTTGCCGACGATCACCCCCGCGAGCCCCACAAAAATCCCGATCCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCa                                                                                                        <  1:423248/144‑1 (MQ=255)                                                ccc‑TAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGTCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCtggtgg                                                                                             >  1:215343/1‑151 (MQ=255)                                                ccc‑AAAAGCAGGCCATACGCTTCGAGATACTAAACCCCAAACCCCAGCAAAAGGAGCGGCTTGCCGACGAACAACACCGCGAGCCCCACAAAAATACCGATCCCCCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCtggtgg                                                                                             <  1:578404/151‑1 (MQ=255)                                                ccc‑AAAACCAGTCCATCCGTTCCGAGATACTAACCCCCAAACACCAGCAAAAGGGGCGGCTGCCCGACGACCACCCCCGCGAGCCCCACAAAAACCCCGATCCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCtggtgg                                                                                             <  2:43104/151‑1 (MQ=255)                                                      aaTCGGGTCACACGCTGCGACATATTTACCCCCAAACCCCAGCAAAGGGGGCGGCCCCCCCACGACCAACCCCGCGCGCCCCACCAAAACCCCGTTCCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTgcgc                                                                                        <  2:566995/151‑1 (MQ=255)                                                                tccGCTTCGAAATATTAAACCCAAAACCCCAGCAAAATGGGCGCCTCGCCAACAACCAACCCCGCGCGCCCCACCAAAATCCCGTTCCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAAc                                                                              <  1:345533/149‑1 (MQ=255)                                                                  cccTTAGAGATACTTACCGCAAAACCCCGGCAAAATGAGCGGCTTCCCGACGATCACCCCCGCGGGCCCCACCAAAATCCCGATCCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGc                                                                            <  1:596067/149‑1 (MQ=255)                                                                          aaTACTTACCGCCAAACCCCAGCGAAAGGAGCGCCTTGCCGACGACAAACCCCGCGGGCCCCACCAAAATCCCGATCCCCCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACctc                                                                    <  1:90884/150‑1 (MQ=255)                                                                                    gCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTaaaaa                                                          >  2:198284/1‑151 (MQ=255)                                                                                                caaaaGGAGCGGCTTGCCGACGATCAACCCCGCGCGCGCCACCAAAATCCCGATCCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTAccc                                                 <  1:318748/146‑1 (MQ=255)                                                                                                                    cgaacacccccGCGAGCCCCACCAAAATCCCGATCCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCcggc                          <  1:357691/142‑1 (MQ=255)                                                                                                                     gaccAACCCCGCGAGCCCCACAAAAATCCCGATCCCGCCCGCCAGTAAACCGGACTTCACACCTaaca                                                                                                            <  2:459973/65‑1 (MQ=255)                                                                                                                     gaTCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTaaca                                                                                                            >  1:459982/1‑68 (MQ=255)                                                                                                                     gaTCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCcggcg                         >  2:543897/1‑151 (MQ=255)                                                                                                                     gaTCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCcggcg                         >  2:595260/1‑151 (MQ=255)                                                                                                                                   cccccacaaAAATCCCGTTCCCCCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGt           <  2:294014/148‑1 (MQ=255)                                                                                                                                           aGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGa   >  1:641922/1‑151 (MQ=255)                                                                                                                                            gAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGAt  >  2:465581/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                             AAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCC‑TAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGAT  >  NC_000913/3131292‑3131580

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