Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,512,729 T→A 5.3% intergenic (‑23/‑313) mntH ← / → nupC manganese/divalent cation transporter/nucleoside (except guanosine) transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,512,7290TA 5.3% 372.2 / 7.0 126intergenic (‑23/‑313)mntH/nupCmanganese/divalent cation transporter/nucleoside (except guanosine) transporter
Reads supporting (aligned to +/- strand):  ref base T (59/59);  new base A (2/5);  total (61/65)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.75e-01

Only 100 of 128 total aligned reads displayed.

GAATATTGGTCGCAAAGTTACCGGGATCGATATAACCAATCGCCGCAATGAACGCAGGTCCCATTAATGCGAGCCTCATCTTGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGATAGACCAGGATTATGACTACGAAAAGATTGATGACTTGAATGTTATGGTTATGTTTAATGTTAGCACATTTAC  >  NC_000913/2512580‑2512876
                                                                                                                                                     |                                                                                                                                                   
gAATATTGGTCGCAAAGTTACCGGGATCGATATAACCAATCGCCGCAATGAACGCAGGTCCCATTAATGCGAGCCTCATCTTGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCt                                                                                                                                                     >  2:1365562/1‑150 (MQ=255)
gAATATTGGTCGCAAAGTTACCGGGATCGATATAACCAATCGCCGCAATGAACGCAGGTCCCATTAATGCGAGCCTCATCTTGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCt                                                                                                                                                     >  1:1751600/1‑150 (MQ=255)
     ttGGTCGCAAAGTTACCGGGATCGATAAAACCAATCGCCGCAATGAACGCAGGTCCCATTAATGCGAGCCTCATCTTGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCt                                                                                                                                                <  2:833394/150‑1 (MQ=255)
       ggTCGCAAAGTTACCGGGATCGATATAACCAATCGCCGCAATGAACGCAGGTCCCATTAATGCGAGCCTCATCTTGAGCGCAGCCCGTCCGCAGCAACTCTCAAATCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTAt                                                                                                                                              >  1:1876074/1‑150 (MQ=255)
        gTCGCAAAGTTACCGGGATCGATATAACCAATCGCCGCAATGAACGCAGGTCCCATTAATGCGAGCCTCATCTTGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATg                                                                                                                                             <  2:588302/150‑1 (MQ=255)
         tCGCAAAGTTACCGGGATCGATATAACCAATCGCCGCAATGAACGCAGGTCCCATTAATGCGAGCCTCATCTTGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGt                                                                                                                                            >  2:1090529/1‑150 (MQ=255)
          cGCAAAGTTACCGGTATCGATATAACCAATCGCCGCAATTAACGCAGGTACGATTAATGAGAGCCTCATCTTGCGCGCCTCCCGTCAGCTGCTACTATCAACGCTATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGtt                                                                                                                                           <  2:1335764/150‑1 (MQ=255)
          cGCAAAGTTACCGGGATCGATATAACCAATCGCCGCAATGAACGCAGGTCCCATTAATGCGAGCCTCATCTTGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGtt                                                                                                                                           <  1:137975/150‑1 (MQ=255)
          cGCAAAGTTACCGGGATCGATATAACAACTCGCCGCAATGAACGCCGGTCCCATTAATGCGAGACTCATGTTGCGCGGCGCCAGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGtt                                                                                                                                           <  2:939817/150‑1 (MQ=255)
           gCAAAGTTACCGGGATCGGTATAACCAATCGCCGCAATGAACGCAGGTCCCATTAATGCGAGCCTCATCTTGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGttt                                                                                                                                          <  1:586049/150‑1 (MQ=255)
            cAAAGTTACCGGGATCGATATAACCAATCGCCGCAATGAACGCAGGTCCCATTAATGCGAGCCTCATCTTGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCAAAGGCTATGtttt                                                                                                                                         <  2:486849/150‑2 (MQ=255)
            cAAAGTTACCGGGATCGATATAACCAATCGCCGCAATGAACGCAGGTCCCATTAATGCGAGCCTCATCTTGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCAAAGGCTATGtttt                                                                                                                                         <  2:697500/150‑2 (MQ=255)
            cAAAGTTACCGGGATCGATATAACCAATCGACGCAATGAACGCAGGTCCCATTAATTCGAGCCTCATCTTGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTc                                                                                                                                         <  1:534422/150‑1 (MQ=255)
              aaGTTACCGGGATCGATATAACCAATCGCCGCAATGAACGCAGGTCCCATTAATGCGAGCCTCATCTTGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCAt                                                                                                                                       <  1:380722/150‑1 (MQ=255)
               aGTTACCGGGATCGATATAACCAATCGCCGCAATGAACGCAGGTCCCATTAATGCGAGCCTCATCTTGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATg                                                                                                                                      <  1:1495376/150‑1 (MQ=255)
               aGTTAACGGGATCGATATAACCAATCGCCGCAATGAACGCAGGTCCCATTAATGCGAGCCTCATCTTGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATg                                                                                                                                      <  2:940596/150‑1 (MQ=255)
                 ttACCGGGATCGATATAACCAATCGCCGCAATGAACGCAGGTCCCATTAATGCGAGCCTCATCTTGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCAAAGGCTATGttttagagg                                                                                                                                    <  1:750566/150‑7 (MQ=255)
                    ccGGGATCGATATAACCAATCGCCGCAATGAACGCAGGTCCCATTAATGCGAGCCTCATCTTGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATg                                                                                                                                 >  1:1134069/1‑150 (MQ=255)
                    ccGGGATCGATATAACCAATCGCCGCAATGAACGCAGGTCCCATTAATGCGAGCCTCATCTTGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGGTTCATGCTATg                                                                                                                                 >  1:491185/1‑150 (MQ=255)
                    ccGGGATCGATATAACCAATCGCCGCAATGAACGCAGGTCCCATTAATGCGAGCCTCATCTTGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTCATCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATg                                                                                                                                 >  2:1737041/1‑150 (MQ=255)
                       ggATCGATATAACCAATCTCCGCAATGAACGCAGGTCCCATAAATGCGAGCCTCATCTTGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCa                                                                                                                              <  1:94710/150‑1 (MQ=255)
                       ggATCGATATAACCAATCGCCGCAATGAACGCAGGTCCCATTAATGCGAGCCTCATCTTGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCa                                                                                                                              >  1:115144/1‑150 (MQ=255)
                          tCGATATAACCAATCGCCGCAATGAACGCAGGTCCCATTAATGCGAGCCTCATCTTGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAAc                                                                                                                           >  1:1864779/1‑150 (MQ=255)
                          tCGATATAACCAATCGCCGCAATGAACGCAGGTCCCATTAATGCGAGCCTCATCTTGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAAc                                                                                                                           >  2:1057332/1‑150 (MQ=255)
                            gATATAACCAATCGCCGCAATGAACGCAGGTCCCATTAATGCGAGCCTCATCTTGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACga                                                                                                                         >  1:1411596/1‑150 (MQ=255)
                             atatAACCAATCGCCGCAATGAACGCAGGTCCCATTAATGCGAGCCTCATCTTGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACgag                                                                                                                        >  2:780198/1‑150 (MQ=255)
                             atatAACCAATCGCCGCAATGAACGCAGGTCCCATTAATGCGAGCCTCATCTTGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACgag                                                                                                                        <  2:124943/150‑1 (MQ=255)
                                taACCAATCGCCGCAATGAACGCAGGTCCCATTAATGCGAGCCTCATCTTGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTCCCTCTAAAACATAGACTTTGCTATGTTTCATGCTATGCCCAACGAGAAt                                                                                                                     >  2:534219/1‑150 (MQ=255)
                                    cAATCGCCGCAATGAACGCAGGTCCCATTAATGCGAGCCTCATCTTGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGAtt                                                                                                                 >  1:1075408/1‑150 (MQ=255)
                                    cAATCGCCGCAATGAACGCAGGTCCCATTAATGCGAGCCTCATCTTGCGCGCCACCCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGACAATGAtt                                                                                                                 >  1:848813/1‑150 (MQ=255)
                                      aTCGCCGCAATGAACGCAGGTCCCATTAATGCGAGCCTCATCTTGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAGAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTAt                                                                                                               <  1:877746/150‑1 (MQ=255)
                                          ccgcAATGAACGCAGGTCCCATTAATGCGAGCCTCATCTTGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCaaa                                                                                                           >  1:515804/1‑150 (MQ=255)
                                          ccgcAATGAACGCAGGTCCCATTAATGCGAGCCTCATCTTGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCaaa                                                                                                           >  2:138870/1‑150 (MQ=255)
                                           ctcactGAACGCAGGTCCCATTAATGCGAGCCTCATCTTGAGCGCAGCCCGTCCGCTGATACTATCAACTCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAAt                                                                                                          <  1:1082903/145‑1 (MQ=255)
                                            gcAATGAACGCAGGTCCCATTAATGCGAGCCTCCTCTTGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCACAtt                                                                                                         >  1:124266/1‑150 (MQ=255)
                                             cAATGAACGCAGGTCCCATTAATGCGAGCCTCATCTTGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTc                                                                                                        >  2:476765/1‑150 (MQ=255)
                                              aaTGAACGCAGGTCCCATTAATGCGAGCCTCATCTTGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCa                                                                                                       >  2:1772694/1‑150 (MQ=255)
                                                  aaCGAAGGTCCCATTAATGCGAGCCTCATATTGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAAGTCAATTa                                                                                                   >  2:518714/1‑150 (MQ=255)
                                                     gCAGGTCCCATTAATGCGAGCCTCATCTTGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAAt                                                                                                <  1:1707644/150‑1 (MQ=255)
                                                           cccATTAATGCGAGCCTCATCTTGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATtg                                                                                          <  2:781777/150‑1 (MQ=255)
                                                           cccATTAATGCGAGCCTCATCTTGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTAGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCAATAAAATGGATtg                                                                                          >  2:474581/1‑150 (MQ=255)
                                                               ttAATGCGAGCCTCATCTTGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTgtggt                                                                                      <  2:852946/150‑1 (MQ=255)
                                                                   tGCGAGCCTCATCTTGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGAtt                                                                                  <  2:791736/150‑1 (MQ=255)
                                                                         ccTCATCTTGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTg                                                                            <  1:138870/150‑1 (MQ=255)
                                                                         ccTCATCTTGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTg                                                                            <  2:572291/150‑1 (MQ=255)
                                                                          cTCATCTTGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCCTCTTGTGCCTCTAAAACATCGCCTTTGCTATGGTTCCTGCTATGCCAAACGAGAATGATTATCAAATTCATTTTAATGGATTGTGGTGATTTCTCTGa                                                                           >  2:278797/1‑150 (MQ=255)
                                                                          cTCATCTTGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGa                                                                           <  2:1134069/150‑1 (MQ=255)
                                                                           tCATCTTGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGAt                                                                          <  2:669015/150‑1 (MQ=255)
                                                                              tCTTGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGATAGa                                                                       >  2:1321957/1‑150 (MQ=255)
                                                                                ttGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGATAGAcc                                                                     <  2:393535/150‑1 (MQ=255)
                                                                                 tGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGATATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGATAGACCa                                                                    >  2:1102970/1‑150 (MQ=255)
                                                                                      gccgccCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGATAGACCAGGAtt                                                               <  1:3257/150‑1 (MQ=255)
                                                                                         gccCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGATAGACCAGGATTATg                                                            <  2:49583/150‑1 (MQ=255)
                                                                                         gccCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGATAGACCAGGATTATg                                                            <  2:734912/150‑1 (MQ=255)
                                                                                             gTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGATAGACCAGGATTATGACTa                                                        <  2:119335/150‑1 (MQ=255)
                                                                                                       aCTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACAAAGCCTTAGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGATAGACCAGGATTATGACTACGAAAAGAtt                                              <  2:639167/150‑1 (MQ=255)
                                                                                                         tctcAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGATAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGATAGACCAGGATTATGACTACGAAAAGATtga                                            <  2:1276157/150‑1 (MQ=255)
                                                                                                            cAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGATAGACCAGGATTATGACTACGAAAATATtgatga                                         >  1:1236511/1‑150 (MQ=255)
                                                                                                               cgcgATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAATTGGATTGTGGTGATTTCTCTGATAGACCAGGATTATGACTACGAAAATATTGATGACtt                                      >  1:1832761/1‑150 (MQ=255)
                                                                                                               cgcgATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGATAGACCAGGATTATGACTACGAAAAGATTGATGACtt                                      <  1:499674/150‑1 (MQ=255)
                                                                                                               cgcgATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGATAGACCAGGATTATGACTACGAAAAGATTGATGACtt                                      >  1:1791350/1‑150 (MQ=255)
                                                                                                                    tAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGATAGACCAGGATTATGACTACGAAAAGATTGATGACTTGAATg                                 <  2:1459732/150‑1 (MQ=255)
                                                                                                                    tAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGATAGACCAGGATTATGACTACGAAAAGATTGATGACTTGAATg                                 <  1:386816/150‑1 (MQ=255)
                                                                                                                     aGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGATAGACCAGGATTATGACTACGAAAAGATTGATGACTTGAAttt                                >  2:829038/1‑148 (MQ=255)
                                                                                                                     aGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATCTAAATGGATTGTGGTGATTTCTCTGATAGACCAGGATTATGACTACGAAAAGATTGATGACTTGAATgt                                >  2:694003/1‑150 (MQ=255)
                                                                                                                      gTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGATAGACCAGGATAATGACTACGAAAAGATTGATGACTTGAATgtt                               >  2:1425031/1‑150 (MQ=255)
                                                                                                                          ttCATCTTGTGCCTCTAAAAAATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGTTTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGATAGACCAGGATTATGACTACGAAACGATTGATGACTTGAATgttatgg                           <  2:1360450/149‑1 (MQ=255)
                                                                                                                          gTCATCTTGTGCCTCTATAACATAGCCTTTGCTATGTTTCATGCTATGCCAATCGAGAATGATTATCTAATTCATTTAATTGGATTGTGGTGATTTCTCTGATAGACCAGGATTATGACTACGAAAAGATTGATGACTTGAATgttatgg                           <  1:474581/150‑1 (MQ=255)
                                                                                                                          gTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGATAGACCAGGATTATGACTACGAAAAGATTGATGACTTGAATgttatgg                           >  1:1104237/1‑150 (MQ=255)
                                                                                                                          gTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGATAGACCAGGATTATGACTACGAAAAGATTGATGACTTGAATgttatgg                           <  2:1751600/150‑1 (MQ=255)
                                                                                                                            cATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGATAGACCAGGATTATGACTACGAAAAGATTGATGACTTGAATgttatggtt                         >  2:960365/1‑150 (MQ=255)
                                                                                                                            cATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGATAGACCAGGATTATGACTACGAAAAGATTGATGACTTGAATgttatggtt                         >  2:1790578/1‑150 (MQ=255)
                                                                                                                             aTCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGATAGACCAGGATTATGACTACGAAAAGATTGATGACTTGATTgttatggtta                        >  1:432128/1‑150 (MQ=255)
                                                                                                                               ctagtgCCTCTAAAACATAGCCGTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGATAGACCAGGATTATGACTACGAAAAGATTGATGACTTGAATgttatggttatg                      <  1:91237/147‑1 (MQ=255)
                                                                                                                                 tgtgCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGATAGACCAGGATTATGACTACGAAAAGATTGATGACTTGAATGTTATGGTTATGtt                    >  1:1309510/1‑150 (MQ=255)
                                                                                                                                  gttccTCTTATACATATCATTTTCTATGTTTCATGCTATGCCATACGAGTATGATAATCAAATTCATTTAAAAGGATTGTGGTGATTGCTTTGATAGACCAGGATTATGACTACGAAAATATTGATGAATTGAATGTTATGGTTATGttt                   <  2:1733702/147‑1 (MQ=255)
                                                                                                                                   tgCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGATAGACCAGGATTATGACTACGAAAAGATTGATGACTTGAATGTTATGGTTATGTTTa                  >  2:1235733/1‑150 (MQ=255)
                                                                                                                                    gCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGATAGACCAGGATTATGACTACGAAAAGTTTGATGACTTGAATGTTATGGTTATGTTTaa                 >  1:965285/1‑150 (MQ=255)
                                                                                                                                     ccTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGATAGACCAGGATTATGACTACGAAAAGATTGATGACTTGAATGTTATGGTTATGTTTAAt                >  2:349004/1‑150 (MQ=255)
                                                                                                                                      gcatgaaaCATAGCAAAGGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGATAGACCAGGATTATGACTACGAAAAGATTGATGACTTGAATGTTATGGTTATGTTTAATg               >  1:1328900/6‑150 (MQ=255)
                                                                                                                                      ctctAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATTATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGATAGACCAGGATTATGACTACGAAAAGATTGATGACTTGAATGTTATGGGTATGTTTAATg               <  1:518714/150‑1 (MQ=255)
                                                                                                                                      ctctAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTGATTTAAATGGATTGTGGTGATTTCTCTGATAGACCAGGATTATGACTACGAAAAGATTGATGACTTGAATGTTATGGTTATGTTTAATg               <  2:4645/150‑1 (MQ=255)
                                                                                                                                      ctctAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGATAGACCAGGATTATGACTACGAAAAGATTGATGACTTGAATGTTATGGTTATGTTTAATg               >  1:1197086/1‑150 (MQ=255)
                                                                                                                                      ctctAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGATAGACCAGGATTATGACTACGAAAAGATTGATGACTTGAATGTTATGGTTATGTTTAATg               <  2:491185/150‑1 (MQ=255)
                                                                                                                                       tctAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGATAGACCAGGATTATGACTACGAAAAGATTGATGACTTGAATGTTATGGTTATGTTTAATGt              >  2:1610066/1‑150 (MQ=255)
                                                                                                                                        ctAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGATAGACCAGGATTATGACTACGAAAAGATTGATGACTTGAATGTTATGGTTATGTTTAATGtt             >  2:1048917/1‑150 (MQ=255)
                                                                                                                                        ctAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGATAGACAAGGATTATGACTACGAAAAGATTGATGACTTTATTGTTATGGTTATGTTTATTGtt             >  1:1245644/1‑150 (MQ=255)
                                                                                                                                        atgaaaCATAGCAAAGGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGATAGACCAGGATTATGACTACGAAAAGATTGATGACTTGAATGTTATGGTTATGTTTAATGtt             >  1:437299/4‑150 (MQ=255)
                                                                                                                                         tAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGATAGACCAGGATTATGACTACGAAAAGATTGATGACTTGAATGTTATGGTTATGTTTAATGTTa            <  2:115144/150‑1 (MQ=255)
                                                                                                                                          aaaaCATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGATAGACCAGGATTATGACTACGAAAAGATTGATGACTTGAATGTTATGGTTATGTTTAATGTTAg           <  1:938537/150‑1 (MQ=255)
                                                                                                                                          aaaaCATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATAATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGATAGACCAGGATTATGACTACGAAAAGATTGATGACTTGAATGTTATGGTTATGTTTAATGTTAg           <  2:507743/150‑1 (MQ=255)
                                                                                                                                           aaaaATAGCTTTTGCTATGTTTCATGGTATTCAAAACGAGAATTATTATCAAATTCATTTTAATGGATTGTGGTGTTTTCTCTTATAGACCAGGATTATTACTACGATAAGATTGATGACTTGAATGTTATGGTTATGTTTAATGTTAGc          <  1:893677/150‑1 (MQ=255)
                                                                                                                                           aaaCATTGCCTTTGCTATGTTTCATGCTATGACAATCGAGAATGATTATCAAATTCATTTATATGGATTGTGGTGATTTCTCTGATAGACCAGGATTATGACTTAGAATAGATTGATGACTTGAATGTTATGGTTATGTTTAATGTTAGc          <  1:935169/150‑1 (MQ=255)
                                                                                                                                            aaCATAGCCTTTGCTATGTTTCATGATATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGATAGACCAGGATTATGAGTACGAGAAGATTGATGACTTGAATGTTATGGTTATGTTTAATGTTAGca         <  1:534219/150‑1 (MQ=255)
                                                                                                                                               aTAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGATAGACCAGGATTATGACTACGAAAAGATTGATGACTTGAATGTTATGGTTATGTTTAATGTTAGCACAt      >  1:427732/1‑150 (MQ=255)
                                                                                                                                                 aGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGATAGACCAGGATTATGACTACGAAAAGATTGATGACTTGAATGTTATGGTTATGTTTAATGTTAGCACAttt    <  1:1693796/150‑1 (MQ=255)
                                                                                                                                                 aGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGATAGACCAGGATTATGACTACGAAAAGATTGATGACTTGAATGTTATGGTTATGTTTAATGTTAGCACAttt    >  2:802213/1‑150 (MQ=255)
                                                                                                                                                  gCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGATAGACCAGGATTATGACTACGAAAAGATTGATGACTTGAATGTTATGGTTATGTTTAATGTTAGCACATTTa   >  2:1691870/1‑150 (MQ=255)
                                                                                                                                                  gCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGATAGACCAGGATTATGACTACGAAAAGATTGATGACTTGAATGTTATGGTTATGTTTAATGTTAGCACATTTa   >  1:615622/1‑150 (MQ=255)
                                                                                                                                                   ccTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGATAGACCAGGATTATGACTACGAAAAGATTGATGACTTGAATGTTATGGTTATGTTTAATGTTAGCACATTTaa  >  1:693074/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GAATATTGGTCGCAAAGTTACCGGGATCGATATAACCAATCGCCGCAATGAACGCAGGTCCCATTAATGCGAGCCTCATCTTGCGCGCCGCCCGTCCGCTGCTACTCTCAACGCGATAGTTCGTCATCTTGTGCCTCTAAAACATAGCCTTTGCTATGTTTCATGCTATGCCAAACGAGAATGATTATCAAATTCATTTAAATGGATTGTGGTGATTTCTCTGATAGACCAGGATTATGACTACGAAAAGATTGATGACTTGAATGTTATGGTTATGTTTAATGTTAGCACATTTAC  >  NC_000913/2512580‑2512876

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