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.0% 642.1 / ‑3.5 197Q716R (CAG→CGG) rpoCRNA polymerase beta prime subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (97/98);  minor base T (2/0);  total (99/98)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.48e-01

Only 100 of 201 total aligned reads displayed.

GGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTC  >  NC_000913/4184984‑4185163
                                                                                          |                                                                                         
 ggCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGt                                                                                        <  1:3354895/93‑1 (MQ=255)
 ggCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGt                                                                                        <  1:2015149/93‑1 (MQ=255)
     gcggcgAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTcc                                                                                    >  2:2975874/1‑93 (MQ=255)
     gcggcgAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTcc                                                                                    >  2:2804082/1‑93 (MQ=255)
      cggcgAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCt                                                                                   >  2:4327044/1‑93 (MQ=255)
      cggcgAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCt                                                                                   >  1:3465573/1‑93 (MQ=255)
      cggcgAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCt                                                                                   <  1:3891306/93‑1 (MQ=255)
          gAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCaa                                                                               <  1:615990/93‑1 (MQ=255)
          gAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCaa                                                                               >  1:4011183/1‑93 (MQ=255)
          gAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCaa                                                                               >  1:2607354/1‑93 (MQ=255)
           aaCGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAAc                                                                              <  1:1695136/93‑1 (MQ=255)
           aaCGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAAc                                                                              <  1:4263393/93‑1 (MQ=255)
           aaCGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAAc                                                                              <  2:264058/93‑1 (MQ=255)
           aaCGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAAc                                                                              <  1:2717689/93‑1 (MQ=255)
            aCGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACa                                                                             <  2:713836/93‑1 (MQ=255)
               aTCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGTTTTCCTTCAACAGCa                                                                          >  1:4226821/1‑93 (MQ=255)
               aTCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCa                                                                          >  2:1186589/1‑93 (MQ=255)
               aTCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCa                                                                          <  1:3066652/93‑1 (MQ=255)
               aTCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCa                                                                          >  1:355467/1‑93 (MQ=255)
               aTCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCa                                                                          >  1:3848463/1‑93 (MQ=255)
               aTCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGCCAGGAAGAGAAGCGGGTTTCCTTCAACAGCa                                                                          >  1:3947495/1‑93 (MQ=255)
                 cGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATc                                                                        >  2:256243/1‑93 (MQ=255)
                  gtgtATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCt                                                                       >  1:2638273/1‑93 (MQ=255)
                  gtgtATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCt                                                                       <  2:4129910/93‑1 (MQ=255)
                     tATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACa                                                                    <  1:78900/93‑1 (MQ=255)
                       tCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACatg                                                                  >  2:1836445/1‑93 (MQ=255)
                        ccAAAGCGATGATGGATAACCTGCAAGCTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACatga                                                                 <  2:1241385/93‑1 (MQ=255)
                        ccAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACatga                                                                 <  1:1034020/93‑1 (MQ=255)
                        ccAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACatga                                                                 <  1:1250988/93‑1 (MQ=255)
                         cAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACatgat                                                                >  2:2014607/1‑93 (MQ=255)
                          aaaGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACatgatg                                                               >  2:1471596/1‑93 (MQ=255)
                           aaGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAATAGCATCTACATGATgg                                                              <  2:4508572/93‑1 (MQ=255)
                            aGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGc                                                             <  1:4683354/93‑1 (MQ=255)
                             gCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGcc                                                            <  2:1462239/93‑1 (MQ=255)
                             gCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGcc                                                            >  2:900188/1‑93 (MQ=255)
                             gCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGcc                                                            >  2:2929765/1‑93 (MQ=255)
                             gCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGcc                                                            <  2:3678070/93‑1 (MQ=255)
                              ggatgatGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCg                                                           <  2:4670711/92‑1 (MQ=255)
                               gatgatGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGa                                                          >  1:1232204/1‑93 (MQ=255)
                               gatgatGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGa                                                          <  2:4609177/93‑1 (MQ=255)
                                atgatgGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGAc                                                         >  2:3545587/1‑93 (MQ=255)
                                atgatgGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGAc                                                         <  1:1731813/93‑1 (MQ=255)
                                 tgatgGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACt                                                        >  1:392396/1‑93 (MQ=255)
                                 tgatgGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACt                                                        >  2:4744787/1‑93 (MQ=255)
                                 tgatgGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACt                                                        >  2:1104126/1‑93 (MQ=255)
                                  gatgGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTc                                                       <  2:2984329/93‑1 (MQ=255)
                                   atgGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTcc                                                      <  2:387457/93‑1 (MQ=255)
                                   atgGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTcc                                                      <  1:188519/93‑1 (MQ=255)
                                   atgGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTcc                                                      <  1:919232/93‑1 (MQ=255)
                                     ggATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAGCAGCATCTACATGATGGCCGACTCCgg                                                    <  1:319725/93‑1 (MQ=255)
                                     ggATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCgg                                                    <  1:4148199/93‑1 (MQ=255)
                                      gATACCCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGt                                                   <  2:1730599/93‑1 (MQ=255)
                                      gATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGt                                                   <  1:752801/93‑1 (MQ=255)
                                      gATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGt                                                   >  1:959490/1‑93 (MQ=255)
                                      gATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGt                                                   >  1:1195891/1‑93 (MQ=255)
                                      gATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGt                                                   >  2:3034982/1‑93 (MQ=255)
                                        tAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTgc                                                 <  1:692867/93‑1 (MQ=255)
                                        tAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTgc                                                 >  2:4073539/1‑93 (MQ=255)
                                        tAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTgc                                                 <  1:2468396/93‑1 (MQ=255)
                                        tAACCTGCAAACTGAAACCGAAATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTgc                                                 <  2:4333230/93‑1 (MQ=255)
                                         aaCCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTgcg                                                >  1:2978471/1‑93 (MQ=255)
                                          aCCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTgcgc                                               >  2:445847/1‑93 (MQ=255)
                                               cAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGt                                          <  1:540428/93‑1 (MQ=255)
                                                aaaCTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGtt                                         >  2:3937548/1‑93 (MQ=255)
                                                aaaCTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGtt                                         <  2:4085879/93‑1 (MQ=255)
                                                 aaCTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTc                                        <  1:1580840/93‑1 (MQ=255)
                                                 aaCTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTc                                        <  2:1626197/93‑1 (MQ=255)
                                                 aaCTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTc                                        <  2:1915479/93‑1 (MQ=255)
                                                  aCTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCt                                       <  1:2675014/93‑1 (MQ=255)
                                                    tGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGc                                     >  1:2956390/1‑93 (MQ=255)
                                                     gAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCg                                    <  2:2279509/93‑1 (MQ=255)
                                                       aaCCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGc                                  >  1:882386/1‑93 (MQ=255)
                                                       aaCCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGc                                  <  1:4397242/93‑1 (MQ=255)
                                                        aCCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGca                                 >  2:633360/1‑93 (MQ=255)
                                                        aCCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGca                                 >  1:3396232/1‑93 (MQ=255)
                                                           gTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAg                              <  1:165418/93‑1 (MQ=255)
                                                             gATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGAt                            >  2:1317862/1‑93 (MQ=255)
                                                             gATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGAt                            <  2:2162896/93‑1 (MQ=255)
                                                              aTTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGAtt                           >  2:2966635/1‑93 (MQ=255)
                                                              aTTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGAtt                           <  1:2564729/93‑1 (MQ=255)
                                                               ttAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTc                          >  1:1170799/1‑93 (MQ=255)
                                                                  aCCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTc                       >  1:2773201/1‑93 (MQ=255)
                                                                      tGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCt                   >  2:246733/1‑93 (MQ=255)
                                                                      tGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCt                   >  1:2342207/1‑93 (MQ=255)
                                                                         tggTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGc                <  1:3018387/92‑1 (MQ=255)
                                                                         tggTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGc                >  2:3169221/2‑93 (MQ=255)
                                                                         tggTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGc                <  2:504681/92‑1 (MQ=255)
                                                                           gTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTg              <  2:1384049/93‑1 (MQ=255)
                                                                            tCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTgg             >  2:3959018/1‑93 (MQ=255)
                                                                              aGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTa           <  2:427288/93‑1 (MQ=255)
                                                                              aGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTa           >  1:487270/1‑93 (MQ=255)
                                                                               ggAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTAt          >  1:681398/1‑93 (MQ=255)
                                                                               ggAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCAAAGATTCGTCAGCTTGCTgggag          >  1:3623256/1‑90 (MQ=255)
                                                                                gAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATg         <  2:1756320/93‑1 (MQ=255)
                                                                                gAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATg         <  1:4642268/93‑1 (MQ=255)
                                                                                 aaGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGc        <  1:1353023/93‑1 (MQ=255)
                                                                                      aaGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGt   >  1:3059806/1‑93 (MQ=255)
                                                                                       aGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTc  >  1:2119512/1‑93 (MQ=255)
                                                                                       aGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTc  >  2:3253328/1‑93 (MQ=255)
                                                                                       aGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTc  >  1:1005157/1‑93 (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: