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

GCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATT‑AAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAT  >  NC_000913/4604209‑4604481
                                                                                                                                       |                                                                                                                                          
gCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCGATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttttttt                                                                                                                                        >  1:39270/1‑140 (MQ=255)
 cGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttcttttt                                                                                                                                       <  1:780059/140‑1 (MQ=255)
   tCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTttttctttttct                                                                                                                                     <  1:418285/140‑1 (MQ=255)
     aTGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGc                                                                                                                                   <  2:279430/140‑1 (MQ=255)
     aTGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGc                                                                                                                                   <  2:13023/140‑1 (MQ=255)
                 ggTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCATGAGaaaggg                                                                                                                       >  1:464755/1‑136 (MQ=255)
                  gTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAATGAg                                                                                                                      >  2:511149/1‑140 (MQ=255)
                    aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAgcg                                                                                                                    <  2:893685/140‑1 (MQ=255)
                       gagaGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTTCATTAGAAGGGAGCGCag                                                                                                                 >  2:914618/1‑140 (MQ=255)
                       gagaGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGGAGTggggggag                                                                                                                 >  2:241645/1‑133 (MQ=255)
                       gagaGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGGGCag                                                                                                                 >  2:682499/1‑140 (MQ=255)
                       gagaGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCag                                                                                                                 >  2:956677/1‑140 (MQ=255)
                       gagaGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCCTTAGAAGTGAGCGCag                                                                                                                 >  2:80343/1‑140 (MQ=255)
                        agagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAATTTTTCTTTTTTTGCCATGAGAAGTGAGCGCaga                                                                                                                >  2:680352/1‑140 (MQ=255)
                         gagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagag                                                                                                               <  2:61821/140‑1 (MQ=255)
                           gTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGTGAGGGGagagag                                                                                                             >  1:264340/1‑140 (MQ=255)
                            tttATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTTATTTTTTCTTTTTTTGCCATGGGAAGGGAGCGCAGAGAgg                                                                                                            >  1:438367/1‑140 (MQ=255)
                             ttATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGTGAGGGCAGAGAggg                                                                                                           >  2:407404/1‑140 (MQ=255)
                             ttATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAggg                                                                                                           <  2:493895/140‑1 (MQ=255)
                               aTTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGGGGGCagggaggggt                                                                                                         >  2:628133/1‑138 (MQ=255)
                                      aaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCAGGGGGGGGTTATTTTg                                                                                                  >  2:874326/1‑140 (MQ=255)
                                      aaGTTCGCGGGTAACCCTTATAAGCTCCATGCCTATACAAGTCAGGAGTCATTTCAGGAAGCCATGTCGCGGATATCGTTTGCGGCGGTAATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTg                                                                                                  <  2:620877/140‑1 (MQ=255)
                                       aGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGc                                                                                                 <  2:702510/140‑1 (MQ=255)
                                                 tAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAAc                                                                                       <  1:407404/140‑1 (MQ=255)
                                                 tAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAAc                                                                                       <  2:191386/140‑1 (MQ=255)
                                                 tAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAAc                                                                                       <  2:891659/140‑1 (MQ=255)
                                                         ttaAGCTCCATTCCTATTCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGAtt                                                                               <  2:874911/139‑1 (MQ=255)
                                                            aGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTGCATGAGAAGAGGGGGCAGGGGGGGGGTATTTTTTCTTGCTGAAATGGGGATT‑TAg                                                                           >  2:926387/1‑140 (MQ=255)
                                                                ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATT‑AAGTTTc                                                                       <  1:291164/140‑1 (MQ=255)
                                                                        ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTTTTTTTTTTTTTTTTCCATGAGAAGGGAGGGCAGAGGGGGGTTTTTTTGCCTGACTGAACTGGCGGTT‑AAGTTTTCGCGTTcc                                                               >  2:341038/1‑140 (MQ=255)
                                                                        ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTCTGCCATGAGAAGTGAGGGCAGAGAGGGGTTATCTTTCCTGACTGAACTGGCGATT‑AAGTTTCCGCGTAcc                                                               >  2:589609/1‑140 (MQ=255)
                                                                         taCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGGTT‑AAGTTTTCCCGTAccc                                                              >  2:358279/1‑140 (MQ=255)
                                                                         taCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGGGAAGTGAGCGCAGAGAGGGGTTATCTTGCCTTACTGAACTGGCGATT‑AAGTTTTCGCGTAccc                                                              >  2:434163/1‑140 (MQ=255)
                                                                         taCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTGCCATTAGAAGTGAGCGCCGAGAGGGGTTATCTTTCCTGACTGAACTGGCGAGT‑AAGTTTCCGCGTAccc                                                              >  2:273423/1‑140 (MQ=255)
                                                                              gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTTTTTTTTTTTTTTTTTCATTGGAAGTGAGCGCAGAGGGGGGGTTTTTTTCCTGACTTAACTGGCGGTT‑AAGTTTTCGCGTACCCCGCGt                                                         >  2:865610/1‑140 (MQ=255)
                                                                              gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTTCCATGAGGAGTGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGATT‑AAGTTTTCGCGTACCCGGCGt                                                         >  2:150920/1‑140 (MQ=255)
                                                                              gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCACTTTTTTTTTTTCTGCCCTGCGAAGAGAACGCGGAGAGGGGGTATCTTTTCTGCCTGTAATGGCGCGT‑AAGTTTCCGCGTACCCcgcgg                                                         >  2:228220/1‑139 (MQ=255)
                                                                              gTCAGGAGTCATTTAAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCAGAGAGGGGTTATCTTTCCTGACTGAACTGGGGATT‑AAGTTTCCCCGTACCCGGCgg                                                         >  2:30448/1‑139 (MQ=255)
                                                                               tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGATT‑AAGTTTCCGCGTACCCcgcggt                                                        >  2:189966/1‑138 (MQ=255)
                                                                                 aGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTTATTTTTTTTTTTTTTTCCATTAGAAGGGAGCGCCGAGAGGGGTTATCTTGCCTGACTGAACTGGCGATT‑AAGTTTTCGCCTACCCCGCGtttt                                                      >  2:525995/1‑139 (MQ=255)
                                                                                   gAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCAGAGGGGGGGTATCTTGCCTGGCTGAAATGGCGATT‑AAAATTTCTCGTACCCCGCGgttttt                                                    >  2:513891/1‑137 (MQ=255)
                                                                                       cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATT‑AAGTTTCCGCGTACCCGGCGTTTAGTTAtt                                                <  1:196286/140‑1 (MQ=255)
                                                                                        aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATT‑AAGTTTCCGCGTACCCGGCGTTTAGTTATTg                                               <  1:511149/140‑1 (MQ=255)
                                                                                          ttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATT‑AAGTTTCCGCGTACCCGGCGTTTAGTTATTGCg                                             <  1:874326/140‑1 (MQ=255)
                                                                                            cAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATT‑AAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGga                                           <  2:199277/140‑1 (MQ=255)
                                                                                            cAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATT‑AAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGga                                           <  2:164627/140‑1 (MQ=255)
                                                                                             aGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCAGAGAGGGGTTATCTTTCCTGACTTAACTGGCGATTAAAGTTTCC‑CGGACCCGGCGGTTTGTTATTTCgggt                                          >  2:383736/1‑138 (MQ=255)
                                                                                             aGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGGGGGGGGTTATCTTGCCTGACTGAAATGGCGGTT‑AAGTTTCCGCGTACCCGGCGTTTTGTTATTTCgggt                                          >  2:675947/1‑138 (MQ=255)
                                                                                                aTGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATT‑AAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatgat                                       <  1:732759/140‑1 (MQ=255)
                                                                                                     aTGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATT‑AAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGatat                                  <  1:189966/140‑1 (MQ=255)
                                                                                                          gcgGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATT‑AAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAg                             <  2:730142/140‑1 (MQ=255)
                                                                                                                                     ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATT‑AAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAt  <  2:438367/140‑1 (MQ=255)
                                                                                                                                     ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATT‑AAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAt  <  2:841120/140‑1 (MQ=255)
                                                                                                                                     ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATT‑AAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAt  <  2:39270/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                          
GCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATT‑AAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAT  >  NC_000913/4604209‑4604481

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