Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,193,565 T→C H70R (CAC→CGC)  mnmA ← tRNA‑specific 2‑thiouridylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,193,5650TC100.0% 965.9 / NA 238H70R (CAC→CGC) mnmAtRNA‑specific 2‑thiouridylase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (129/109);  total (129/109)

Only 100 of 244 total aligned reads displayed.

TGCACAGAATATCCGGATTCGGCGTGCGACCGGCTTTATATTCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGTGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAGATCCGCTGCCGCTGTGCAATATTCCTCACCGTCGTCTTCTTCCCAGTTCTT  >  NC_000913/1193467‑1193662
                                                                                                  |                                                                                                 
tgCACAGAATATCCGGATTCGGCGTGCGACCGGCTTTATATTCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTgcgc                                                                                                 >  1:1558434/1‑101 (MQ=255)
tgCACAGAATATCCGGATTCGGCGTGCGACCGGCTTTATATTCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTgcgc                                                                                                 <  1:1558025/101‑1 (MQ=255)
 gCACAGAATATCCGGATTCGGCGTGCGACCGGCTTTATATTCGGCAAGAAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCa                                                                                                <  1:1558083/101‑1 (MQ=255)
  cacaGAATATCCGGATTCGGCGTGCGACCGGCTTTATATTCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAg                                                                                               >  1:1558440/1‑101 (MQ=255)
  cacaGAATATCCGGATTCGGCGTGCGACCGGCTTTATATTCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAg                                                                                               >  2:1558439/1‑101 (MQ=255)
   acaGAATATCCGGATTCGGCGTGCGACCGGCTTTATATTCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGt                                                                                              >  1:1558443/1‑101 (MQ=255)
    caGAATATCCGGATTCGGCGTGCGACCGGCTTTATATTCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGtt                                                                                             >  1:1558445/1‑101 (MQ=255)
    caGAATATCCGGATTCGGCGTGCGACCGGCTTTATATTCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGtt                                                                                             >  2:1558446/1‑101 (MQ=255)
      gAATATCCGGATTCGGCGTGCGACCGGCTTTATATTCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCa                                                                                           <  2:1557983/101‑1 (MQ=255)
      gAATATCCGGATTCGGCGTGCGACCGGCTTTATATTCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCa                                                                                           <  1:1558148/101‑1 (MQ=255)
          atCCGGATTCGGCGTGCGACCGGCTTTATATTCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGc                                                                                       <  1:1557811/101‑1 (MQ=255)
          atCCGGATTCGGCGTGCGACCGGCTTTATATTCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGc                                                                                       >  2:1558452/1‑101 (MQ=255)
          atCCGGATTCGGCGTGCGACCGGCTTTATATTCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGc                                                                                       >  1:1558451/1‑101 (MQ=255)
          atCCGGATTCGGCGTGCGACCGGCTTTATATTCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGc                                                                                       <  2:1558087/101‑1 (MQ=255)
           tCCGGATTCGGCGTGCGACCGGCTTTATATTCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGcc                                                                                      >  1:1558454/1‑101 (MQ=255)
           tCCGGATTCGGCGTGCGACCGGCTTTATATTCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGcc                                                                                      <  1:1557943/101‑1 (MQ=255)
             cGGATTCGGCGTGCGACCGGCTTTATATTCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGa                                                                                    <  1:1558090/101‑1 (MQ=255)
             cGGATTCGGCGTGCGACCGGCTTTATATTCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGa                                                                                    <  1:1558038/101‑1 (MQ=255)
              ggATTCGGCGTGCGACCGGCTTTATATTCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAg                                                                                   >  2:1558459/1‑101 (MQ=255)
              ggATTCGGCGTGCGACCGGCTTTATATTCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAg                                                                                   <  2:1558031/101‑1 (MQ=255)
              ggATTCGGCGTGCGACCGGCTTTATATTCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAg                                                                                   >  1:1558461/1‑101 (MQ=255)
              ggATTCGGCGTGCGACCGGCTTTATATTCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAg                                                                                   >  1:1558458/1‑101 (MQ=255)
                aTTCGGCGTGCGACCGGCTTTATATTCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCt                                                                                 >  1:1558465/1‑101 (MQ=255)
                 ttCGGCGTGCGACCGGCTTTATATTCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCtt                                                                                >  2:1558466/1‑101 (MQ=255)
                  tCGGCGTGCGACCGGCTTTATATTCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCa                                                                                           <  2:1558055/89‑1 (MQ=255)
                   cGGCGTGCGACCGGCTTTATATTCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGt                                                                              <  1:1557977/101‑1 (MQ=255)
                    ggCGTGCGACCGGCTTTATATTCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTc                                                                             >  1:1558472/1‑101 (MQ=255)
                    ggCGTGCGACCGGCTTTATATTCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTc                                                                             >  2:1558471/1‑101 (MQ=255)
                      cGTGCGACCGGCTTTATATTCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGc                                                                           >  1:1558475/1‑101 (MQ=255)
                      cGTGCGACCGGCTTTATATTCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGc                                                                           >  2:1558476/1‑101 (MQ=255)
                        tGCGACCGGCTTTATATTCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAg                                                                         <  1:1558052/101‑1 (MQ=255)
                         gCGACCGGCTTTATATTCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGa                                                                        >  1:1558479/1‑101 (MQ=255)
                            aCCGGCTTTATATTCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAg                                                                     <  2:1558360/101‑1 (MQ=255)
                            aCCGGCTTTATATTCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAg                                                                     <  2:1558191/101‑1 (MQ=255)
                             ccGGCTTTATATTCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGc                                                                    >  2:1558486/1‑101 (MQ=255)
                             ccGGCTTTATATTCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGc                                                                    >  1:1558484/1‑101 (MQ=255)
                              cGGCTTTATATTCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGcc                                                                   >  2:1558490/1‑101 (MQ=255)
                              cGGCTTTATATTCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGcc                                                                   >  1:1558489/1‑101 (MQ=255)
                                 cTTTATATTCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTgg                                                                >  2:1558494/1‑101 (MQ=255)
                                 cTTTATATTCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTgg                                                                >  1:1558495/1‑101 (MQ=255)
                                   ttATATTCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGc                                                              <  2:1558170/101‑1 (MQ=255)
                                    tataTTCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCa                                                             <  2:1558045/101‑1 (MQ=255)
                                     atatTCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCAt                                                            >  1:1558504/1‑101 (MQ=255)
                                     atatTCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCAt                                                            >  1:1558500/1‑101 (MQ=255)
                                     atatTCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCAt                                                            >  2:1558498/1‑101 (MQ=255)
                                       atTCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCa                                                          >  1:1558507/1‑101 (MQ=255)
                                         tCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGc                                                        <  1:1558012/101‑1 (MQ=255)
                                          cGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGcc                                                       <  2:1558259/101‑1 (MQ=255)
                                          cGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGcc                                                       >  1:1558509/1‑101 (MQ=255)
                                           ggCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCa                                                      >  2:1558510/1‑101 (MQ=255)
                                           ggCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCa                                                      >  1:1558512/1‑101 (MQ=255)
                                           ggCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCa                                                      <  2:1557929/101‑1 (MQ=255)
                                            gCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAg                                                     <  2:1558155/101‑1 (MQ=255)
                                            gCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAg                                                     <  2:1557990/101‑1 (MQ=255)
                                             cAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAGa                                                    <  2:1558076/101‑1 (MQ=255)
                                             cAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAGa                                                    >  1:1558516/1‑101 (MQ=255)
                                             cAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAGa                                                    <  1:1558130/101‑1 (MQ=255)
                                               aGGAACAGTTCGAAGACGGTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCa                                                          <  2:1557821/93‑1 (MQ=255)
                                                ggAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAGATcc                                                 <  1:1558015/101‑1 (MQ=255)
                                                   aCAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAGATccgct                                              <  2:1558199/101‑1 (MQ=255)
                                                   aCAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAGATccgct                                              <  2:1558097/101‑1 (MQ=255)
                                                   aCAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAGATccgct                                              >  2:1558517/1‑101 (MQ=255)
                                                     aGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAGATccgctgc                                            >  1:1558520/1‑101 (MQ=255)
                                                        tCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAGATccgctgccgc                                         >  1:1558523/1‑101 (MQ=255)
                                                         cGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAGATccgctgccgct                                        <  1:1558234/101‑1 (MQ=255)
                                                          gAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAGATccgctgccgctg                                       >  1:1558527/1‑101 (MQ=255)
                                                          gAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAGATccgctgccgctg                                       <  2:1558406/101‑1 (MQ=255)
                                                            aGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAGATCCGCTGCCGCtgtg                                     <  1:1558099/101‑1 (MQ=255)
                                                             gACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAGATCCGCTGCCGCTGTGc                                    <  1:1558208/101‑1 (MQ=255)
                                                               cGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAGATCCGCTGCCGCTGTGCaa                                  >  2:1558529/1‑101 (MQ=255)
                                                                gTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAGATCCGCTGCCGCtgt                                      <  2:1558149/96‑1 (MQ=255)
                                                                gTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAGATCCGCTGCCGCTGTGCAat                                 <  1:1557819/101‑1 (MQ=255)
                                                                 ttGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAGATCCGCTGCCGCTGTGCAata                                >  1:1558533/1‑101 (MQ=255)
                                                                 ttGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAGATCCGCTGCCGCTGTGCAata                                >  1:1558531/1‑101 (MQ=255)
                                                                  tGTCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAGATCCGCTGCCGCTGTGCAatat                               <  1:1558249/101‑1 (MQ=255)
                                                                    tCCCAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAGATCCGCTGCCGCTGTGCAATATTc                             >  1:1558537/1‑101 (MQ=255)
                                                                       cAGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAGATCCGCTGCCGCTGTGCAATATTCCTc                          <  2:1558248/101‑1 (MQ=255)
                                                                        aGTACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAGATCCGCTGCCGCTGTGCAATATTCCTCa                         >  1:1558539/1‑101 (MQ=255)
                                                                          tACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAGATCCGCTGCCGCTGTGCAATATTCCTCAcc                       >  2:1558541/1‑101 (MQ=255)
                                                                          tACTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAGATCCGCTGCCGCTGTGCAATATTCCTCAcc                       >  1:1558543/1‑101 (MQ=255)
                                                                            cTCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAGATCCGCTGCCGCTGTGCAATATTCCTCACcgt                     <  1:1558201/101‑1 (MQ=255)
                                                                             tCGGCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAGATccgc                                               >  2:1558546/1‑74 (MQ=255)
                                                                               ggCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAGATCCGCTGCCGCTGTGCAATATTCCTCACcgtcgt                  <  2:1558141/101‑1 (MQ=255)
                                                                               ggCAGCAAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAGATCCGCTGCCGCTGTGCAATATTCCTCACcgtcgt                  >  1:1558549/1‑101 (MQ=255)
                                                                                gcagcaAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAGATCCGCTGCCGCTGTGCAATATTCCTCACCgtcgtc                 >  2:1558550/1‑101 (MQ=255)
                                                                                 cagcaAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAGATCCGCTGCCGCTGTGCAATATTCCTCACCGTCGtct                <  1:1558163/101‑1 (MQ=255)
                                                                                 cagcaAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAGATCCGCTGCCGCTGTGCAATATTCCTCACCGTCGtct                <  1:1558119/101‑1 (MQ=255)
                                                                                    caAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAGATCCGCTGCCGCTGTGCAATATTCCTCACCGTCGtcttct             >  1:1558553/1‑101 (MQ=255)
                                                                                    caAAGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAGATCCGCTGCCGCTGTGCAATATTCCTCACCGTCGtcttct             <  1:1558103/101‑1 (MQ=255)
                                                                                     aaaGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAGATCCGCTGCCGCTGTGCAATATTCCTCACCGTCGTcttctt            >  2:1558556/1‑101 (MQ=255)
                                                                                       aGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAGATCCGCTGCCGCTGTGCAATATTCCTCACCGTCGTCTTCTTcc          >  1:1558558/1‑101 (MQ=255)
                                                                                       aGTTAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAGATCCGCTGCCGCTGTGCAATATTCCTCACCGTCGTCTTCTTcc          >  1:1558561/1‑101 (MQ=255)
                                                                                         ttAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAGATCCGCTGCCGCTGTGCAATATTCCTCACCGTCGTCTTCTTCCCa        <  2:1558521/101‑1 (MQ=255)
                                                                                          tAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAGATCCGCTGCCGCTGTGCAATATTCCTCACCGTCGTCTTCTTCCCAg       <  1:1557870/101‑1 (MQ=255)
                                                                                          tAACGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAGATCCGCTGCCGCTGTGCAATATTCCTCACCGTCGTCTTCTTCCCAg       <  2:1558106/101‑1 (MQ=255)
                                                                                           aaCGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAGATCCGCTGCCGCTGTGCAATATTCCTCACCGTCGTCTTCTTCCCAGt      >  2:1558563/1‑101 (MQ=255)
                                                                                           aaCGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAGATCCGCTGCCGCTGTGCAATATTCCTCACCGTCGTCTTCTTCCCAGt      <  1:1558137/101‑1 (MQ=255)
                                                                                            aCGGTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAGATCCGCTGCCGCTGTGCAATATTCCTCACCGTCGTCTTCTTCCCAGtt     >  1:1558565/1‑101 (MQ=255)
                                                                                              ggTGCGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAGATCCGCTGCCGCTGTGCAATATTCCTCACCGTCGTCTTCTTCCCAGttct   <  1:1558265/101‑1 (MQ=255)
                                                                                               gtgcGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAGATCCGCTGCCGCTGTGCAATATTCCTCACCGTCGTCTTCTTCCCAGttctt  <  2:1558287/101‑1 (MQ=255)
                                                                                                  |                                                                                                 
TGCACAGAATATCCGGATTCGGCGTGCGACCGGCTTTATATTCGGCAAGGAACAGTTCGAAGACGTTGTCCCAGTACTCGGCAGCAAAGTTAACGGTGTGCAGTTCAATGCCGAGCTTGTCGCAGACAGCCTGGGCATCAGCCAGATCCGCTGCCGCTGTGCAATATTCCTCACCGTCGTCTTCTTCCCAGTTCTT  >  NC_000913/1193467‑1193662

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