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,0740AG99.4% 1096.1 / ‑3.4 347Q716R (CAG→CGG) rpoCRNA polymerase beta prime subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (182/163);  minor base T (2/0);  total (184/163)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.55e-01

Only 100 of 361 total aligned reads displayed.

GGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTC  >  NC_000913/4184984‑4185163
                                                                                          |                                                                                         
gggCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGCAGAGAAGCggg                                                                                         >  1:6501987/1‑93 (MQ=255)
gggCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCggg                                                                                         >  1:4675803/1‑93 (MQ=255)
gggCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCggg                                                                                         >  1:1574716/1‑93 (MQ=255)
gggCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCggg                                                                                         >  2:1914688/1‑93 (MQ=255)
 ggCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGt                                                                                        <  1:6020353/93‑1 (MQ=255)
  gCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGtt                                                                                       >  1:7918939/1‑93 (MQ=255)
  gCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGtt                                                                                       <  1:4126078/93‑1 (MQ=255)
   cTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGttt                                                                                      <  1:2182021/93‑1 (MQ=255)
    tGCGGCGAACGATCGTGTATCCAGAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTc                                                                                     <  2:5634637/93‑1 (MQ=255)
    tGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTc                                                                                     >  2:6264195/1‑93 (MQ=255)
      cggcgAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCt                                                                                   <  2:2009538/93‑1 (MQ=255)
      cggcgAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCt                                                                                   <  1:3613917/93‑1 (MQ=255)
      cggcgAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACCGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCt                                                                                   <  2:4671681/93‑1 (MQ=255)
       ggcgAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCtt                                                                                  <  2:238658/93‑1 (MQ=255)
           aaCGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGAGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAAc                                                                              <  1:7027217/93‑1 (MQ=255)
             cGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAg                                                                            <  1:3691444/93‑1 (MQ=255)
              gATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGc                                                                           >  2:753795/1‑93 (MQ=255)
               aTCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCa                                                                          >  1:7690362/1‑93 (MQ=255)
               aTCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCa                                                                          <  2:7645994/93‑1 (MQ=255)
                tCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGGTTCCTTCAACAGCAt                                                                         >  2:4297335/1‑93 (MQ=255)
                 cGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATc                                                                        <  1:1943970/93‑1 (MQ=255)
                  gtgtATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCt                                                                       >  1:4328742/1‑93 (MQ=255)
                  gtgtATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCt                                                                       <  2:1417624/93‑1 (MQ=255)
                  gtgtATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCt                                                                       <  2:4054665/93‑1 (MQ=255)
                   tgtATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTa                                                                      >  2:5764435/1‑93 (MQ=255)
                   tgtATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTa                                                                      >  1:5194795/1‑93 (MQ=255)
                   tgtATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTa                                                                      >  1:6215268/1‑93 (MQ=255)
                   tgtATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTa                                                                      >  2:88797/1‑93 (MQ=255)
                    gtATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTAc                                                                     >  2:7594114/1‑93 (MQ=255)
                     tATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACa                                                                    >  2:4741676/1‑93 (MQ=255)
                     tATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACa                                                                    >  1:4597466/1‑93 (MQ=255)
                     tATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACa                                                                    <  1:1300870/93‑1 (MQ=255)
                     tATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACa                                                                    <  2:5541082/93‑1 (MQ=255)
                       tCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACatg                                                                  >  2:6450723/1‑93 (MQ=255)
                       tCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACatg                                                                  >  2:7071626/1‑93 (MQ=255)
                         cAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACatgat                                                                <  2:7764602/93‑1 (MQ=255)
                          aaaGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACatgatg                                                               <  1:1394251/93‑1 (MQ=255)
                            aGCTATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGc                                                             <  1:5697646/93‑1 (MQ=255)
                            aGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGc                                                             >  2:5091228/1‑93 (MQ=255)
                            aGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGc                                                             <  1:25648/93‑1 (MQ=255)
                             gCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGcc                                                            >  2:2337224/1‑93 (MQ=255)
                             gCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGcc                                                            >  2:7879535/1‑93 (MQ=255)
                             gCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGcc                                                            <  2:256763/93‑1 (MQ=255)
                              cGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCg                                                           >  2:1570996/1‑93 (MQ=255)
                               gatgatGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGa                                                          >  2:1657764/1‑93 (MQ=255)
                                atgatgGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGAc                                                         <  1:4616454/93‑1 (MQ=255)
                                atgatgGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGAc                                                         <  1:4063278/93‑1 (MQ=255)
                                 tgatgGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACt                                                        <  1:7061175/93‑1 (MQ=255)
                                 tgatgGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACt                                                        <  2:3643382/93‑1 (MQ=255)
                                  gatgGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTc                                                       >  2:5827394/1‑93 (MQ=255)
                                   atgGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTcc                                                      <  1:7581519/93‑1 (MQ=255)
                                     ggATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCgg                                                    <  2:5189823/93‑1 (MQ=255)
                                      gATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGt                                                   <  2:2685968/93‑1 (MQ=255)
                                        tAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTgc                                                 >  2:6565614/1‑93 (MQ=255)
                                          aCCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTgcgc                                               >  2:7336226/1‑93 (MQ=255)
                                          aCCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTgcgc                                               >  2:330225/1‑93 (MQ=255)
                                          aCCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTgcgc                                               >  1:1046086/1‑93 (MQ=255)
                                          aCCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTgcgc                                               >  1:147987/1‑93 (MQ=255)
                                            cTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGt                                             <  2:6816296/93‑1 (MQ=255)
                                             tGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTg                                            >  2:357338/1‑93 (MQ=255)
                                                 aaCTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTc                                        <  2:7363722/93‑1 (MQ=255)
                                                  aCTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCt                                       <  1:1031834/93‑1 (MQ=255)
                                                  aCTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGGTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCt                                       >  1:1222695/1‑93 (MQ=255)
                                                   cTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTg                                      >  2:2145572/1‑93 (MQ=255)
                                                    tGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGc                                     >  1:7438942/1‑93 (MQ=255)
                                                     gAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCg                                    >  1:6770759/1‑93 (MQ=255)
                                                      aaaaCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCgg                                   <  2:6018632/93‑1 (MQ=255)
                                                       aaCCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGc                                  >  1:5947377/1‑93 (MQ=255)
                                                       aaCCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGc                                  >  1:3256497/1‑93 (MQ=255)
                                                        aCCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGca                                 >  2:1091133/1‑93 (MQ=255)
                                                           gTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAg                              >  1:2507145/1‑93 (MQ=255)
                                                           gTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAg                              >  1:4811506/1‑93 (MQ=255)
                                                            tGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGa                             <  2:2288362/93‑1 (MQ=255)
                                                            tGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGa                             >  1:4731560/1‑93 (MQ=255)
                                                            tGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGa                             <  1:7267763/93‑1 (MQ=255)
                                                             gATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGGCTCCGGTGCGCGTGGTTCTGCGGCACAGAt                            <  1:2114813/93‑1 (MQ=255)
                                                              aTTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGAtt                           <  2:46053/93‑1 (MQ=255)
                                                              aTTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGAtt                           >  2:3113214/1‑93 (MQ=255)
                                                                tAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCg                         <  1:3929748/93‑1 (MQ=255)
                                                                 aaCCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGt                        <  1:2326996/93‑1 (MQ=255)
                                                                  aCCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTc                       >  2:391877/1‑93 (MQ=255)
                                                                  aCCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTc                       >  1:2834014/1‑93 (MQ=255)
                                                                       gatggTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCtt                  <  1:2374830/90‑1 (MQ=255)
                                                                        atggTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTg                 >  1:5361059/3‑93 (MQ=255)
                                                                         tggTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGc                <  2:6880735/92‑1 (MQ=255)
                                                                            tCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTgg             <  1:5135556/93‑1 (MQ=255)
                                                                            tCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTgg             <  1:3415692/93‑1 (MQ=255)
                                                                             cAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGt            >  2:221110/1‑93 (MQ=255)
                                                                              aGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTa           <  1:3174036/93‑1 (MQ=255)
                                                                              aGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTa           <  2:6188688/93‑1 (MQ=255)
                                                                               ggAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTAt          <  1:3330928/93‑1 (MQ=255)
                                                                                gAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATg         >  1:5504943/1‑93 (MQ=255)
                                                                                gAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATg         <  2:5449991/93‑1 (MQ=255)
                                                                                gAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATg         <  1:4095550/93‑1 (MQ=255)
                                                                                   gagaAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGt      >  1:3059644/1‑93 (MQ=255)
                                                                                   gagaAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGt      >  2:3834222/1‑93 (MQ=255)
                                                                                   gagaAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGt      <  1:1704183/93‑1 (MQ=255)
                                                                                   gagaAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGt      >  2:2441531/1‑93 (MQ=255)
                                                                                     gaAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTgg    <  2:4241219/93‑1 (MQ=255)
                                                                                     gaAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTgg    >  2:5008712/1‑93 (MQ=255)
                                                                                          |                                                                                         
GGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTC  >  NC_000913/4184984‑4185163

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