Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,967,132 A→G H239R (CAC→CGC)  rho → transcription termination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,967,1320AG99.2% 447.9 / ‑4.0 123H239R (CAC→CGC) rhotranscription termination factor
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (67/55);  total (68/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

Only 100 of 126 total aligned reads displayed.

ACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTG  >  NC_000913/3966985‑3967278
                                                                                                                                                   |                                                                                                                                                  
aCATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGt                                                                                                                                                 <  1:770359/151‑1 (MQ=255)
  aTTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTc                                                                                                                                               >  3:58495/1‑151 (MQ=255)
    tGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGATGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCa                                                                                                                                              <  4:210447/150‑1 (MQ=255)
    tGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTc                                                                                                                                               >  2:1166245/1‑149 (MQ=255)
       tCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGtt                                                                                                                                          >  1:1003138/1‑151 (MQ=255)
        cAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTg                                                                                                                                         <  1:1013679/151‑1 (MQ=255)
         agagCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGc                                                                                                                                        >  2:882389/1‑151 (MQ=255)
          gagCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCg                                                                                                                                       >  1:826204/1‑151 (MQ=255)
             cATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGa                                                                                                                                     >  4:282192/1‑150 (MQ=255)
              aTTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGaa                                                                                                                                    >  4:425075/1‑150 (MQ=255)
                tGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAAt                                                                                                                                  <  2:495607/150‑1 (MQ=255)
                    tACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTg                                                                                                                              >  2:621788/1‑150 (MQ=255)
                    tACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTg                                                                                                                              >  1:1207163/1‑150 (MQ=255)
                     aCAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGAt                                                                                                                            <  3:84257/151‑1 (MQ=255)
                      cAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATc                                                                                                                           <  1:934865/151‑1 (MQ=255)
                        accaccCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCga                                                                                                                         <  2:221312/151‑1 (MQ=255)
                            cccGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCgaga                                                                                                                       >  2:986701/1‑149 (MQ=255)
                             ccGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAgg                                                                                                                    >  1:161193/1‑151 (MQ=255)
                             ccGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAg                                                                                                                     <  1:261310/150‑1 (MQ=255)
                              cGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGc                                                                                                                   <  1:730873/151‑1 (MQ=255)
                               ggATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGcc                                                                                                                  >  3:143795/1‑151 (MQ=255)
                                 aTTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCaa                                                                                                                <  1:92856/151‑1 (MQ=255)
                                 aTTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCaa                                                                                                                <  1:1198217/151‑1 (MQ=255)
                                    gtgtgCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAAc                                                                                                              <  1:1099424/150‑1 (MQ=255)
                                    gtgtgCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAAc                                                                                                              <  3:425069/150‑1 (MQ=255)
                                    gtgtgCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACg                                                                                                             >  2:947081/1‑151 (MQ=255)
                                        gCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCt                                                                                                          <  1:833341/150‑1 (MQ=255)
                                          tGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGt                                                                                                       >  1:301789/1‑151 (MQ=255)
                                           gATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGtt                                                                                                      <  2:760468/151‑1 (MQ=255)
                                           gATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGtt                                                                                                      >  1:626645/1‑151 (MQ=255)
                                             tGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTg                                                                                                     >  2:167539/1‑150 (MQ=255)
                                               gTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGc                                                                                                  <  1:361461/151‑1 (MQ=255)
                                                 tCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGca                                                                                                 <  2:348705/150‑1 (MQ=255)
                                                   tgctgATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGcaca                                                                                               >  2:1249498/1‑150 (MQ=255)
                                                   tgctgATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGcaca                                                                                               >  2:1190572/1‑150 (MQ=255)
                                                    gctgATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGcac                                                                                                <  4:389276/148‑1 (MQ=255)
                                                     ctgATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGcaca                                                                                               >  1:10277/1‑148 (MQ=255)
                                                      tgATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaaga                                                                                            >  3:84194/1‑150 (MQ=255)
                                                      tgATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaaga                                                                                            >  1:471573/1‑150 (MQ=255)
                                                      tgATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaaga                                                                                            >  2:1176200/1‑150 (MQ=255)
                                                         tCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaagaa                                                                                           >  2:1259471/1‑148 (MQ=255)
                                                          cgacgaACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaagaaaga                                                                                        <  2:1129025/150‑1 (MQ=255)
                                                          cgacgaACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGAc                                                                                       >  2:659930/1‑151 (MQ=255)
                                                           gacgaACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACg                                                                                      <  1:1194982/151‑1 (MQ=255)
                                                              gaACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTa                                                                                   <  3:374497/151‑1 (MQ=255)
                                                                 cGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATc                                                                                 >  1:46832/1‑150 (MQ=255)
                                                                 cGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATc                                                                                 >  4:149323/1‑150 (MQ=255)
                                                                     cGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCt                                                                            >  3:412513/1‑151 (MQ=255)
                                                                     cGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCt                                                                            <  1:1098500/151‑1 (MQ=255)
                                                                      ggAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTg                                                                           <  3:97748/151‑1 (MQ=255)
                                                                        aagaagTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCt                                                                         <  3:273396/151‑1 (MQ=255)
                                                                        aagaagTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCt                                                                         >  2:503208/1‑151 (MQ=255)
                                                                         agaagTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTc                                                                        >  2:126149/1‑151 (MQ=255)
                                                                           aagTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCg                                                                       >  1:605417/1‑150 (MQ=255)
                                                                            agTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGAc                                                                     >  2:1068146/1‑151 (MQ=255)
                                                                              tAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTc                                                                   >  1:638824/1‑151 (MQ=255)
                                                                               aaCCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTcc                                                                  >  3:86421/1‑151 (MQ=255)
                                                                               aaCCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTcc                                                                  <  1:1226638/151‑1 (MQ=255)
                                                                               aaCCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTca                                                                  <  2:822779/151‑2 (MQ=255)
                                                                                  tgagaTGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATc                                                               >  2:876202/2‑151 (MQ=255)
                                                                                  cGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATc                                                               >  2:778140/1‑151 (MQ=255)
                                                                                    agaTGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCAc                                                             <  4:10441/151‑1 (MQ=255)
                                                                                     gaTGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACt                                                            <  4:83605/151‑1 (MQ=255)
                                                                                       tGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCg                                                          >  4:202586/1‑151 (MQ=255)
                                                                                        gCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCg                                                          <  2:36285/150‑1 (MQ=255)
                                                                                               cTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTgg                                                     >  2:293255/1‑148 (MQ=255)
                                                                                                 ggTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGgcgcg                                                 <  3:415025/150‑1 (MQ=255)
                                                                                                  gTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGgcgcgcg                                               <  1:473418/151‑1 (MQ=255)
                                                                                                  gTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGgcgcgcg                                               <  1:950903/151‑1 (MQ=255)
                                                                                                     aaaGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCtt                                            >  2:979426/1‑151 (MQ=255)
                                                                                                         gTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTaca                                         >  2:104988/1‑150 (MQ=255)
                                                                                                           gAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTacaaca                                      <  4:22188/151‑1 (MQ=255)
                                                                                                            aaGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAacac                                     >  3:334229/1‑151 (MQ=255)
                                                                                                              gttgttGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCg                                   <  1:33625/151‑1 (MQ=255)
                                                                                                                   tgCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCgttgtt                              >  2:1066721/1‑151 (MQ=255)
                                                                                                                     cTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTcc                            >  3:280175/1‑151 (MQ=255)
                                                                                                                      ttCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTTCTCGACTCCATCACTCGTCTTGCGCGCGCTTACAACACCGTTGTTcc                            >  4:427522/1‑150 (MQ=255)
                                                                                                                        cTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGc                         >  1:907236/1‑151 (MQ=255)
                                                                                                                          aCCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCgg                          <  2:932278/148‑1 (MQ=255)
                                                                                                                           ccTTTGGCGAACCCGCAGCTCGCCGCGTTTAGGTTGCGGAAATGGTGATCGAGAAGGCCAACCGCCTGGTTGAGCTCAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGGGCGCTTACAACACCGTTGTTCCggag                        <  2:484474/149‑3 (MQ=255)
                                                                                                                           ccTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGt                       <  1:545112/150‑1 (MQ=255)
                                                                                                                            cTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTc                      >  3:239183/1‑150 (MQ=255)
                                                                                                                             tttGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTc                      >  2:135331/1‑149 (MQ=255)
                                                                                                                              ttGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAgg                   >  2:710368/1‑151 (MQ=255)
                                                                                                                               tGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGt                  >  2:769001/1‑151 (MQ=255)
                                                                                                                                gACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTa                 <  1:10573/151‑1 (MQ=255)
                                                                                                                                  cGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTaa                >  1:386224/1‑150 (MQ=255)
                                                                                                                                   gAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAg              <  1:1105300/151‑1 (MQ=255)
                                                                                                                                    aaCCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAgt             <  2:1059723/151‑1 (MQ=255)
                                                                                                                                     aCCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGCGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAgtg            >  2:418822/1‑151 (MQ=255)
                                                                                                                                       ccGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGtt          >  4:22903/1‑151 (MQ=255)
                                                                                                                                       ccGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGtt          <  1:1029898/151‑1 (MQ=255)
                                                                                                                                        cGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTg         <  1:1282363/151‑1 (MQ=255)
                                                                                                                                          cATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGAc       <  4:24676/151‑1 (MQ=255)
                                                                                                                                            tctcGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCg     >  2:237671/1‑151 (MQ=255)
                                                                                                                                             ctcGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCgg    >  4:362976/1‑151 (MQ=255)
                                                                                                                                               tgCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCggtg  >  2:679255/2‑151 (MQ=255)
                                                                                                                                               cgccgcGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCggtg  <  1:327945/151‑1 (MQ=255)
                                                                                                                                               cgccgcGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCggtg  <  2:32015/151‑1 (MQ=255)
                                                                                                                                               cgccgcGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCga    <  2:879275/149‑2 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGTG  >  NC_000913/3966985‑3967278

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: