Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,489,174	0	G	T	39.6%	64.7 / 18.3	48	G315C (GGT→TGT)	aldA	aldehyde dehydrogenase A, NAD‑linked
Reads supporting (aligned to +/- strand):  ref base G (6/23);  new base T (19/0);  total (25/23)
Fisher's exact test for biased strand distribution p-value = 1.61e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

TCTTGAACTGGCAGTCAAAGCCATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCA‑ATTTGGTAACCCCGCTGAACGCAACGACA‑TTGCGATGGGG‑CCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAG‑CAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCGGTGGCAAAGCGGTAGAGGGGAAAGGATATT  >  NC_000913/1489026‑1489321                                                                                                                                                      |                                                                                                                                                       tCTTGAACTGGCAGTCAAAGCCATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCA‑ATTTTGt                                                                                                                                                      >  1:90046/1‑151 (MQ=255)    tGAACTGGCAGTCAAAGCCATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCA‑ATTTg                                                                                                                                                        <  1:221929/146‑1 (MQ=255)    tGAACTGGCAGTCAAAGCCATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCA‑ATTTGGTaa                                                                                                                                                    <  1:977264/150‑1 (MQ=255)     gAACTGGCAGTCAAAGCCATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTTT‑ATTTTGTaa                                                                                                                                                    >  2:861806/1‑149 (MQ=255)      aaCTGGCAGTCAAAGCCATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCA‑ATTTGGTAAccc                                                                                                                                                 <  1:32462/151‑1 (MQ=255)       aCTGGCAGTCAAAGCCATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCA‑ATTTGGTAAcccc                                                                                                                                                <  1:909713/151‑1 (MQ=255)          ggCAGTCAAAGCCATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCA‑ATTTGGTAACCCCGCt                                                                                                                                             <  1:1013602/151‑1 (MQ=255)             aGTCAAAGCCATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTTA‑ATTTTGTAAcccc                                                                                                                                                >  1:213000/1‑145 (MQ=255)             aGTCAAAGCCATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCA‑ATTTGGTAACCCCGCTGaa                                                                                                                                          <  1:445672/151‑1 (MQ=255)                 aaaGCCATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCC‑ATTTTTTAACCCCCCCGAAc                                                                                                                                         >  2:492670/1‑148 (MQ=255)                  aaGCCATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCA‑ATTTGGTAACCCCGCTGAACGCaa                                                                                                                                     <  1:34656/151‑1 (MQ=255)                         cGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCA‑ATTTTGTAACCCCGCTGAACGCAAc                                                                                                                                    >  1:690586/1‑145 (MQ=255)                            tGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCA‑ATTTGGTAACCCCGCTGAACGCAACGACA‑TTGc                                                                                                                           <  1:492742/150‑1 (MQ=255)                                   cgcgTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCA‑ATTTGGTAACCCCGCTGAACGCAACGACA‑TTGCGATGGGG‑c                                                                                                                  <  1:998686/151‑1 (MQ=255)                                          ttAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCC‑ATTTTGTAACCCCGCCGAACGCAACGACA‑TTGCGGTGGGG‑CCGTTGAt                                                                                                           >  1:75557/1‑151 (MQ=255)                                               aGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCACTTTTGGGAACCCC‑CTGAAAGCCACGAAC‑TTGCCATTGGG‑CCGTTGttt                                                                                                          >  1:1006401/1‑147 (MQ=255)                                                gTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCA‑ATTTTGTAACCCCCCTGGACGCAACGACA‑TTTCGGTGGGG‑CCGTTTATTAAagc                                                                                                     >  2:217787/1‑151 (MQ=255)                                                 tGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTTA‑ATTTTGTAAACCCCCTGAACGCAACGACA‑ATGGGGTGGGG‑CGCTTGTTTaa                                                                                                        >  2:11998/1‑147 (MQ=255)                                                  gggCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCA‑ATTTGGTAACCCCGCTGAACGCAACGACA‑TTGCGATGGGG‑CCGTTGATTAAcgcc                                                                                                    <  2:475385/150‑1 (MQ=255)                                                     cAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCA‑ATTTGGTAACCCCGCTGAACGCAACGACA‑TTGCGATGGGG‑CCGTTGATTAACGCCgcgg                                                                                                <  2:578647/151‑1 (MQ=255)                                                      aaGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTTT‑AATTTGTAAACCCGCCGGACCCCACGGCA‑ATGAGACGGGG‑GCGTTTGTTTACGCCgcgg                                                                                                >  2:465223/1‑150 (MQ=255)                                                      aaGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCA‑ATTTGGTAACCCCGCTGAACGCAACGACA‑TTGCGATGGGG‑CCGTTGATTAACGCCgcggc                                                                                               <  2:100631/151‑1 (MQ=255)                                                          gtgtAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCA‑ATTTGGTAACCCCGCTGAACGCAACGACA‑TTGCGATGGGG‑CCGTTGATTAACGCCGCGGCGCTg                                                                                           <  2:756807/151‑1 (MQ=255)                                                            gtAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCA‑ATTTTGTAAACCCCCTCAACGCAACGACA‑TTTCGATGGGG‑CCGGTGGTTTACGCCccgg                                                                                                >  2:566748/1‑144 (MQ=255)                                                                  gtgCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCA‑ATTTGGTAACCCCGCTGAACGCAACGACA‑TTGCGATGGGG‑CCGTTGATTAACGCCGCGGCGCTGGAAAggg                                                                                    <  1:735034/150‑1 (MQ=255)                                                                     cAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCA‑ATTTGGTAACCCCGCTGAACGCAACGACA‑TTGCGATGGGG‑CCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGa                                                                                <  2:796370/151‑1 (MQ=255)                                                                         aCGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCA‑ATTTGGTAACCCCGCTGAACGCAACGACA‑TTGCGATGGGG‑CCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAG‑Caa                                                                           <  1:181523/151‑1 (MQ=255)                                                                           gtgtTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCA‑ATTTTGTAACCCCGCTGAACGCAACGACA‑TTGCGATGGGG‑CCGGTGATTAACGCCGCGGCGCTGGAAAGGGGCGAG‑Caaa                                                                          >  1:505225/1‑150 (MQ=255)                                                                           gtgtTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCA‑ATTTTGTAAACCCCCTGAACGCAACGACA‑TTGGGATGGGG‑GCGGTTATTAACCCCCCGGCGCGGGGAAGGGGCGAg                                                                               >  2:497621/1‑146 (MQ=255)                                                                              tttATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTTA‑ATTTTGTAACCCCCCTGAACGCAACGACA‑TTTCGATGGGG‑CCGTTGATTAACGCCGCGGGGCTGGGAAAGGGCGGG‑Caaaaa                                                                        >  2:824310/1‑149 (MQ=255)                                                                               ttATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCA‑ATTTTGTAACCCCGCTGAACGCAACGACA‑ATGCGGTGGGG‑CCGTTGATTAACGCCGCGGGGCTGGAAAGGGTCGAG‑Caaaaa                                                                        >  1:337021/1‑148 (MQ=255)                                                                                tATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTTA‑AATTGGTAACCCCCCCGGACGCAACGACA‑TTGGGGGGGGG‑CCGTTTATTAAACCCGCCGGGCGGGGAAGGGGCGAG‑Caaaaa                                                                        >  2:347407/1‑147 (MQ=255)                                                                                         aaaGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTTA‑ATTTTGGAACCCCCGCGGACGCCACAACA‑TTGTGGAGGGG‑CCGTTGTTTAACGCCGCGGGGGTGGAAAAGGGCGAG‑CAAAAAATGgcgcgcgcgg                                                           >  1:841489/1‑149 (MQ=255)                                                                                         aaaGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCA‑ATTTTGTAACCCCCCTGAACGCAACGACA‑ATTCGATGGGG‑GCGTTGATTTACGCCGCGGCGCTGGGAAGGGGCGAG‑CAAAAAAAGgcg                                                                  >  2:67563/1‑144 (MQ=255)                                                                                             gCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCA‑ATTTGGTAACCCCGCTGAACGCAACGACA‑TTGCGATTGGG‑CCGTTGATTAACGCCCCGGCGCTGGGAAGGGGCGCG‑CAAAAA‑TGGCGCGCGCAg                                                           >  1:799199/1‑146 (MQ=255)                                                                                               aTTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCA‑ATTTGGTAACCCCGCTGAACGCAACGACA‑TTGCGATGGGG‑CCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAG‑CAAAAAGTGGCGCGCGCAGTagaag                                                     <  2:668062/151‑1 (MQ=255)                                                                                                  tATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGCTTT‑ATTTGGGAACCCCCCCCAACGCAAC‑ACACTTGCGGTGGGGTCCGTT‑ATTAACAGCACCGCGCTCGGAAAGGTCGAG‑CCAAAACAGGCCCGCGCAGt                                                          >  1:296232/1‑143 (MQ=255)                                                                                                   aTGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTTA‑AATTGGTAAACCCGCTGGACGCAACCACC‑TTGCGATGGGG‑GCGGTTGTTTACACCGCCGCGGCTGAAAGGGGGGGC‑CAAAAAATGgggcgcgccccgaaaaa                                                    >  1:788991/1‑139 (MQ=255)                                                                                                   aTGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGGTTA‑AATTTGTAACCCCCCCGGACGCAAAGCCC‑TTTCGATGGGG‑GCGGTGTTTTAAACCGCGGCGCTGGGGAGGGGCGGG‑CGAAAAACTGCGCGCGCCGGagaaaga                                                   >  2:132917/1‑146 (MQ=255)                                                                                                     gATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCA‑ATTTTGTAAACCCGCCGAAATCAACGACA‑TTTCCGTGGGG‑CCGTTGTTTATCGCCGCGGCGGCGGGAAGGGTGGAGCCAAAAAGT‑GTGGGGGCCGGAGAAGAAgggg                                               >  2:886901/1‑151 (MQ=255)                                                                                                                cAATCGGCTGGGTGAAGCGATGCAGGCGGTTTA‑ATTTTGTAACCCCCCCGAAAGCAACGACA‑TTGCGATGGGG‑GCGTTTATTAACACCGCGGCGCTGGGAAGGGGCGAG‑CAAAAAATGGCGCGCGCAGTAGAAGAAGGGGCGAGAGGGGCg                                    >  1:96740/1‑151 (MQ=255)                                                                                                                         gggTGAAGCGATGCAGGCGGTTCA‑ATTTGGTAACCCCGCTGAACGCAACGACA‑TTGCGATGGGG‑CCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAG‑CAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCGGTGGc                           <  2:420833/151‑1 (MQ=255)                                                                                                                              aaGCGATGCAGGCGGTTCA‑ATTTGGTAACCCCGCTGAACGCAACGACA‑TTGCGATGGGG‑CCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAG‑CAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCGGTGGCAAAGc                      <  1:426936/151‑1 (MQ=255)                                                                                                                                      cAGGCGGTTCA‑ATTTGGTAACCCCGCTGAACGCAACGACA‑TTGCGATGGGG‑CCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAG‑CAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCGGTGGCAAAGCGGTAGAgg              <  1:94201/151‑1 (MQ=255)                                                                                                                                      cAGGCGGTTCA‑ATTTGGTAACCCCGCTGAACGCAACGACA‑TTGCGATGGGG‑CCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAG‑CAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCGGTGGCAAAGCGGTAGAgg              <  1:492688/151‑1 (MQ=255)                                                                                                                                       aGGCGGTTCA‑ATTTGGTAACCCCGCTGAACGCAACGACA‑TTGCGATGGGG‑CCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAG‑CAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCGGTGGCAAAGCGGTaga                <  2:1014059/148‑1 (MQ=255)                                                                                                                                        ggCGGTTCA‑ATTTGGTAACCCCGCTGAACGCAACGACA‑TTGCGATGGGG‑CCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAG‑CAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCGGTGGCAAAGCGGTAGAggg             <  1:493351/150‑1 (MQ=255)                                                                                                                                           ggTTCA‑ATTTGGTAACCCCGCTGAACGCAACGACA‑TTGCGATGGGG‑CCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAG‑CAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCGGTGGCAAAGCGGTAGAGGGGaaa         <  1:348448/151‑1 (MQ=255)                                                                                                                                                   tttGGTAACCCCGCTGAACGCAACGACA‑TTGCGATGGGG‑CCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAG‑CAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCGGTGGCAAAGCGGTAGAGGGGAAAGGAtatt  >  1:962597/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                       TCTTGAACTGGCAGTCAAAGCCATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCA‑ATTTGGTAACCCCGCTGAACGCAACGACA‑TTGCGATGGGG‑CCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAG‑CAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCGGTGGCAAAGCGGTAGAGGGGAAAGGATATT  >  NC_000913/1489026‑1489321

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