Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,604,3500CT46.0% 29.4 / 16.3 37S64F (TCT→TTT) bglJbgl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (5/15);  new base T (17/0);  total (22/15)
Fisher's exact test for biased strand distribution p-value = 3.12e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTA  >  NC_000913/4604201‑4604498
                                                                                                                                                     |                                                                                                                                                    
aGAAAAATGCGTCATGAGTAGTATCGGTAGTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctttttc                                                                                                                                                      <  2:80579/150‑1 (MQ=255)
                agtagtATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAgc                                                                                                                                     <  2:353066/151‑1 (MQ=255)
                     taTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGGgcg                                                                                                                                    >  1:690972/1‑147 (MQ=255)
                        cGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATTAGAAGTGAGCGCagag                                                                                                                               >  1:159236/1‑149 (MQ=255)
                             ttGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTTCCCTGAGAAGGGAGCGCAGGGAGGGGtt                                                                                                                        >  2:576408/1‑151 (MQ=255)
                               gagaGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGGCAGGAGTCATTTCAGGATGCCATGGCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTAt                                                                                                                      <  2:174830/151‑1 (MQ=255)
                                agagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTAt                                                                                                                      <  2:460380/150‑1 (MQ=255)
                                 gagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGGGAAGGGAggggaggggggggtt                                                                                                                        >  1:296709/1‑133 (MQ=255)
                                        ttCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCAGAGAGGGATTATCTTGCct                                                                                                               >  1:190293/1‑149 (MQ=255)
                                          cAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctgac                                                                                                            <  2:436111/150‑1 (MQ=255)
                                           aGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGCCATGAGAAAGGAGGGCCGAGAGGGGGTATCTTGCctgactg                                                                                                          >  1:370020/1‑151 (MQ=255)
                                           aGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctgact                                                                                                           >  1:712982/1‑150 (MQ=255)
                                             aaaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTTCctgactga                                                                                                         >  1:706358/1‑150 (MQ=255)
                                                gTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTg                                                                                                     <  1:178794/151‑1 (MQ=255)
                                                            cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGCCATGAGGAGAGGGGGGGGAGAGGGTATATTTTGCCTGACTGGACTGGCGGTTAAGttt                                                                                         >  1:155058/1‑151 (MQ=255)
                                                             ccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGATTAAggt                                                                                          >  1:238005/1‑147 (MQ=255)
                                                                         cATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGc                                                                            <  1:182451/151‑1 (MQ=255)
                                                                            aCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGt                                                                          <  1:139928/150‑1 (MQ=255)
                                                                               tataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCTCTTTTTTTTTTTCTTCCCTGTGGAGGGGGGGCCGAGGGGGGTTATTTTGCCCGACTCAAATGGGGTTTTAGGTTTcgcggtaccggtggtt                                                                        >  1:569448/1‑137 (MQ=255)
                                                                                 taCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTGGCATGAGAAGGGAGCGCAGCGGGGGATTTTTTTTCCTGACTGAACTGGCGGTTAAGTTTTCGCGTACCCGGCGGTTAGtt                                                                    >  2:65741/1‑151 (MQ=255)
                                                                                   ccAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTAt                                                                  <  1:652402/151‑1 (MQ=255)
                                                                                     aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTCTTCCATGAGAAGTGAGCGCAGAGGGGGGTTTTCTTGCCTGACTGAAATGGGGATTAAGTTTCCGCGTACCCCGCGGTTTGTtttt                                                                 >  1:507878/1‑150 (MQ=255)
                                                                                     aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGACATTTTTTTTTTTTTTGCGCTGGGAAAGGAGCGCAGGGGGGGATTTTTTTGGCTGACTTAACTGGGGATTAAGGTTGCGGGGACCCgggggtttgt                                                                     >  2:706881/1‑138 (MQ=255)
                                                                                       tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGTCCTGTGGAGGGGGGGCTGAGGGGGGTTATTTTTGCTTACTGAACTGGCGGTTTAGTTTTCGCGGACCCGGCGGTTTGTTATTGCg                                                              >  2:348965/1‑151 (MQ=255)
                                                                                          ggAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGTGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGGTTTAGGTTTCGCGCACCCGGCGTTTTGTTAtttggg                                                             >  1:316106/1‑145 (MQ=255)
                                                                                          ggAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgat                                                           <  1:204391/151‑1 (MQ=255)
                                                                                               cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTCTTTTTTTGCCCTGAGAAGTGAGCGCCGAGGGGGGGTTTTTTGCCTGTCTTAACTGGCGAATAAGATTTCGCGGACCCCGCGGTTTGTTTTTGCGggtgat                                                        >  2:267855/1‑149 (MQ=255)
                                                                                                    cAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATAtt                                                  <  1:53192/150‑1 (MQ=255)
                                                                                                      ggATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGaa                                               <  1:593344/151‑1 (MQ=255)
                                                                                                        aTGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTCTTTTTTTGCCATTAGAAGGGAGCGCAGAGAGGGATTTTTTTGCCTGACTGAAATGGCGATTAAGGTTCCGCGTACCCGGCGTTTAGTTTTTGCGGGTGATGATATTTAAGc                                             >  2:483344/1‑151 (MQ=255)
                                                                                                        aTGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCCTGAGAAGGGAGCGCAGAGAGGGGTTTTCTTTCCCGACTGAACTGGCGATTAAAGTTCCGCGTACCCGGGGTTTTGTTATTTGGGGTGGTGATATTGaa                                               >  1:19092/1‑149 (MQ=255)
                                                                                                          gCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTc                                           <  1:74885/151‑1 (MQ=255)
                                                                                                            cATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGGGAGCGCAGAGAGGGGTTATTTTTCCTGACTTAACTGGCGATTAAGTTTCCGCGTACCCCGCGTTTTGTTTTTTCGGGTGGTGATATTTAAGCTCgg                                         >  1:199460/1‑151 (MQ=255)
                                                                                                                 cgcgGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGTTGATATTGAAGCTCGGCTGa                                     <  1:114024/150‑1 (MQ=255)
                                                                                                                  gcgGATCTCGTTTGCTGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAGCTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAATCTCGGCTGa                                     <  2:264258/149‑1 (MQ=255)
                                                                                                                   cgGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGGAGTGAGGGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGATTAAGGTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGAtt                                   >  1:634634/1‑150 (MQ=255)
                                                                                                                    ggATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTTCCATGAGAAGGGAGCGCAGAGAGGGGTTATTTTGTCTGACTGAACTGGCGAGTAAGTTTTCGCGTACCCGGCGTTTTGTTATTGCGGGTGATGATATTTAAGCTCGGCTGGtttg                                 >  2:294428/1‑149 (MQ=255)
                                                                                                                                    ggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCa                 <  1:380073/151‑1 (MQ=255)
                                                                                                                                           ttttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTgg          >  2:584924/1‑151 (MQ=255)
                                                                                                                                            tttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTgg          <  1:23659/150‑1 (MQ=255)
                                                                                                                                                    tctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGAGTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCCCCGCTGGACGGTGTa  <  2:163359/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTA  >  NC_000913/4604201‑4604498

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