Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,604,3440CT34.2% 53.1 / 16.5 41S62F (TCT→TTT) bglJDNA‑binding transcriptional regulator BglJ
Reads supporting (aligned to +/- strand):  ref base C (5/22);  new base T (14/0);  total (19/22)
Fisher's exact test for biased strand distribution p-value = 3.30e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCA  >  NC_000913/4604209‑4604471
                                                                                                                                       |                                                                                                                                
gCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttttttt                                                                                                                              >  2:573486/1‑140 (MQ=255)
gCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttttttt                                                                                                                              >  1:440226/1‑140 (MQ=255)
gCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctttt                                                                                                                              >  1:625439/1‑140 (MQ=255)
 cGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttcttttt                                                                                                                             <  1:573486/140‑1 (MQ=255)
 cGTCATGAGTAGTATCGGTATAGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttcttttt                                                                                                                             <  2:476323/140‑1 (MQ=255)
  gTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAttttttttttttt                                                                                                                            >  1:949177/1‑139 (MQ=255)
   tCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTttttctttttct                                                                                                                           <  1:874831/140‑1 (MQ=255)
    cATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTTAttttttttttttttt                                                                                                                          >  2:874831/1‑137 (MQ=255)
     aTGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGc                                                                                                                         <  2:787375/140‑1 (MQ=255)
     aTGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGc                                                                                                                         <  2:200678/140‑1 (MQ=255)
          tagtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATga                                                                                                                    <  1:277690/140‑1 (MQ=255)
           agtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGc                                                                                                                         <  2:836785/134‑1 (MQ=255)
           agtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGc                                                                                                                         >  1:836785/1‑134 (MQ=255)
               tCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGt                                                                                                               <  2:949177/140‑1 (MQ=255)
                 ggTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTttttctttttct                                                                                                                           >  1:776618/1‑126 (MQ=255)
                 ggTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTttttctttttct                                                                                                                           <  2:776618/126‑1 (MQ=255)
                  gTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGTCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGTGAg                                                                                                            >  1:270701/1‑140 (MQ=255)
                    aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAgcg                                                                                                          >  2:819127/1‑140 (MQ=255)
                         gagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagag                                                                                                     <  2:625439/140‑1 (MQ=255)
                          agTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTTCCCTGAGAAGTGAGCGgagaga                                                                                                    >  2:54848/1‑140 (MQ=255)
                            tttATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGTGAGGGCAGAGAgg                                                                                                  >  2:539457/1‑140 (MQ=255)
                            tttATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGTGAGCGCAGAGAgg                                                                                                  >  1:6081/1‑140 (MQ=255)
                             ttATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAggg                                                                                                 <  2:190954/140‑1 (MQ=255)
                             ttATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAggg                                                                                                 <  1:290875/140‑1 (MQ=255)
                             ttATTCAGAAAGTTTGCGGGGAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCACGTTTGCGGCGGTGCTGTTTTGTTTTTCTGCGATGAGAAGCGGGCGCGGAGCggg                                                                                                 >  2:650196/1‑140 (MQ=255)
                                  cAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCACGAGAAGTGAGGGCAGAGAGGGATTAt                                                                                            >  1:532340/1‑140 (MQ=255)
                                      aaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATCTTTTCTTTTTTTCCCATGCAAAATGAGCGCAAAGAGGGGTTAttttt                                                                                        >  1:422087/1‑139 (MQ=255)
                                                           aaGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCGGAGAGGGATTATTTTTCCTGACTGGACT‑GGCGATTaa                                                                  >  2:269481/1‑140 (MQ=255)
                                                               tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGttt                                                              >  2:68911/1‑140 (MQ=255)
                                                                 cATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCTATGGTAAGGGATCGCAGAGCTGGATTTTTTTGCCTTGCTGCACT‑TGCGATTAAGGTTcc                                                            >  2:578669/1‑140 (MQ=255)
                                                                   tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAATGAGCGCAGAGAGGGGTTATTTTTCCTGTCTGAACT‑GGCGATTAAATTTTcgc                                                          >  2:493693/1‑140 (MQ=255)
                                                                         taCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGAGGGCGCAGAGAGGGGTTATCTTTCCCGACTGAAATGGGGGTTTAAGTTTcccggaacc                                                     >  2:698891/1‑133 (MQ=255)
                                                                         taCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTAccc                                                    <  2:680499/140‑1 (MQ=255)
                                                                          aCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTACCCg                                                   <  1:536018/140‑1 (MQ=255)
                                                                              gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTTCCCCGAGAAGTGAGCGCAGAGAGGGATTTTCTTTCCTGACTGAACT‑GGCGATTAAGTTTCCCCGTTCCCcgcgg                                               >  2:102184/1‑139 (MQ=255)
                                                                              gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTACCCGGCGt                                               <  1:426479/140‑1 (MQ=255)
                                                                               tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTCTTTTTCTTTCATGAGAAGTGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACT‑GGCGGTTAAGTTTTCGCGTACCCcgcggt                                              >  2:18486/1‑138 (MQ=255)
                                                                               tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGGGAAAGGAGCGCAGAGAGGGGTTTTCTTGCCTGACTGAACT‑GGCTATTAAGTTTCCGCGGAccccggggt                                              >  2:592363/1‑134 (MQ=255)
                                                                                   gAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCAGGGAGGGGTTATCTTGCCTGACTGAAAT‑GGCGATTAAGTTTTCCCGTACCCCGCGGTTAGt                                          >  2:806246/1‑140 (MQ=255)
                                                                                   gAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGCCATTTTTTCTTTTTCTTCCCCGAGAAGTGAGCGCCGAGAGGGGTTATCTTTCCCGACCGAACC‑GGCGCTTAAGTTTCCCCCTACCCCGCGCTTTGt                                          >  2:426479/1‑140 (MQ=255)
                                                                                     gTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTa                                        <  1:698891/140‑1 (MQ=255)
                                                                                     gTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTa                                        <  1:650196/140‑1 (MQ=255)
                                                                                       cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTAtt                                      <  1:493693/140‑1 (MQ=255)
                                                                                        aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCCTGAGAAGTGAGCGCAGAGAGGGATTTTCTTGCCTGACTGAACT‑GGCGGTTAAATTTCCGCGGACCCGGCGTTTAGTttttt                                     >  2:246586/1‑139 (MQ=255)
                                                                                        aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTg                                     <  1:102184/140‑1 (MQ=255)
                                                                                          ttCAGGATGTCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCg                                   <  2:270701/140‑1 (MQ=255)
                                                                                                aTGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCTATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatgat                             <  1:592363/140‑1 (MQ=255)
                                                                                                 tGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatgatg                            <  1:959643/140‑1 (MQ=255)
                                                                                                  gCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGTGggtgatga                           >  2:888916/1‑140 (MQ=255)
                                                                                                   ccATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGGGAGCGCAGAGGGGGGTTATCTTTCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTACCCGGCGGTTAGTTATTTCGgatggtgat                          >  2:59060/1‑140 (MQ=255)
                                                                                                    cATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTTCCATGAGAAGTGAGCGCATAGAGGGGTTATCTTGTCTGGCTTAACT‑GGCGATTAAGTTTTCCCGTACCCGGCGGTTTGTTTTTTTGGATGGTGgta                         >  2:790147/1‑140 (MQ=255)
                                                                                                            ggATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCt                 <  1:59060/140‑1 (MQ=255)
                                                                                                                    tttGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGAt         <  1:197785/140‑1 (MQ=255)
                                                                                                                        cggcggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGt     <  2:620968/140‑1 (MQ=255)
                                                                                                                           cggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCa  <  1:819127/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                
GCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCA  >  NC_000913/4604209‑4604471

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