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.3% 982.3 / NA 303Q716R (CAG→CGG) rpoCRNA polymerase beta prime subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (161/140);  minor base T (1/1);  total (162/141)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 312 total aligned reads displayed.

GGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTC  >  NC_000913/4184984‑4185163
                                                                                          |                                                                                         
gggCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCggg                                                                                         <  2:2869775/93‑1 (MQ=255)
gggCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCggg                                                                                         <  1:811393/93‑1 (MQ=255)
 ggCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGt                                                                                        <  1:3788171/93‑1 (MQ=255)
 ggCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGt                                                                                        >  2:2081278/1‑93 (MQ=255)
  gCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGtt                                                                                       <  1:1299616/93‑1 (MQ=255)
     gcggcgAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTcc                                                                                    >  1:5587941/1‑93 (MQ=255)
      cggcgAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCt                                                                                   <  2:2008211/93‑1 (MQ=255)
      cggcgAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCt                                                                                   <  1:2951439/93‑1 (MQ=255)
      cggcgAACGATCGTGTATCCAAAGAGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCt                                                                                   <  2:5352429/93‑1 (MQ=255)
        gcgAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTc                                                                                 <  1:4909182/93‑1 (MQ=255)
        gcgAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTc                                                                                 >  1:1762395/1‑93 (MQ=255)
        gcgAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTc                                                                                 >  2:3907430/1‑93 (MQ=255)
         cgAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCa                                                                                <  2:6130416/93‑1 (MQ=255)
          gAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCaa                                                                               >  1:1571331/1‑93 (MQ=255)
              gATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGc                                                                           >  2:6411738/1‑93 (MQ=255)
               aTCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCa                                                                          >  1:3335471/1‑93 (MQ=255)
               aTCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCa                                                                          >  2:1175992/1‑93 (MQ=255)
                 cGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATc                                                                        <  1:5903034/93‑1 (MQ=255)
                  gtgtATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCt                                                                       <  1:1995405/93‑1 (MQ=255)
                  gtgtATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCt                                                                       >  1:1093504/1‑93 (MQ=255)
                   tgtATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTa                                                                      >  1:6689133/1‑93 (MQ=255)
                   tgtATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTa                                                                      <  2:4434854/93‑1 (MQ=255)
                    gtATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTAc                                                                     >  2:1586084/1‑93 (MQ=255)
                    gtATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTAc                                                                     >  2:857518/1‑93 (MQ=255)
                     tATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTTAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACa                                                                    <  1:3171846/93‑1 (MQ=255)
                     tATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACa                                                                    <  2:4499219/93‑1 (MQ=255)
                       tCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACatg                                                                  >  2:743963/1‑93 (MQ=255)
                        ccAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACatga                                                                 <  2:4261794/93‑1 (MQ=255)
                         cAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACatgat                                                                >  1:1367213/1‑93 (MQ=255)
                             gCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGCTTCCTTCAACAGCATCTACATGATGGcc                                                            >  2:3716213/1‑93 (MQ=255)
                              cGATGATGGATACCCTGCAACCTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGAGTTTCCTTCAACAGCATCTACATGATGGCCg                                                           <  1:2408536/93‑1 (MQ=255)
                              cGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCg                                                           <  1:57035/93‑1 (MQ=255)
                               gatgatGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGa                                                          >  2:3409091/1‑93 (MQ=255)
                                 tgatgGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACt                                                        >  2:4766271/1‑93 (MQ=255)
                                 tgatgGATAACCTGCAAACTGAAACCGTGATTAACCGGGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCGACATGATGGCCGAcg                                                        >  1:6209982/1‑92 (MQ=255)
                                   atgGATAACCTGCAAACTGAATCCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGGTTGCTTCAACAGCATCTACATGATGGCCGACTcc                                                      >  2:3160394/1‑93 (MQ=255)
                                     ggATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCgg                                                    >  2:2414153/1‑93 (MQ=255)
                                     ggATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCgg                                                    >  1:2839225/1‑93 (MQ=255)
                                     ggATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCgg                                                    >  1:958702/1‑93 (MQ=255)
                                     ggATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCgg                                                    <  1:6478640/93‑1 (MQ=255)
                                       aTAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGCAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACCCAGGTg                                                  >  1:5938145/1‑93 (MQ=255)
                                       aTAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTg                                                  >  1:5032640/1‑93 (MQ=255)
                                         aaCCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTgcg                                                >  2:4975539/1‑93 (MQ=255)
                                         aaCCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGACGACTCCCGTgcg                                                >  2:464332/1‑93 (MQ=255)
                                          aCCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCTGGTTTCCTTCAACAGCATCTACATGATGGCAGACTCCGGTgcgc                                               >  1:2277360/1‑93 (MQ=255)
                                           ccTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGcgcg                                              >  2:104658/1‑93 (MQ=255)
                                            cTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGcgcgg                                             >  2:3813188/1‑92 (MQ=255)
                                            cTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGt                                             >  2:6034370/1‑93 (MQ=255)
                                            cTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGt                                             <  1:4774801/93‑1 (MQ=255)
                                               cAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGt                                          <  1:1729592/93‑1 (MQ=255)
                                                 aaCTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTc                                        <  2:5435553/93‑1 (MQ=255)
                                                 aaCTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCAGACTCCGGTGCGCGTGGTTc                                        >  2:2707618/1‑93 (MQ=255)
                                                   cTGAAGCCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTg                                      <  1:2716656/93‑1 (MQ=255)
                                                   cTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTg                                      >  2:2743672/1‑93 (MQ=255)
                                                    tGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGc                                     >  2:5628956/1‑93 (MQ=255)
                                                    tGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGc                                     >  1:4453052/1‑93 (MQ=255)
                                                    tGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGc                                     <  1:2557093/93‑1 (MQ=255)
                                                    tGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGc                                     <  2:4047433/93‑1 (MQ=255)
                                                    tGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGc                                     >  1:6584875/1‑93 (MQ=255)
                                                     gAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCg                                    >  2:5308496/1‑93 (MQ=255)
                                                     gAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCg                                    >  1:6274333/1‑93 (MQ=255)
                                                       aaCCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGc                                  >  2:6262256/1‑93 (MQ=255)
                                                       aaCCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGc                                  <  2:1371286/93‑1 (MQ=255)
                                                       aaCCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGc                                  >  1:2489695/1‑93 (MQ=255)
                                                        aCCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGca                                 >  2:2480194/1‑93 (MQ=255)
                                                        aCCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGca                                 <  2:1789823/93‑1 (MQ=255)
                                                           gTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAg                              <  1:1011696/93‑1 (MQ=255)
                                                            tGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGa                             <  1:3607816/93‑1 (MQ=255)
                                                            tGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGa                             >  2:3648621/1‑93 (MQ=255)
                                                             gATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGAt                            <  2:1057785/93‑1 (MQ=255)
                                                             gATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGAt                            <  2:1936752/93‑1 (MQ=255)
                                                             gATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGAt                            <  1:2105253/93‑1 (MQ=255)
                                                               ttAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGGACAGATTc                          >  1:3290117/1‑93 (MQ=255)
                                                               ttAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTc                          >  2:6465137/1‑93 (MQ=255)
                                                                tAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCg                         >  2:5188079/1‑93 (MQ=255)
                                                                  aCCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTc                       >  2:5463272/1‑93 (MQ=255)
                                                                  aCCGTGATGCTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTc                       <  1:5352119/93‑1 (MQ=255)
                                                                  aCCATGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAGCAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTc                       >  2:1507531/1‑93 (MQ=255)
                                                                    cGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAg                     >  2:3328404/1‑93 (MQ=255)
                                                                    cGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAg                     <  2:5740266/93‑1 (MQ=255)
                                                                    cGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAg                     >  2:1330410/1‑93 (MQ=255)
                                                                     gTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGc                    >  2:5871650/1‑93 (MQ=255)
                                                                     gTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGc                    >  1:1869057/1‑93 (MQ=255)
                                                                      tGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCt                   <  1:431538/93‑1 (MQ=255)
                                                                       gatggTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCtt                  >  2:1867419/4‑93 (MQ=255)
                                                                         tggTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGc                <  1:4597986/92‑1 (MQ=255)
                                                                         tggTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGc                <  2:5587941/92‑1 (MQ=255)
                                                                          ggTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCt               >  1:3067018/1‑93 (MQ=255)
                                                                            tCAGGCAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCACCTTGCTgg             >  1:5403014/1‑93 (MQ=255)
                                                                             cAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGt            <  2:2646291/93‑1 (MQ=255)
                                                                             cAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGt            >  2:4190222/1‑93 (MQ=255)
                                                                              aGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTa           >  1:1933239/1‑93 (MQ=255)
                                                                               ggAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTAt          >  1:1503296/1‑93 (MQ=255)
                                                                               ggAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTAt          >  2:2950554/1‑93 (MQ=255)
                                                                               ggAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTAt          <  1:4526064/93‑1 (MQ=255)
                                                                               ggAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTAt          <  1:148312/93‑1 (MQ=255)
                                                                               ggAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTAt          >  1:3232950/1‑93 (MQ=255)
                                                                                   gagaAGCGTGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGt      <  2:2793729/93‑1 (MQ=255)
                                                                                     gaAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTgg    >  2:6614183/1‑93 (MQ=255)
                                                                                       aGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTc  >  2:5110690/1‑93 (MQ=255)
                                                                                          |                                                                                         
GGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTC  >  NC_000913/4184984‑4185163

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