Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,131,430	0	G	A	60.0%	5.3 / 21.7	20	L327L (CTG→TTG)	yghQ	putative inner membrane polysaccharide flippase
Reads supporting (aligned to +/- strand):  ref base G (6/2);  new base A (0/12);  total (6/14)
Fisher's exact test for biased strand distribution p-value = 7.22e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.77e-05

AAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCC‑TAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGAT  >  NC_000913/3131287‑3131580                                                                                                                                                 |                                                                                                                                                       aaatgaaatgACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCC‑TAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATAc                                                                                                                                                   <  1:118028/149‑1 (MQ=255)      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)                                                                                                                                                 |                                                                                                                                                       AAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCC‑TAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGAT  >  NC_000913/3131287‑3131580

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