Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,604,754 T→A 6.1% S199T (TCT→ACT)  bglJ → bgl operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,604,7540TA 6.1% 464.0 / 12.7 164S199T (TCT→ACT) bglJbgl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base T (76/77);  new base A (6/4);  total (82/82)
Fisher's exact test for biased strand distribution p-value = 7.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.37e-01

Only 100 of 165 total aligned reads displayed.

AGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTCCCCGGCAGTGCTTTCATTTCAGCGTACAATCGCCACATTGCT  >  NC_000913/4604607‑4604903
                                                                                                                                                   |                                                                                                                                                     
  ccGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGa                                                                                                                                                   >  1:1697008/1‑150 (MQ=255)
  ccGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGa                                                                                                                                                   >  1:748493/1‑150 (MQ=255)
         ttAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAggg                                                                                                                                            >  2:1170406/1‑150 (MQ=255)
         ttAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTACTGACGCAggg                                                                                                                                            >  1:1598413/1‑150 (MQ=255)
            aGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGgttg                                                                                                                                         >  2:1423766/1‑150 (MQ=255)
             gCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGgttgt                                                                                                                                        <  1:565781/150‑1 (MQ=255)
                cGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTgg                                                                                                                                     >  2:231786/1‑150 (MQ=255)
                  aCGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAg                                                                                                                                   <  2:892388/150‑1 (MQ=255)
                   cGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTCGTTGGAgg                                                                                                                                  <  2:1075105/150‑1 (MQ=255)
                    ggAGAGAGAAATATTGCGCTTTTTGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGc                                                                                                                                 >  2:1477787/1‑150 (MQ=255)
                          agaAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGa                                                                                                                           <  1:1030846/150‑1 (MQ=255)
                           gaAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGat                                                                                                                          >  2:113918/1‑150 (MQ=255)
                           gaAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGat                                                                                                                          <  1:791975/150‑1 (MQ=255)
                               tatTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTc                                                                                                                      >  2:12692/1‑150 (MQ=255)
                                  tGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAACCTGGGCGTCAGTTCTGACGCAGGGTAGTTGGCGGCCGCAGATAATCTgg                                                                                                                   >  2:1531712/1‑149 (MQ=255)
                                  tGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGt                                                                                                                   <  1:43870/150‑1 (MQ=255)
                                  tGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATAGACGAGCAGCTTAAACGCAATATCAAACCGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGt                                                                                                                   <  1:1156461/150‑1 (MQ=255)
                                    cgcTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTta                                                                                                                 <  2:428483/150‑1 (MQ=255)
                                    cgcTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTta                                                                                                                 >  1:668033/1‑150 (MQ=255)
                                     gcTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTtat                                                                                                                >  2:948931/1‑150 (MQ=255)
                                     gcTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATAATCTGTtat                                                                                                                >  2:662467/1‑150 (MQ=255)
                                     gcTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGCGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTtat                                                                                                                <  2:684656/150‑1 (MQ=255)
                                      cTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACTCAATATCAAAACGATACGTGCACATAAATTTAATGTTATGTCGAAACTGGGCGTCATTTCAGACGCAGGGTTGTTGGACGCCGCAGATATTCTGTtatt                                                                                                               >  2:613567/1‑149 (MQ=255)
                                      cTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTtatg                                                                                                               >  2:399179/1‑150 (MQ=255)
                                         tATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTtatgtat                                                                                                            <  2:1395578/150‑1 (MQ=255)
                                          aTGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTtatgtatg                                                                                                           <  2:443736/150‑1 (MQ=255)
                                           tGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGc                                                                                                          <  1:73511/150‑1 (MQ=255)
                                            gTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCg                                                                                                         <  2:179153/150‑1 (MQ=255)
                                            gTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCg                                                                                                         <  2:148566/150‑1 (MQ=255)
                                              cgcgTGTCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACACAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGc                                                                                                       <  1:836442/150‑1 (MQ=255)
                                              cgcgTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGc                                                                                                       <  2:1374679/150‑1 (MQ=255)
                                              cgcgTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGc                                                                                                       <  1:384841/150‑1 (MQ=255)
                                              cgcgTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGAACTGACGccca                                                                                                                                             <  2:1178392/112‑4 (MQ=255)
                                              cgcgTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGAACTGACGccca                                                                                                                                             >  1:1178392/1‑109 (MQ=255)
                                                  tGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTg                                                                                                   >  2:290294/1‑150 (MQ=255)
                                                  tGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTg                                                                                                   <  2:665226/150‑1 (MQ=255)
                                                   ggCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGc                                                                                                  <  1:1007248/150‑1 (MQ=255)
                                                    gCTACTCAATGACACAAATTGCCGAGCAGCTTACACGCAATATCCAAACGATCCGGGCACATAACTTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCg                                                                                                 <  1:1531712/150‑1 (MQ=255)
                                                    gCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCg                                                                                                 <  1:450815/150‑1 (MQ=255)
                                                        cTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAAc                                                                                             >  1:473935/1‑150 (MQ=255)
                                                          cAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACaa                                                                                           >  1:973972/1‑150 (MQ=255)
                                                           aaTGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAg                                                                                          >  1:582637/1‑150 (MQ=255)
                                                                cacaAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAAtg                                                                                     >  2:543815/1‑150 (MQ=255)
                                                                cacaAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAAtg                                                                                     >  2:810485/1‑150 (MQ=255)
                                                                cacaAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAAtg                                                                                     >  2:526538/1‑150 (MQ=255)
                                                                  caAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCg                                                                                                 <  1:908988/136‑1 (MQ=255)
                                                                  caAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCg                                                                                                 >  2:908988/1‑136 (MQ=255)
                                                                  caAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAAtgtg                                                                                   >  2:1115304/1‑150 (MQ=255)
                                                                   aaaTTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATgtgt                                                                                  >  2:1001614/1‑150 (MQ=255)
                                                                    aaTTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGtt                                                                                 >  1:1639183/1‑150 (MQ=255)
                                                                    aaTTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTCCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGtt                                                                                 <  2:1598413/150‑1 (MQ=255)
                                                                       tGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCa                                                                              >  2:1041822/1‑150 (MQ=255)
                                                                        gCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCAt                                                                             >  2:303990/1‑150 (MQ=255)
                                                                        gCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCAt                                                                             >  1:1381498/1‑150 (MQ=255)
                                                                           gAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATccc                                                                          <  2:466399/150‑1 (MQ=255)
                                                                            agcagcTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCt                                                                         <  2:1330168/150‑1 (MQ=255)
                                                                             gcagcTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTa                                                                        <  1:1099121/150‑1 (MQ=255)
                                                                              cagcTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTAt                                                                       >  2:400038/1‑150 (MQ=255)
                                                                                      aCGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGc                                                               <  1:1648246/150‑1 (MQ=255)
                                                                                       cGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCa                                                              >  1:1253725/1‑150 (MQ=255)
                                                                                       cGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCa                                                              <  2:1639183/150‑1 (MQ=255)
                                                                                         cAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCatg                                                            >  2:574428/1‑150 (MQ=255)
                                                                                         cAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCatg                                                            >  2:1003143/1‑150 (MQ=255)
                                                                                         cAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCatg                                                            >  2:405256/1‑150 (MQ=255)
                                                                                            tatCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCatgatg                                                         >  2:753394/1‑150 (MQ=255)
                                                                                             atCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACACAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGc                                                        >  1:976431/1‑150 (MQ=255)
                                                                                                aaaaCGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCgg                                                     >  2:815195/1‑150 (MQ=255)
                                                                                                aaaaCGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTAATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCgg                                                     <  2:55978/150‑1 (MQ=255)
                                                                                                 aaaCGATCCGTGCACATAAATTTAATGTGATGTCGAATCTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCTTGATGCCggg                                                    <  1:252293/150‑1 (MQ=255)
                                                                                                 aaaCGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCggg                                                    <  1:944363/150‑1 (MQ=255)
                                                                                                    cGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGttt                                                 <  1:1477787/150‑1 (MQ=255)
                                                                                                      ttCCTTGCACATAAATTTAATGTGATGTCGAAACTGTGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTAc                                               <  1:272238/149‑1 (MQ=255)
                                                                                                      aTCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTAc                                               <  1:1355028/150‑1 (MQ=255)
                                                                                                      aTCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTAc                                               >  2:1355028/1‑150 (MQ=255)
                                                                                                         cGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTc                                            >  2:426969/1‑150 (MQ=255)
                                                                                                         cGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTc                                            >  2:504418/1‑150 (MQ=255)
                                                                                                          gTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTcc                                           <  1:1001614/150‑1 (MQ=255)
                                                                                                           tGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTccc                                          >  2:1197108/1‑150 (MQ=255)
                                                                                                           ggCACATAAATTTAATGTGAGGTCGAAAGTGGGGGTCAGTTCGGACGCAGGGGTGTTGGCGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGGACTGTGGTGCATCCCTATTAATCCTCATGATGCCGGGTTTACTTccc                                          <  1:1304736/149‑1 (MQ=255)
                                                                                                                aTAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTCCCCGGCa                                     <  2:768508/150‑1 (MQ=255)
                                                                                                                aTAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCATATTAATCCGCATGATGCCGGGTTTACTTCCCCGGCa                                     >  1:921666/1‑150 (MQ=255)
                                                                                                                 tAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTCCCCGGCAg                                    <  2:724182/150‑1 (MQ=255)
                                                                                                                         aTGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTCCCCGGCAGTGCtttca                            <  1:370025/150‑1 (MQ=255)
                                                                                                                           ttgATGTAGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTCCCCGGCAGTGCtttcatt                          <  1:176658/149‑1 (MQ=255)
                                                                                                                           gtgATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTCCCCGGCAGTGCtttcatt                          >  1:415942/1‑150 (MQ=255)
                                                                                                                            tgATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTCCCCGGCAGTGCtttcattt                         <  2:1715022/150‑1 (MQ=255)
                                                                                                                            ggATGTCGAAACTGGGCGTCAGTTCTGACGCTGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCTGGTTTACTTCCCCGGCAGTGCtttcattt                         <  1:605578/149‑1 (MQ=255)
                                                                                                                              aTGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTCCCCGGCAGTGCtttcatttca                       >  2:363487/1‑150 (MQ=255)
                                                                                                                               tGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTCCCCGGCAGTGCTTTCATTTCAg                      <  2:473935/150‑1 (MQ=255)
                                                                                                                               tGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTCCCCGGCAGTGCTTTCATTTCAg                      >  2:746184/1‑150 (MQ=255)
                                                                                                                               tGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTCCCCGGCAGTGCTTTCATTTCAg                      <  2:802209/150‑1 (MQ=255)
                                                                                                                                ggcGAAACTGGGCGGCAGGTCTGACGCAGGGGTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTCCCCGGCAGTGCTTTCATTTCAGc                     <  1:113918/148‑1 (MQ=255)
                                                                                                                                gTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTGGCATCCCTATTAATCATCATGATGACGGGTTTACTTCCCCGGCAGTGAATTCATTTCAGc                     >  1:63236/1‑150 (MQ=255)
                                                                                                                                        cctgCGTCAGAACTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTCCCCGGCAGTGCTTTCATTTCAGCGTACAATc             >  2:581183/4‑150 (MQ=255)
                                                                                                                                         gggCGTCAGTTCTGACGCATGGTTGTTGGAGGCCGCAGATATTCTGTAATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTCCCCGGCAGTGCTTTCATTTCAGCGTACAATCg            <  2:62390/150‑1 (MQ=255)
                                                                                                                                         gggCGTCAGTTCTGACGCAGGGGGGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCTTGATGCCGGGTTTACTTCCCCGGCAGTGCTTTCATTTCAGCGTACAATCg            <  2:779773/150‑1 (MQ=255)
                                                                                                                                         ctgCGTCAGAACTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTCCCCGGCAGTGCTTTCATTTCAGCGTACAATCg            >  1:190255/3‑150 (MQ=255)
                                                                                                                                          tgCGTCAGAACTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTCCCCGGCAGTGCTTTCATTTCAGCGTACAATCGc           >  2:649029/2‑150 (MQ=255)
                                                                                                                                             gTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTCCCCGGCAGTGCTTTCATTTCAGCGTACAATCGCcac        <  1:290294/150‑1 (MQ=255)
                                                                                                                                                aGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTCCCCGGCAGTGCTTTCATTTCAGCGTACAATCGCCACAtt     >  2:929375/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTCCCCGGCAGTGCTTTCATTTCAGCGTACAATCGCCACATTGCT  >  NC_000913/4604607‑4604903

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