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.5% 1280.1 / NA 401Q716R (CAG→CGG) rpoCRNA polymerase beta prime subunit
Reads supporting (aligned to +/- strand):  ref base A (1/0);  major base G (184/214);  minor base T (0/2);  total (185/216)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 417 total aligned reads displayed.

GGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTC  >  NC_000913/4184984‑4185163
                                                                                          |                                                                                         
 ggCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGt                                                                                        >  1:4830583/1‑93 (MQ=255)
    tGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTc                                                                                     >  1:1895672/1‑93 (MQ=255)
    tGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTc                                                                                     >  1:640845/1‑93 (MQ=255)
     gcggcgAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTcc                                                                                    <  2:8720615/93‑1 (MQ=255)
     gcggcgAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTcc                                                                                    >  2:9410413/1‑93 (MQ=255)
      cggcgAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCt                                                                                   <  1:2090259/93‑1 (MQ=255)
      cggcgAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCt                                                                                   <  2:8356143/93‑1 (MQ=255)
        gcgAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTc                                                                                 >  2:5427092/1‑93 (MQ=255)
        gcgAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGGCAGGAAGAGAAGCGGGTTTCCTTc                                                                                 <  2:7051527/93‑1 (MQ=255)
          gAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCaa                                                                               >  1:8589259/1‑93 (MQ=255)
           aaCGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAAc                                                                              >  1:10020099/1‑93 (MQ=255)
           aaCGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAAc                                                                              <  2:4702547/93‑1 (MQ=255)
           aaCGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAAc                                                                              <  2:8246545/93‑1 (MQ=255)
           aaCGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAAc                                                                              <  2:5035685/93‑1 (MQ=255)
            aCGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACa                                                                             <  1:2418829/93‑1 (MQ=255)
            aCGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACa                                                                             <  2:3839427/93‑1 (MQ=255)
              gATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGc                                                                           <  2:7589528/93‑1 (MQ=255)
              gATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGc                                                                           >  2:6685912/1‑93 (MQ=255)
               aTCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCa                                                                          >  2:3972093/1‑93 (MQ=255)
                tCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCAt                                                                         <  2:4848863/93‑1 (MQ=255)
                   tgtATCCAAAGCGATGATGGATAACCTGCAATCTGAAACCGTGATTAGCCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTa                                                                      <  1:5565706/93‑1 (MQ=255)
                   tgtATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTa                                                                      >  1:4135648/1‑93 (MQ=255)
                   tgtATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTa                                                                      >  2:8197763/1‑93 (MQ=255)
                   tgtATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTa                                                                      <  2:20824/93‑1 (MQ=255)
                   tgtATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTa                                                                      >  2:1740529/1‑93 (MQ=255)
                   tgtATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTa                                                                      <  2:8412495/93‑1 (MQ=255)
                    gtATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTAc                                                                     <  1:110180/93‑1 (MQ=255)
                    gtATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTAc                                                                     >  1:4417213/1‑93 (MQ=255)
                    gtATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTAc                                                                     <  1:7362877/93‑1 (MQ=255)
                     tATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACa                                                                    >  1:6149803/1‑93 (MQ=255)
                     tATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACa                                                                    >  1:6875653/1‑93 (MQ=255)
                     tATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACa                                                                    <  1:4013355/93‑1 (MQ=255)
                        ccAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACatga                                                                 <  2:7295192/93‑1 (MQ=255)
                        ccAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACatga                                                                 <  2:1782773/93‑1 (MQ=255)
                        ccAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACatga                                                                 >  2:6143414/1‑93 (MQ=255)
                        ccAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACatga                                                                 <  2:6262453/93‑1 (MQ=255)
                         cAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACatgat                                                                <  1:9820074/93‑1 (MQ=255)
                           aaGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATgg                                                              <  1:5296641/93‑1 (MQ=255)
                           aaGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATgg                                                              >  1:8061164/1‑93 (MQ=255)
                              cGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCg                                                           <  1:9408398/93‑1 (MQ=255)
                                  gatgGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTc                                                       >  1:3363580/1‑93 (MQ=255)
                                  gatgGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGCACGACTc                                                       >  1:1452333/1‑93 (MQ=255)
                                    tgGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCg                                                     >  2:1881850/1‑93 (MQ=255)
                                    tgGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCg                                                     >  2:6388072/1‑93 (MQ=255)
                                    tgGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCg                                                     >  2:4425106/1‑93 (MQ=255)
                                    tgGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGGTTCCTTCAACAGCATCTACATGATGGCCGACTCCg                                                     >  2:92756/1‑93 (MQ=255)
                                     ggATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCgg                                                    >  1:2905230/1‑93 (MQ=255)
                                      gATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCggg                                                   >  2:7988779/1‑92 (MQ=255)
                                       aTAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTg                                                  >  1:2205003/1‑93 (MQ=255)
                                           ccTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGcgcg                                              >  1:9628638/1‑93 (MQ=255)
                                           ccTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGACAAGCTGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGcgcg                                              <  1:3224009/93‑1 (MQ=255)
                                            cTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGt                                             <  1:7182617/93‑1 (MQ=255)
                                            cTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGt                                             <  2:4112111/93‑1 (MQ=255)
                                            cTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGt                                             >  1:5662553/1‑93 (MQ=255)
                                            cTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGt                                             >  2:5239648/1‑93 (MQ=255)
                                              gCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTgg                                           <  1:6626528/93‑1 (MQ=255)
                                                  aCTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATTTACATGATGGCCGACTTCGGTGCGCGTGGTTCt                                       >  1:8983686/1‑93 (MQ=255)
                                                  aCTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCt                                       <  1:4629208/93‑1 (MQ=255)
                                                  aCTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCt                                       <  1:2767058/93‑1 (MQ=255)
                                                   ctaaaaCCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTg                                      <  2:4305608/90‑1 (MQ=255)
                                                     gAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCg                                    >  2:2843464/1‑93 (MQ=255)
                                                     gAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCg                                    >  2:9537630/1‑93 (MQ=255)
                                                        aCCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGca                                 <  1:1028881/93‑1 (MQ=255)
                                                            tGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAgc                             >  2:3604999/1‑92 (MQ=255)
                                                            tGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGa                             <  2:1631853/93‑1 (MQ=255)
                                                             gATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGAt                            <  1:9379703/93‑1 (MQ=255)
                                                             gATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGAt                            <  2:5661236/93‑1 (MQ=255)
                                                               ttAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTc                          >  1:8776918/1‑93 (MQ=255)
                                                               atAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTc                          <  1:3127586/92‑1 (MQ=255)
                                                                tAACCGTGATGGTCAGGAAGAGGAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTct                         >  2:9827155/1‑92 (MQ=255)
                                                                tAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCg                         >  1:4983742/1‑93 (MQ=255)
                                                                 aaCCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGt                        <  2:2761957/93‑1 (MQ=255)
                                                                   ccGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCa                      <  2:2653466/93‑1 (MQ=255)
                                                                   ccGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCa                      <  1:3700097/93‑1 (MQ=255)
                                                                   ccGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCCTGGTTCTGCGGACCAGAGTCGTCa                      >  2:3410854/1‑93 (MQ=255)
                                                                    cGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAg                     >  2:7841491/1‑93 (MQ=255)
                                                                     gTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGc                    >  2:10040622/1‑93 (MQ=255)
                                                                     gTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGc                    <  2:6828686/93‑1 (MQ=255)
                                                                      tGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCt                   <  1:1401060/93‑1 (MQ=255)
                                                                      tGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCt                   <  1:2661223/93‑1 (MQ=255)
                                                                        atggTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTg                 >  1:5104503/3‑93 (MQ=255)
                                                                         tggTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGc                <  1:3561556/92‑1 (MQ=255)
                                                                          ggTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCt               >  2:7629617/1‑93 (MQ=255)
                                                                          ggTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCt               <  2:1122360/93‑1 (MQ=255)
                                                                          ggTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCCACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCt               <  2:8847544/93‑1 (MQ=255)
                                                                           gTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTg              >  2:2203357/1‑93 (MQ=255)
                                                                            tCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTgg             <  1:8245433/93‑1 (MQ=255)
                                                                            tCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTgg             >  1:7671168/1‑93 (MQ=255)
                                                                            gcAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTgg             <  1:7830906/92‑1 (MQ=255)
                                                                              aGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTa           <  1:3819299/93‑1 (MQ=255)
                                                                              aGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTa           <  2:3150715/93‑1 (MQ=255)
                                                                               ggAAGAGAAGCGGGGTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTAt          >  1:9170248/1‑93 (MQ=255)
                                                                                gAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATg         <  1:6771588/93‑1 (MQ=255)
                                                                                  agagAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCg       >  2:1337767/1‑93 (MQ=255)
                                                                                  agagAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCg       >  1:124564/1‑93 (MQ=255)
                                                                                    agaAGCGTGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTg     <  1:8496318/93‑1 (MQ=255)
                                                                                    agaAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTg     <  2:5570971/93‑1 (MQ=255)
                                                                                     gaAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTgg    <  1:9103828/93‑1 (MQ=255)
                                                                                       aGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTc  <  2:738323/93‑1 (MQ=255)
                                                                                       aGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTc  <  2:9088613/93‑1 (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: