Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,480,371	0	T	G	55.0%	15.8 / 26.1	20	G390G (GGT→GGG)	dsdA	D‑serine dehydratase
Reads supporting (aligned to +/- strand):  ref base T (3/6);  new base G (11/0);  total (14/6)
Fisher's exact test for biased strand distribution p-value = 2.17e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.83e-03

TTGGTCGCGCATCAGGCTTTGTCGGGCGGGCAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTATGG‑CCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAAC‑AGATGCAC‑GGTTTCAGCGCAGAACAACTGCGTAATACCACTCATCTGGTGTGGGCGACGGGAGGTGGAATGGTGCCGGAAGAAGAGATGAATCAATAT  >  NC_000913/2480226‑2480512                                                                                                                                                  |                                                                                                                                                 ttGGTCGCGCATCAGGCTTTGTCGGGCGGGCAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTAt                                                                                                                                                <  1:59243/148‑1 (MQ=255)               ggCTTTGTCGGGCGGGCAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGGTGG‑CCGGGCCACcgcgcg                                                                                                                              >  2:133927/1‑151 (MQ=255)                       cgggcgggCAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTATGG‑CCGGACCTCAGCGCGTGTGTGCa                                                                                                                      <  2:404347/151‑1 (MQ=255)                                            gctgGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTATGG‑CCGGGCCCCACCCCGGGGGGGGTTTCGTAAGTTACCAAC‑AAATg                                                                                                >  2:232335/1‑151 (MQ=255)                                                    gCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGGATTCGTCTTGAACCTTCGGCACTGGCGGGGATGG‑GCGGGCCTCAGCGCCTGGGGGCATCACTAAATTaacaacccaaacaagggtt                                                                                         >  1:224559/1‑138 (MQ=255)                                                      ttCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGATGG‑CCCGGGCCCCGGGCGGGTGTGGATCAGAAAAATACCAAC‑AAAGGGAC‑GGcttta                                                                                     >  2:83913/1‑149 (MQ=255)                                                       tCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGATGG‑GCGGAGCCCCGCGGCTGTGGGGGTCAGGAAAGTAACAACAAAATGAAC‑GGttttg                                                                                     >  1:79360/1‑149 (MQ=255)                                                       tCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGATGG‑GCGGACCTCAGCGCGTGGGTGCATCCGCAAGTTACCAAC‑AGaaggcagggtttag                                                                                     >  1:343475/1‑138 (MQ=255)                                                        cTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATGCGTCTTGAACCTTCGGCACTGGCGGGTATGG‑CCGGACCTCAGCGCGTGTGTGCa                                                                                                                      <  2:81379/118‑1 (MQ=255)                                                          ataCCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGATGG‑CCGGGCCCCCGGGCGGGTGGGCGGCAGGAAGTTAACAAC‑AGATGCAC‑GGgtttcgcgg                                                                                 >  2:9474/1‑145 (MQ=255)                                                                              aTGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGATGGCCCGGACC‑CCGCGCGGGTGGGCCTCATTAAAGTATCCAA‑AAAAGGAC‑GGGTTTTGCGCCGAAAAACCGCGGAATAcc                                                              >  2:410620/1‑150 (MQ=255)                                                                              aTGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGATGG‑GCGGGCCTCCGCGCGGGTGTGCATCAGTAAATTAACAAC‑AAATAGAC‑GGGTTCAGCGCCGAAAAACTGCGCAATAACa                                                             >  1:327014/1‑151 (MQ=255)                                                                                  aTGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGATGG‑CCGGAGCCCCACCCGTGTGTGCCTCAGTAAAGTATCAAC‑AAATGGAC‑GGGTTTAGCCCCGGACCACTGCGGAATACCAccc                                                          >  1:192084/1‑148 (MQ=255)                                                                                          cTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTATGG‑CCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAAC‑AGATGCAC‑GGTTTCAGCGCAGAACAACTGCGTAATACCACTCATCTGgtg                                                  <  1:98939/150‑1 (MQ=255)                                                                                          cTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTATGG‑CCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAAC‑AGATGCAC‑GGTTTCAGCGCAGAACAACTGCGTAATACCACTCATCTGgtg                                                  <  2:327137/150‑1 (MQ=255)                                                                                                           ggAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGAAGG‑GCGGACCCCCGCGCGGGTGTGCATCCGGAAGTTAACCAC‑AGAAGGACGGGTTTCA‑CGCCGGAAAAATGCGGAATACCACCCAACCGGTGGGGGGGGAGGGAGggggg                                >  1:73754/1‑150 (MQ=255)                                                                                                            gaagaaGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGATGG‑CCGGAACCCACCGCGTGGGTGCATCAGAAAGTTAACCAA‑AAATGGAC‑GGGTTCAGCGCAGAACAAATGCGGAATAACACCCCTCTTGTGTGGGGGGCGggggggggga                               >  2:333000/1‑149 (MQ=255)                                                                                                                          tCTTGAACCTTCGGCACTGGCGGGTATGG‑CCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAAC‑AGATGCAC‑GGTTTCAGCGCAGAACAACTGCGTAATACCACTCATCTGGTGTGGGCGACGGGAGGTGGAATGGTGCCGgaaga                  <  1:133896/150‑1 (MQ=255)                                                                                                                                  cTTCGGCACTGGCGGGTATTG‑CCGGACCCCAGCGACTGGGGGCATCAGTAAGTTTCCCAC‑AAAAGCAC‑GGGTTAAGCGCAGAACAACTTCGGCATACCCCTCATCTCGCGTGGGCGGCGGGGGGGGGAAAGGTGGCGGAAGCAAAGATGaa         >  1:163001/1‑151 (MQ=255)                                                                                                                                            ggCGGGTATGG‑CCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAAC‑AGATGCAC‑GGTTTCAGCGCAGAACAACTGCGTAATACCACTCATCTGGTGTGGGCGACGGGAGGTGGAATGGTGCCGGAAGAAGAGATGAATCAatat  >  2:216402/1‑148 (MQ=255)                                                                                                                                                  |                                                                                                                                                 TTGGTCGCGCATCAGGCTTTGTCGGGCGGGCAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTATGG‑CCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAAC‑AGATGCAC‑GGTTTCAGCGCAGAACAACTGCGTAATACCACTCATCTGGTGTGGGCGACGGGAGGTGGAATGGTGCCGGAAGAAGAGATGAATCAATAT  >  NC_000913/2480226‑2480512

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