Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,604,3440CT47.6% 33.1 / 19.6 42S62F (TCT→TTT) bglJbgl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (5/17);  new base T (20/0);  total (25/17)
Fisher's exact test for biased strand distribution p-value = 1.08e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGG  >  NC_000913/4604197‑4604494
                                                                                                                                                   |                                                                                                                                                      
tCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCtttttttttt                                                                                                                                                     >  2:657008/1‑151 (MQ=255)
 cATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCtttttttttttt                                                                                                                                                    >  2:727856/1‑151 (MQ=255)
 cATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttcttt                                                                                                                                                     <  1:228685/150‑1 (MQ=255)
   tAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTtttttctttttt                                                                                                                                                  >  1:514122/1‑150 (MQ=255)
       aaaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGcc                                                                                                                                              >  2:704166/1‑151 (MQ=255)
         aaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGc                                                                                                                                               >  1:726971/1‑148 (MQ=255)
           tGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAtttttttttttttttccctg                                                                                                                                           >  1:280590/1‑142 (MQ=255)
            gCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATga                                                                                                                                          <  2:421215/150‑1 (MQ=255)
             cGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTTCCATg                                                                                                                                           >  1:577917/1‑148 (MQ=255)
                cATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGt                                                                                                                                     <  2:490557/151‑1 (MQ=255)
                         taTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGGGAgcg                                                                                                                                >  2:143084/1‑147 (MQ=255)
                          aTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagag                                                                                                                           <  1:488118/151‑1 (MQ=255)
                                 ttGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGGCATGTGGAGGGAGAGCAGAGGGGGGtt                                                                                                                    >  1:495232/1‑151 (MQ=255)
                                         ttATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTTAGGAGGGAGCGCAGCGGGGGGttttttttg                                                                                                             >  1:114887/1‑143 (MQ=255)
                                          tATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGcc                                                                                                            <  2:582598/150‑1 (MQ=255)
                                           aTTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctg                                                                                                          <  2:463475/151‑1 (MQ=255)
                                              cAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctgact                                                                                                       <  2:637494/151‑1 (MQ=255)
                                               aGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGACCATTTTTCTTTTTTTCCCCTGGGCAGGGCGCGCAGAGAGGGATTATCTTGCctgaa                                                                                                        >  1:547405/1‑148 (MQ=255)
                                                gAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctgactga                                                                                                     <  2:679582/151‑1 (MQ=255)
                                                 aaaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGTGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGaa                                                                                                    <  2:492042/151‑1 (MQ=255)
                                                  aaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGGTTATCTTGCCTGACTGaa                                                                                                    >  1:207419/1‑150 (MQ=255)
                                                    gTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTg                                                                                                 <  1:362405/151‑1 (MQ=255)
                                                       tGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGc                                                                                               <  2:473150/150‑1 (MQ=255)
                                                                     ataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGATTATTTTGCCTGACTGAACTGGCGGTTAAGTTTTcgcg                                                                                >  2:556380/1‑151 (MQ=255)
                                                                            ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTTCGCGGACCCgg                                                                         >  1:345231/1‑151 (MQ=255)
                                                                            ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCgg                                                                         <  2:510210/151‑1 (MQ=255)
                                                                            ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGGGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGATTAAGTTTTCGCGTACCCgg                                                                         >  2:162622/1‑151 (MQ=255)
                                                                             cATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGATTAAGGTTTCGCGTAccc                                                                           >  2:122000/1‑148 (MQ=255)
                                                                               tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGGGAAATGGGCGCCGGGAGGGGTTTTCTTGGCTGACTTAAACGGCGATTAAGGTTCCGCGGACCCcgcg                                                                       >  1:244671/1‑150 (MQ=255)
                                                                               tACATATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCCTGAGAAGTGAGCGCAGAGAGGGGTTTTCTTTGCTTACTGAACTGGCGGTTAAGATTTCGCGGACCCgg                                                                         >  1:666173/1‑148 (MQ=255)
                                                                                 ccTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGGTTAAGTTTTCGCGTACCCGGCGgtt                                                                    >  1:335809/1‑151 (MQ=255)
                                                                                    ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAg                                                                  <  2:60350/150‑1 (MQ=255)
                                                                                      aCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGtt                                                                <  2:282892/150‑1 (MQ=255)
                                                                                      aCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGtt                                                                <  2:144828/150‑1 (MQ=255)
                                                                                      aCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGtt                                                                <  1:417772/150‑1 (MQ=255)
                                                                                       ccAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGGAGTGAGCGGCGAGGGGGGTTATCTTGTCTGTCTGGAATGGGGGTTTAGTTTTCGCGGACCCCGCGGTTTGtt                                                                >  1:221985/1‑149 (MQ=255)
                                                                                            taGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTTATTTTTTTTTTTTCTGCCATGAGGAGGGAGCGGAGGGGGGGGTTATTTTTTCCTACTGAAACGGCGGTTAAGGTTTCGCGTACCCCGCGTTTTTTTATTGCgg                                                         >  1:680202/2‑151 (MQ=255)
                                                                                              ggAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgat                                                       <  2:630519/151‑1 (MQ=255)
                                                                                                             tGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGGAGGGAGCGCAGAGAGGGGTTTTCTTGCCTTACTGAACTGGCGGTTTAGTTTTCGCGGAACCGGGGGTTTGTTATTGCGGATGATGATATTTaaac                                         >  1:663513/1‑148 (MQ=255)
                                                                                                                  tGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGTGAGGGCAGAGAGGGGTTATTTTTCCTGACTGAACTGGCGATTAAGTTTTCGCGTACCCGGCGGTTAGTTTTTTCGGGTGATGATATTGAAGCTCGGc                                    >  2:78328/1‑150 (MQ=255)
                                                                                                                            ctcGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGGGAGCGCAGAGAGGGATTATCTTGCCTGACTGAAATGGCGATTAAGTTTTCGCGTACCCGGCGTTTTGTTATTGCGGGTGATGATATCGAAGCTCGGCTGATTGGCTc                          >  2:656524/1‑150 (MQ=255)
                                                                                                                             tcGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGGTTATCTTTCCTTTCTGAACTGGGGGTTAAGTTTTCGCGTACCCGGCGTTTTGTTTTTGGGGGTGATGATATTGAAGCTCGGCTGGTTGGTTCAt                        >  1:17107/1‑151 (MQ=255)
                                                                                                                              cGTTTGCGGCGGTCATTTTTTCTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTAGCGATTAAGTTTCCGCGTACCCGGCGTTTTGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTTAtt                       >  2:163544/1‑151 (MQ=255)
                                                                                                                              cGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCAt                        <  2:588044/150‑1 (MQ=255)
                                                                                                                                tttGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGTGAGCGCAGAGAGGGATTTTTTTTCCTGACTGAACTGGCGGTTTAGTTTTCGCGTACCCGGCGTTTTGTTTTTGTGGGTGATGATATTTAAGCTCCGGTGCTTGGTTCAtttt                     >  1:397554/1‑149 (MQ=255)
                                                                                                                                tttGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGGAGGGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGGGATTAAGTTTTCCCGTACCCGGCGTTTTGTTATTGCGGGTGATGGTATTGAAGCTCGGGTGATTGGTTCATTGt                     >  2:459817/1‑151 (MQ=255)
                                                                                                                                             ttttttCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCt        >  2:164655/1‑151 (MQ=255)
                                                                                                                                                   ctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCAGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACgg  >  1:167481/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
TCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGG  >  NC_000913/4604197‑4604494

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