Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913501,5720AC50.0% 42.6 / 12.0 40V556G (GTG→GGG) ybaLinner membrane putative NAD(P)‑binding transporter
Reads supporting (aligned to +/- strand):  ref base A (17/3);  new base C (0/20);  total (17/23)
Fisher's exact test for biased strand distribution p-value = 2.57e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCG  >  NC_000913/501422‑501719
                                                                                                                                                        |                                                                                                                                                   
aTCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑‑ACCCCGTca                                                                                                                                                     >  2:684532/1‑151 (MQ=255)
 tCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑‑ACCCCGTca                                                                                                                                                     >  1:742485/1‑150 (MQ=255)
  cGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑‑ACCCCGTcac                                                                                                                                                    >  2:329248/1‑150 (MQ=255)
   gTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑‑ACCCCGTcacc                                                                                                                                                   >  1:502540/1‑150 (MQ=255)
         ttATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑‑ACCCCGTCACCACCTCAc                                                                                                                                            <  2:546860/151‑1 (MQ=255)
                    tCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑‑ACCCCGTCACCACCTCACCCGCCGGTgg                                                                                                                                  >  2:497549/1‑150 (MQ=255)
                    tCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑‑ACCCCGTCACCACCTCACCCGCCGGTgg                                                                                                                                  >  1:44904/1‑150 (MQ=255)
                         ggCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑‑ACCCCGTCACCACCTCACCCGCCGGTGGCGttt                                                                                                                             >  1:606334/1‑150 (MQ=255)
                         ggCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑‑ACCCCGTCACCACCTCACCCGCCGGTGGCGttt                                                                                                                             >  1:289126/1‑150 (MQ=255)
                              gATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCcagca                                                                                                                       >  2:241770/1‑151 (MQ=255)
                              gATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCcagca                                                                                                                       >  2:575475/1‑151 (MQ=255)
                              gATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCcagca                                                                                                                       >  1:391959/1‑151 (MQ=255)
                                                 ttttggtttttGGATCTCGCTGCTGAATGGGCCGTGGGCGCGTTTTTTCAGGCCTAAATTTTTTCCATTTTTTTCATCCGGCAACCCCCTTTT‑‑CCCCTGCACCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGc                                                                                                    <  2:169221/143‑1 (MQ=255)
                                                     ggCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGAtt                                                                                                >  1:713094/1‑151 (MQ=255)
                                                                      gCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCa                                                                                >  1:73642/1‑150 (MQ=255)
                                                                       cttGAGGCGACGCGGGCGCGTCTTTCCATGCTTCAATTTTTTTCATTTTTTCCA‑CCGGCCACCCCCGTTTTCCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCAc                                                                               <  1:151535/149‑1 (MQ=255)
                                                                           ttGCGCCGCTGGCGCGTCTTACCAGGCCTACATTTTTTTCATATTTTCCACCCGGCACCCCCCTTTT‑‑CCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTAcc                                                                           <  2:91441/149‑1 (MQ=255)
                                                                           aTGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTAc                                                                            >  2:293778/1‑149 (MQ=255)
                                                                              cgccgttGGCGCGTCTTACCATCCCTACATATTTTTCATATTTTACACCCGCCACCCCCCGTTT‑‑CCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGAt                                                                       <  1:853748/145‑1 (MQ=255)
                                                                               caCGCGGGCGCGTTTTATGACACCTACATTTTTTTCATATTTTACACCCGGAAACCCCCTTTA‑‑CCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGAtt                                                                      <  1:60681/150‑1 (MQ=255)
                                                                                    cggCGCTTTTTATCATGCCTACTTATTTTTCATATTTTACACCCGCCAACCACCTTTT‑‑CCCCCGTCACCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcac                                                                 <  1:264731/150‑1 (MQ=255)
                                                                                          tttttATCAGCCCTACATTTTTTTCATTTTTTCCACCCGCCACCCCCCTTTC‑‑CCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGtt                                                           <  1:497338/148‑1 (MQ=255)
                                                                                             ttATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGCTTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCgg                                                        >  2:954545/1‑151 (MQ=255)
                                                                                              tcccttcccTACATTTTTTTCATTTTTCCCACCCGCCACCCCCCGTTT‑‑CCCCCCCCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGt                                                       <  2:998140/144‑1 (MQ=255)
                                                                                              taccATGCCTACTTATTTTTCATTTTTTACTGCCGGCAACCGCCGTTT‑‑CCCCGGTCCCCACCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGt                                                       <  1:347247/148‑1 (MQ=255)
                                                                                              tATCTTCCCTACATTTTTTTCATTTTTTACCCCCGCCACCCACCGTTT‑‑CCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGt                                                       <  1:1049286/151‑1 (MQ=255)
                                                                                              tATCATGCCTACATATTTTTCATATTTTACATCCGGTAACCACCGTTT‑‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATGTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGt                                                       >  1:1035787/1‑151 (MQ=255)
                                                                                                 cATGCCTACTTTTTTTCCATTTTTTCCCCCCGCCACCCCCCTTTT‑‑CCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGat                                                    <  1:116432/151‑1 (MQ=255)
                                                                                                 cATGCCTACATTTTTTTAATTTTTTACACCCGGCACCCACCGTTT‑‑CCCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGat                                                    <  1:879498/151‑1 (MQ=255)
                                                                                                   tGCCTCAATTTTTTTCATTTTTTCCACCCGGCACCCCCCGTTT‑‑CCCCCCCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGatat                                                  <  2:567499/151‑1 (MQ=255)
                                                                                                   tGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGatat                                                  >  2:464847/1‑151 (MQ=255)
                                                                                                     ccTACTTATTTTTCATTTTTTACACCCGGCACCCCCCTTTT‑‑CCCCCGCCCCCACCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATAc                                                <  1:584491/151‑1 (MQ=255)
                                                                                                     ccTACTTATTTTACTTATTTTACATCCGGCACCCCCCGTTT‑‑CCCCCTCCCCCACCTCGCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATAc                                                <  1:716383/151‑1 (MQ=255)
                                                                                                      cTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑‑ACCCCGTCACCACCTCACCCGCCGGGGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACg                                               >  2:548946/1‑151 (MQ=255)
                                                                                                       tACATATTTTTATTATTTTACACCGGGCACCCCCCGTTT‑‑CCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGc                                              <  2:459803/151‑1 (MQ=255)
                                                                                                                 tCATTTTTTCCACCCGCCACCCACCGTTT‑‑CCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCg                                    <  1:306579/151‑1 (MQ=255)
                                                                                                                                    cccccctttt‑‑CCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaa                 <  2:906464/144‑1 (MQ=255)
                                                                                                                                        cctttt‑‑CCCCCGCCCCCTCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaataat             <  2:502756/148‑1 (MQ=255)
                                                                                                                                        ccGTTT‑‑CCCCCGTCCCCACCACACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaataat             <  1:386283/151‑1 (MQ=255)
                                                                                                                                        ccGTTC‑‑CCCCCGTCACCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaataat             <  1:399207/151‑1 (MQ=255)
                                                                                                                                         cGTTT‑‑CCCCCTCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAAtc            <  1:426277/151‑1 (MQ=255)
                                                                                                                                         cGTTT‑‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCCCGCCCATCACTACCTGATTCGCACCACGTTCGGCGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaataat             >  2:559085/1‑150 (MQ=255)
                                                                                                                                                     gTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCg  >  2:662869/1‑151 (MQ=255)
                                                                                                                                                     gTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCg  >  1:86410/1‑151 (MQ=255)
                                                                                                                                                        |                                                                                                                                                   
ATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCG  >  NC_000913/501422‑501719

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