Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913501,5640AC50.0% 80.1 / 20.6 58*559E (TAA→GAA) ybaLinner membrane putative NAD(P)‑binding transporter
Reads supporting (aligned to +/- strand):  ref base A (26/3);  new base C (0/29);  total (26/32)
Fisher's exact test for biased strand distribution p-value = 3.30e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATA  >  NC_000913/501421‑501706
                                                                                                                                                |                                                                                                                                              
gATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTc                                                                                                                                         >  1:966138/1‑151 (MQ=255)
       aaGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCaccacc                                                                                                                                   >  2:333860/1‑150 (MQ=255)
       aaGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCt                                                                                                                                  >  1:747652/1‑151 (MQ=255)
              cgtcgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACccgc                                                                                                                           >  2:547490/1‑151 (MQ=255)
                          ggCCAGGTAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCCCCGGTGGCGTTTc                                                                                                               >  2:963890/1‑151 (MQ=255)
                              aGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCcagc                                                                                                           >  2:351460/1‑151 (MQ=255)
                                 tAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGt                                                                                                        >  2:1002744/1‑151 (MQ=255)
                                        tccccgccgcATCTGGCTTTTGCCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATCCCTACATTTTTTTCATATTTTACACCCGGCACCCCCCGTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCa                                                                                                    <  1:182538/144‑1 (MQ=37)
                                           ccgccgcATTGGGTATTTGCCTCCATCTGCCTGATGCAGCGCGGCCGCATCTTATCATGCCACCATTTTTTTTATTTTTACCTCCGGGAACCCACGTTTTCCCCCGTCCCCCCCACACCCGCCGGTGGCGTTTCCAGCAGTTcc                                                                                                     <  2:520907/143‑1 (MQ=255)
                                           aCGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATgg                                                                                              >  2:241986/1‑151 (MQ=255)
                                           aCGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATgg                                                                                              >  1:811759/1‑151 (MQ=255)
                                              ccgttgcTGGCATTTCTCTCTCGATACCGGTTGCGACGGTGGCGGGTCTTAACAAGCCTCCATATTTTACATATTTTTCATCCGGCAACCTCCTTTTCCCCCGTCCCCCCCTCTCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGt                                                                                             <  1:1198479/143‑1 (MQ=255)
                                              ccgcATCTGGCTTTTGGCTCTCGATGCCTGATGCGACGCTGGCGCTTCTTATCATGCCTACTTTTTTTTCATTTTTTACATCCGGCAACCCCCTTTTCCCCCGTCACCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTAc                                                                                           <  1:447647/151‑1 (MQ=38)
                                                 cATCTGTCTTTTGTCCCTCGTCGCCGGAGGCGCCGCCGGCGCGTTTTCTCGTCCCTAATTTTTTTTCATTTTTTTCATCCGCCAACCACTGTTGCCCCTGCCGCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATgg                                                                                              <  2:1201909/145‑1 (MQ=255)
                                                   tCTGGCTTTTGGCTCTCGATGCCTGATGCGACGCTGGCGCTTTTTATCATCCCTACATTTTTTTCATTTTTTCCACCCGGCAACCCCCTTTTACCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACggg                                                                                        <  2:1170954/149‑1 (MQ=255)
                                                    cTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTAACCCGTCACCACCTCACCCGCCGGCGGCGTTTCCAGCAGATCCAGCATGGTACGGGCg                                                                                      >  2:792507/1‑150 (MQ=255)
                                                     tGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGAt                                                                                    >  1:362765/1‑151 (MQ=255)
                                                      ggCTTTTGGCTCTCTATCCTGGATGCGACGCTGGCCCGTTTTATCATGCCTACTTTTTTTCCATTTTTTACACCCGCCACCCCCCGTTACCCCCTCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGAt                                                                                    <  2:79697/150‑1 (MQ=255)
                                                      ggCTTTTGGCTCTCGATGCCTGATGCGACGCGGGCGCTTCTTATCATGCCTACATATTTTTCATATTTTACACCCGGCACCCACCTTTTCCCCCGTCACCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGAt                                                                                    <  1:845602/150‑1 (MQ=255)
                                                      ggCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGAtt                                                                                   >  1:1224315/1‑151 (MQ=255)
                                                          tttGTCTCCCGATGCCTTACGCGACGCGGGCGCGTCTTATCTGCCCTACATATTTTTAATATTTTACACCCGCCAACCACCGTTTCCCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGAttt                                                                                  <  1:1025726/148‑1 (MQ=255)
                                                          tttGGCTCTCGTTGCCTGATGCGCCGCGGGCGCGTCTTATCATGCCTACATTTTTTTCATTTTTTCCACCCGGCACCCCCCTTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCa                                                                                <  1:125289/150‑1 (MQ=255)
                                                          tttGGCTCTCGATGCCTGATGCGCCGCGGGCGCTTCTTTTCATCCCTACATTTTTTTCATTTTTTACACCCGGCACCCCCCTTTTCCCCCGCCCCCCCCGCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTc                                                                                 <  2:11358/149‑1 (MQ=255)
                                                                  tcGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGccc                                                                       >  2:963978/1‑151 (MQ=255)
                                                                        ccTGATGCGACGCGGGCGCGTTTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCCCCTTTTCCCCCGCCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCAc                                                                  <  1:345796/150‑1 (MQ=255)
                                                                          tGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTAc                                                               >  1:893082/1‑151 (MQ=255)
                                                                              gCGACGCGGGGGGTTTTTATCATGCCTACATATTTTTCATATTTTACACCCGGCACCCCCCGTTTCCCCCGCCACCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTg                                                            <  2:793154/150‑1 (MQ=255)
                                                                                gACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGAtt                                                         >  1:1128184/1‑151 (MQ=255)
                                                                                      ggCGCTTTTTATCATGCCTACATTTTTTTCATATTTTACACCCGGCAACCCCCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGc                                                      <  1:541485/148‑1 (MQ=255)
                                                                                      ggCGCGTTTTATCATCCCTACTTTTTTTTAATTTTTTACACCCGCCACCCCCCTTTTACCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcac                                                    <  2:730483/150‑1 (MQ=255)
                                                                                        cgctttttATCATGCCTACATATTTTTAATTTTTTCCATCCGCCAACCCCCGTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcacca                                                  <  1:582016/144‑1 (MQ=255)
                                                                                               tATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGt                                          >  1:743163/1‑151 (MQ=255)
                                                                                               tATCATGCCTACATATTTTTCATATTTTACACCCGGCACCCCCCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTc                                             <  1:862260/148‑1 (MQ=255)
                                                                                               tATCATGCCTACATATTTTTCATATTTTACACCCGGCACCCCCCGTTTCCCCCGTCCCCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTc                                             <  2:468504/148‑1 (MQ=255)
                                                                                                aTCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTg                                         >  1:39453/1‑151 (MQ=255)
                                                                                                 tCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTg                                         >  1:8469/1‑150 (MQ=255)
                                                                                                 tCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTg                                         >  1:534636/1‑150 (MQ=255)
                                                                                                  cATGCCTACATTTTTTTCATATTTTCCACCCGGCACCCCCCTTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGt                                          <  1:918004/148‑1 (MQ=255)
                                                                                                      ccTACATTTTTTTCATTTTTTACATCCGGCACCCACCTTTTACCCCGTCACCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGata                                      <  2:656535/148‑1 (MQ=255)
                                                                                                         aCATTTTTTTCATATTTTACACCCGACAACCCCCGTTTCCCCCTTCACCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATAc                                   <  1:364875/148‑1 (MQ=255)
                                                                                                          cATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGcc                                >  2:582377/1‑150 (MQ=255)
                                                                                                              tttttCATTTTTTACCTCGGGCAACCCCCGTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACg                                  <  1:911214/144‑1 (MQ=255)
                                                                                                              tttttCATATTTTCCATCCGGCAACCCCCTTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGcc                                <  1:675793/146‑1 (MQ=255)
                                                                                                                 ttctttttttccccccGGCACCCCCCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCAATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCAt                         <  2:188934/136‑1 (MQ=255)
                                                                                                                    atatTTTACATCCGCCACCCCCCTTTCCCCCCGCCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGt                      <  1:1092294/150‑1 (MQ=255)
                                                                                                                     tatTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTc                     >  1:537697/1‑150 (MQ=255)
                                                                                                                       ttttACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCAt                   >  2:881282/1‑150 (MQ=255)
                                                                                                                       tttaccacccGCCACCCACCGTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTc                     <  2:11179/140‑1 (MQ=255)
                                                                                                                        tttACCCCCGGCACCCCCCTTTTCCCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATa                  <  2:288781/150‑1 (MQ=255)
                                                                                                                          tACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAAt                >  2:174722/1‑150 (MQ=255)
                                                                                                                          tACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATg               >  1:730805/1‑151 (MQ=255)
                                                                                                                          tACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATg               >  2:278265/1‑151 (MQ=255)
                                                                                                                               ccGGCACCCCCCTTTTCCCCCGTCACCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgcgc          <  1:915026/151‑1 (MQ=255)
                                                                                                                               ccGCCACCCCCCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgcg           <  1:327015/150‑1 (MQ=255)
                                                                                                                                cGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGcgcg         >  1:539700/1‑151 (MQ=255)
                                                                                                                                    caccccccttttCCCCCGTCACCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGcgcg         <  2:741004/139‑1 (MQ=255)
                                                                                                                                     accccccttttCCCCCATCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGc      <  2:951021/142‑1 (MQ=255)
                                                                                                                                        ccccGTTTCCCCCGCCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaata  <  2:202919/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                              
GATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATA  >  NC_000913/501421‑501706

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: