Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,489,174	0	G	T	50.0%	63.4 / 32.8	54	G315C (GGT→TGT)	aldA	aldehyde dehydrogenase A, NAD‑linked
Reads supporting (aligned to +/- strand):  ref base G (6/21);  new base T (27/0);  total (33/21)
Fisher's exact test for biased strand distribution p-value = 1.14e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAACTGGCAGTCAAAGCCATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACG‑ACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCGGTGGCAAAGCGGTAGAGGGGAAAGGAT  >  NC_000913/1489029‑1489318                                                                                                                                                  |                                                                                                                                                  tGAACTGGCAGTCAAAGCCATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTaa                                                                                                                                               <  2:803298/150‑1 (MQ=255)  gAACTGGCAGTCAAAGCCATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGGTCAATTTTGt                                                                                                                                                 >  1:103733/1‑147 (MQ=255)   aaCTGGCAGTCAAAGCCATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAAccc                                                                                                                                            <  2:930048/151‑1 (MQ=255)       ggCAGTCAAAGCCATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGc                                                                                                                                         <  2:128300/150‑1 (MQ=255)       ggCAGTCAAAGCCATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCt                                                                                                                                        <  2:8955/151‑1 (MQ=255)          aGTCAAAGCCATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTTGTAACCCCGCTGa                                                                                                                                      >  2:1007363/1‑150 (MQ=255)            tCAAAGCCATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTTGTAACCCCGCTGAACg                                                                                                                                   >  2:717788/1‑151 (MQ=255)              aaaGCCATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTTGTAAACCCGCTGaa                                                                                                                                     >  2:1210233/1‑147 (MQ=255)                 gCCATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTTGTAACCCCCCTGAACGCaa                                                                                                                                >  2:536925/1‑149 (MQ=255)                   cATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTTGTAACCCCGCTGAACACaa                                                                                                                                >  1:773710/1‑147 (MQ=255)                     tCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTTGGAACCCCGCTGAACg                                                                                                                                   >  1:772814/1‑142 (MQ=255)                           aTTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTTCATTTTGTAACCCCGCTGAACGCAAAG‑ACATTGCg                                                                                                                     >  1:1203196/1‑149 (MQ=255)                             tCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACG‑ACATTGCGATg                                                                                                                  <  2:940284/150‑1 (MQ=255)                               aCGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTTGTAAACCCCCCGAACGCAACG‑ACATTGCGATgggg                                                                                                               >  1:1256823/1‑151 (MQ=255)                                     cATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACG‑ACATTGCGATGGGGCCGTTg                                                                                                         <  2:293549/151‑1 (MQ=255)                                            aGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTTGTAACCCCCCTGAACGCAACG‑ACATTGCGATGGGGGCGTTTAtt                                                                                                      >  1:20161/1‑147 (MQ=255)                                            aGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACG‑ACATTGCGATGGGGCCGTTGATTAAc                                                                                                   <  2:177702/150‑1 (MQ=255)                                              tGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAAATTTGTTACCCCGCTGAACACCACC‑ACAATTCGCTGGGGCCCGTGATTAAcgcc                                                                                                >  1:111329/1‑151 (MQ=255)                                               gggCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACG‑ACATTGCGATGGGGCCGTTGATTAAcgcc                                                                                                <  1:660791/150‑1 (MQ=255)                                                ggCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACG‑ACATTGCGATGGGGCCGTTGATTAAcgccgc                                                                                              <  2:1128846/151‑1 (MQ=255)                                                      tgtgtAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTTGTAACCCCCCTGAACGCAACG‑ACATTGCGGTGGGGCCGTTGGTTTAAACCccggcg                                                                                          >  2:943309/1‑149 (MQ=255)                                                       gtgtAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAAACCCCCTGAACGCAACG‑ACATTGCGGTGGGGCCGTTGATTAAAGCCccggcg                                                                                          >  1:754068/1‑148 (MQ=255)                                                        tgtAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTTGTAACCCCCCTGAACGCAACG‑ACATTGCGATGGGGGCGGTGATTTAAGCCGCGGCGCTgg                                                                                      >  1:1205389/1‑151 (MQ=255)                                                           aaCTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTTCATTTTGTAACCCCGCTGAACGGAACG‑ACATTTGGATGGGGCCGTTGGTTAACGCCGCGGCGCTGGaa                                                                                    >  2:98272/1‑150 (MQ=255)                                                           aaCTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTTAATTTTGGAACCCCCCCGAACGCAACG‑ACATTTCGATGGGGGCGTTGGTTAACGCCGCGGCGCTggga                                                                                    >  1:95070/1‑148 (MQ=255)                                                            aCTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTTAATTTTGTAAACCCCCTGAACGCCACG‑ACATTGCGGTGGGGCCGGTGGTTAACGCCccgg                                                                                            >  1:612115/1‑141 (MQ=255)                                                               gtgCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACG‑ACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAgg                                                                                 <  2:1190194/149‑1 (MQ=255)                                                               gtgCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACG‑ACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGt                                                                               <  2:110007/151‑1 (MQ=255)                                                                     aaCGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACG‑ACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCa                                                                         <  2:891288/151‑1 (MQ=255)                                                                     aaCGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACG‑ACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCa                                                                         <  1:514540/151‑1 (MQ=255)                                                                       cGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTTGTAACCCCCCTGAACGCCACG‑ACATTGCGATGGGGCCGTTGATTAACACCGCCGCGCTGGAAAGGGTCGGGCaa                                                                        >  2:653415/1‑150 (MQ=255)                                                                        gtgtTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTTAATTTTGTAACCCCCCCGAAAGCAACG‑ACATTGCGATGGGGCCGTTTGTTTACGCCGCGGCGCTGGGAAAGGGCGAGCaaa                                                                       >  2:1214652/1‑150 (MQ=255)                                                                         tgtTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTTGTAACCCCCCCTAACGCAACG‑ACCTTGCGATGGGGCCGTTTATTAAAGCCGCGGCGCTGGGAAGGGGCGAGCaaaa                                                                      >  2:773617/1‑150 (MQ=255)                                                                           tttATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGGAACCCCGCTGAACGCAACG‑ACATTGCGGTGGGGGCGGTGATTTAAGCCGCGGCGCTGGGAAGGGGCGAGCaaaaa                                                                     >  1:992153/1‑149 (MQ=255)                                                                                gTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACG‑ACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCAAAAAGTGgcg                                                               <  1:1005906/150‑1 (MQ=255)                                                                                           cATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTTGTAACCCCCCTGGACGCAACG‑ACATTTCGATGGGGCCGTTGGTTAAAGCCCCGGCGCTGGGAAAGGGCGGGCCAAAAAGGgcgcgcgccgga                                                      >  1:1160485/1‑144 (MQ=255)                                                                                             tttATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTTAATTTTGTAACCCCCCCGAACCCCACGAACATT‑CGGTGGGGGCGGTGTTTTAAGCCGCGGCGCTGGCAAGGGTGGGGCaaaaaagggggcg                                                             >  1:1079176/1‑131 (MQ=255)                                                                                                    tCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTTGTAACCCCCCCGGACGCAACG‑GCATTGCGATGGGGGCGTTTGTTTACGCCCCGGCGCTGGGAAAGGTCGGGCCAAAAGGGGCGCCGGCCGGAGAAGAAgggg                                            >  1:564095/1‑149 (MQ=255)                                                                                                        ttCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTTAAATTTGTAAACCCGCTGAAAGCAAAG‑GCATTTCGATGGGGCCGGTGATTAACGCCGCGGCGCCGGGAAAGGGCGGGCCAAAAATGGCGCGCGCCGGAGAAGAAAGGGCgag                                        >  2:65073/1‑149 (MQ=255)                                                                                                            tCAATCGGCTGGGTGAAGCGATGCAGGCGGGTCAATTTTGTAAACCCCCTGAACGCAAAC‑ACATTTCGGTGGGGCCGTTGATTAAAGCCGCCGCGCTGGGAAGGGTCGAGCAAAAAAAGGCGCGCGCCGTAGAAGAAGGGGGgagggggg                                   >  2:42110/1‑145 (MQ=255)                                                                                                                  ggCTGGGTGAAGCGATGCAGGCGGTTCAATTTTGTAACCCCGCTGAACGCAACG‑ACAATGCGGAGGGGCCGTTGTTTAACGCCGCGGCGCTGGAAAGGGGCGCGGCAAAAGTGGTGCGCGCCGTAAAAGAAAGGGGGggag                                      >  1:826814/1‑141 (MQ=255)                                                                                                                  ggCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACG‑ACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCgg                            <  2:490828/151‑1 (MQ=255)                                                                                                                        gTGAAGCGATGCAGGCGGTTCCATTTTGTAAACCCCCTGGAAGCAACG‑CCATTGCGATGGGGGCGCTGGTTAACACCCCCGCGCTGGAAAGGGTCGGGCAAAAAAAGGCGCGCGCCGGAGAAGAAGGGGGGAGGGGGGGGGTCGGTgg                         >  1:25665/1‑148 (MQ=255)                                                                                                                        gTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACG‑ACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCGGTGGCaa                      <  1:911664/151‑1 (MQ=255)                                                                                                                              cGATGCAGGCGGTTCCATTTTGTAAACCCCCTGAACGCAACG‑ACATTGCGATGGGGCCGTTGATTAACGCCCCGGCGCTGGGAAGGGTCGAGCAAAAAGTGGCGCGCGCAGGAGAAGAAGGGGGGAGGGGGGCGGTCGGGGGGAAAGCggg                >  1:49481/1‑150 (MQ=255)                                                                                                                                  gCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACG‑ACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCGGTGGCAAAGCGGTagag            <  1:569181/151‑1 (MQ=255)                                                                                                                                    aGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACG‑ACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCGGTGGCAAAGCGGTAGAggg          <  1:236333/151‑1 (MQ=255)                                                                                                                                     ggCGGTTCAATTTGGTAACCCCGCTGAACGCAACG‑ACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCGGTGGCAAAGCGGTAGAgggg         <  1:337866/151‑1 (MQ=255)                                                                                                                                      gCGGTTCAATTTGGTAACCCCGCTGAACGCAACG‑ACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCGGTGGCAAAGCGGTAGAGGGGa        <  2:372529/151‑1 (MQ=255)                                                                                                                                      gCGGTTCAATTTGGTAACCCCGCTGAACGCAACG‑ACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCGGTGGCAAAGCGGTAGAGGGGa        <  2:224808/151‑1 (MQ=255)                                                                                                                                       cGGTTCAATTTGGTAACCCCGCTGAACGCAACG‑ACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCGGTGGCAAAGCGGTAGAGGGGa        >  1:567418/1‑150 (MQ=255)                                                                                                                                       cGGTTCAATTTGGTAACCCCGCTGAACGCAACG‑ACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCGGTGGCAAAGCGGTAGAGGGGa        >  2:354305/1‑150 (MQ=255)                                                                                                                                           tCAATTTGGTAACCCCGCTGAACGCAACG‑ACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCGGTGGCAAAGCGGTAGAGGGGAAAgg    >  2:1035802/1‑150 (MQ=255)                                                                                                                                             aaTTTGGTAACCCCGCTGAACGCAACG‑ACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCAAAAAGCGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCGGTGGCAAAGCGGTAGAGGGGAAAGGat  >  2:892307/1‑150 (MQ=255)                                                                                                                                                  |                                                                                                                                                  TGAACTGGCAGTCAAAGCCATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACG‑ACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCGGTGGCAAAGCGGTAGAGGGGAAAGGAT  >  NC_000913/1489029‑1489318

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