Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,185,074 A→G Q716R (CAG→CGG)  rpoC → RNA polymerase beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,185,0740AG100.0% 708.1 / NA 224Q716R (CAG→CGG) rpoCRNA polymerase beta prime subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (109/115);  total (109/115)

Only 100 of 236 total aligned reads displayed.

GGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTC  >  NC_000913/4184984‑4185163
                                                                                          |                                                                                         
gggCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCggg                                                                                         >  2:1268543/1‑93 (MQ=255)
gggCTGCGGCGAACGATCGTGTATCCAAACCGATGATGCATAACCGGCAAACTGAACCCGTGATTAACCGTGATGGTCAGGAAGAGAAGCggg                                                                                         <  2:857019/93‑1 (MQ=255)
 ggCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGt                                                                                        <  2:3908098/93‑1 (MQ=255)
  gCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGtt                                                                                       <  1:3927990/93‑1 (MQ=255)
     gcggcgAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTcc                                                                                    <  1:632256/93‑1 (MQ=255)
      cggcgAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCt                                                                                   <  1:1854390/93‑1 (MQ=255)
      cggcgAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCt                                                                                   >  2:170291/1‑93 (MQ=255)
      cggcgAACGAGCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCt                                                                                   <  2:1075435/93‑1 (MQ=255)
       ggcgAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGGTTCAtt                                                                                  >  1:4558731/1‑93 (MQ=255)
        gcgAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTc                                                                                 <  2:1790701/93‑1 (MQ=255)
        gcgAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTc                                                                                 <  2:4067793/93‑1 (MQ=255)
        gcgAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTc                                                                                 <  2:4029647/93‑1 (MQ=255)
        gcgAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTc                                                                                 <  1:3325932/93‑1 (MQ=255)
            aCGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACa                                                                             >  1:154212/1‑93 (MQ=255)
            aCGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACa                                                                             <  2:1149974/93‑1 (MQ=255)
             cGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAg                                                                            >  1:447218/1‑93 (MQ=255)
             cGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAg                                                                            >  1:4050246/1‑93 (MQ=255)
             cGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAg                                                                            <  1:2460191/93‑1 (MQ=255)
              gATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGc                                                                           <  1:1613151/93‑1 (MQ=255)
               aTCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCa                                                                          <  2:2697675/93‑1 (MQ=255)
                 cGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATc                                                                        <  1:2883280/93‑1 (MQ=255)
                  gtgtATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCt                                                                       <  2:1532032/93‑1 (MQ=255)
                  gtgtATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCt                                                                       <  1:2864480/93‑1 (MQ=255)
                  gtgtATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCt                                                                       >  2:131129/1‑93 (MQ=255)
                   tgtATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTa                                                                      >  2:1630689/1‑93 (MQ=255)
                   tgtATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTa                                                                      >  2:1762351/1‑93 (MQ=255)
                    gtATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTAc                                                                     <  1:2522380/93‑1 (MQ=255)
                    gtATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTAc                                                                     >  1:1693962/1‑93 (MQ=255)
                     tATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACa                                                                    >  2:542344/1‑93 (MQ=255)
                     tATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACa                                                                    <  1:2758584/93‑1 (MQ=255)
                     tATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACa                                                                    >  2:2977155/1‑93 (MQ=255)
                     tATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACa                                                                    <  2:3800056/93‑1 (MQ=255)
                     tATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACa                                                                    <  2:2265163/93‑1 (MQ=255)
                      aTCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACat                                                                   >  1:3727738/1‑93 (MQ=255)
                       tCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACatg                                                                  >  1:2312938/1‑93 (MQ=255)
                        ccAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACatga                                                                 <  1:1056800/93‑1 (MQ=255)
                        ccAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACatga                                                                 >  1:4122371/1‑93 (MQ=255)
                         cAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACatgat                                                                <  2:4127283/93‑1 (MQ=255)
                          aaaGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACatgatg                                                               >  2:2694796/1‑93 (MQ=255)
                          aaaGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACatgatg                                                               >  2:2746125/1‑93 (MQ=255)
                          aaaGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACatgatg                                                               >  2:601378/1‑93 (MQ=255)
                          aaaGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACatgatg                                                               >  2:4430578/1‑93 (MQ=255)
                              cGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCg                                                           <  2:1879169/93‑1 (MQ=255)
                              cGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCg                                                           <  2:4018274/93‑1 (MQ=255)
                               gatgatGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGa                                                          >  1:2107632/1‑93 (MQ=255)
                               gatgatGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGGTTCCTTCAACAGCATCTACATGATGGCCGa                                                          >  2:3640657/1‑93 (MQ=255)
                                atgatgGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGAc                                                         >  2:4211252/1‑93 (MQ=255)
                                atgatgGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGAc                                                         >  1:3227471/1‑93 (MQ=255)
                                atgatgGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGAc                                                         <  2:732518/93‑1 (MQ=255)
                                 tgatgGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACt                                                        >  2:2194469/1‑93 (MQ=255)
                                 tgatgGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACt                                                        >  2:365045/1‑93 (MQ=255)
                                  gatgGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTc                                                       <  2:2534044/93‑1 (MQ=255)
                                     ggATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCgg                                                    <  2:3444173/93‑1 (MQ=255)
                                      gATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGt                                                   <  1:929450/93‑1 (MQ=255)
                                      gATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGt                                                   >  1:2380140/1‑93 (MQ=255)
                                      gATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGt                                                   >  1:2076941/1‑93 (MQ=255)
                                       aTAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTg                                                  >  1:4506185/1‑93 (MQ=255)
                                        tAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTgc                                                 <  1:3536464/93‑1 (MQ=255)
                                        tAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTgc                                                 >  2:4109920/1‑93 (MQ=255)
                                          aCCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTgcgc                                               >  2:2760390/1‑93 (MQ=255)
                                            ctcaaaaCGGAAAGGGTGATTAACGGGGAGGGTCGGGAAGAGAAGCGGTTTTCCTTAAACAGCATCTACATGATGGCCGACTCCGGTGCGCGt                                             <  2:1222176/89‑1 (MQ=255)
                                             tGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTg                                            >  1:399856/1‑93 (MQ=255)
                                              gCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTgg                                           >  1:1262775/1‑93 (MQ=255)
                                              gCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTgg                                           >  2:2383104/1‑93 (MQ=255)
                                              gCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGCGCGCGTgg                                           >  1:50029/1‑93 (MQ=255)
                                               cAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGt                                          <  2:3033427/93‑1 (MQ=255)
                                               cAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGt                                          <  2:322767/93‑1 (MQ=255)
                                                aaaCTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGtt                                         >  1:887565/1‑93 (MQ=255)
                                                 aaCTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTc                                        <  1:719971/93‑1 (MQ=255)
                                                 aaCTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTc                                        >  1:1217699/1‑93 (MQ=255)
                                                  aCTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGCTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCt                                       <  2:65042/93‑1 (MQ=255)
                                                   cTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGCTTCCTTCAACAGCATCTACATGATAGCCGACTCCGGTGCGCGTGGTTCTg                                      >  2:88761/1‑93 (MQ=255)
                                                    tGAACCCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGc                                     <  1:149849/93‑1 (MQ=255)
                                                     gAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCg                                    >  2:2622666/1‑93 (MQ=255)
                                                     gAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGGTCTGCg                                    <  2:3903582/93‑1 (MQ=255)
                                                         ccGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGcac                                <  1:1167490/93‑1 (MQ=255)
                                                         ccGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGcac                                <  1:4327669/93‑1 (MQ=255)
                                                          cGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGcaca                               <  2:2784833/93‑1 (MQ=255)
                                                            tGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGa                             >  1:3792405/1‑93 (MQ=255)
                                                              aTTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGAtt                           <  2:2035751/93‑1 (MQ=255)
                                                                   ccGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCa                      >  2:3427209/1‑93 (MQ=255)
                                                                      tGATGGTCAGGCAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCt                   >  2:3929573/1‑93 (MQ=255)
                                                                      tGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCt                   >  1:563946/1‑93 (MQ=255)
                                                                        atggTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTg                 >  2:19343/3‑93 (MQ=255)
                                                                          ggTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCt               <  1:1518831/93‑1 (MQ=255)
                                                                          ggTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCt               >  2:751899/1‑93 (MQ=255)
                                                                           gTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTg              <  1:812737/93‑1 (MQ=255)
                                                                            tCAGGAAGAGAAGCGGTTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTgg             <  1:3491158/93‑1 (MQ=255)
                                                                            tCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTgg             >  2:1662814/1‑93 (MQ=255)
                                                                               ggAAGAGAAGCGGCTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTAt          <  2:1506613/93‑1 (MQ=255)
                                                                                gAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATg         <  2:746656/93‑1 (MQ=255)
                                                                                  agagAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCg       <  2:2343053/93‑1 (MQ=255)
                                                                                  agagAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCg       >  1:3461035/1‑93 (MQ=255)
                                                                                   gagaAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGt      <  2:1959506/93‑1 (MQ=255)
                                                                                     gaAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTgg    >  1:4238865/1‑93 (MQ=255)
                                                                                     gaAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTgg    >  1:269956/1‑93 (MQ=255)
                                                                                     gaAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTgg    >  2:1443180/1‑93 (MQ=255)
                                                                                      aaGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGt   <  2:4536998/93‑1 (MQ=255)
                                                                                       aGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGTTc  <  2:1197396/93‑1 (MQ=255)
                                                                                       aGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTc  <  1:309942/93‑1 (MQ=255)
                                                                                          |                                                                                         
GGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTC  >  NC_000913/4184984‑4185163

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