Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,131,423	0	G	C	55.6%	15.1 / 26.0	18	A329G (GCG→GGG)	yghQ	putative inner membrane polysaccharide flippase
Reads supporting (aligned to +/- strand):  ref base G (5/3);  new base C (0/10);  total (5/13)
Fisher's exact test for biased strand distribution p-value = 6.54e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.11e-03

GCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCC‑TAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCG  >  NC_000913/3131273‑3131557                                                                                                                                                        |                                                                                                                                       gCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCC‑TAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAgcg                                                                                                                                        >  1:360118/1‑151 (MQ=255)               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)                                                                                                                                                        |                                                                                                                                       GCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCC‑TAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCG  >  NC_000913/3131273‑3131557

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