Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 122,235 T→C R48R (CGT→CGC pdhR → transcriptional regulator for pyruvate dehydrogenase complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913122,2350TC99.4% 564.5 / ‑5.1 176R48R (CGT→CGCpdhRtranscriptional regulator for pyruvate dehydrogenase complex
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (76/99);  total (77/99)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

Only 100 of 183 total aligned reads displayed.

TGGAGTTTTTGATCCTCGAAGGCACTCTCCGCCCGGGCGAAAAACTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGTCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCTTGTTGCTTCGTCGCCAGGGTGGCGGCACTTTTGTCCAGAGCAGCCTATG  >  NC_000913/122147‑122327
                                                                                        |                                                                                            
tGGAGTTTTTGATCATCGAAGGCACTCTCCGCCCGGGCGAAAAACTGCCACCGGAACGCGAACTGGCAGAACAGTTTGACGTCTCCCGCCcct                                                                                          >  2:2165058/1‑93 (MQ=255)
tGGAGTTTTTGAGCCTCGAAGGCACTCTCCGCCCGGGCGAAAAACTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCcct                                                                                          <  1:368838/93‑1 (MQ=255)
 tgAGGTTTTTATCCTCCAAGGCACTCTCCGCCCCGGGGAAAAAATCCCACCGGAAAGCGAACTGGCAAAAAAAGTTGACGTCTCCcctccctc                                                                                         >  2:1968134/2‑93 (MQ=255)
  gAGTTTTTGATCCTCGAAGGCACTCTCCGCCCGGGCGAAAACCTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCcctcc                                                                                        <  1:3669425/93‑1 (MQ=255)
   aGTTTTTGATCCTCGAAGGCACTCTCCGCCCGGGCGAAAAACTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCcctcct                                                                                       <  1:80325/93‑1 (MQ=255)
   aGTTTTTGATCATCGAAGGCACTCTCCGCCCGGGCGAAAAACTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCcctcct                                                                                       <  2:2788595/93‑1 (MQ=255)
    gTTTTTGATCCTCGAAGGCACTCTCCGCCCGGGCGAAAAACTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCtt                                                                                      <  2:434707/93‑1 (MQ=255)
      ttttGATCCTCGAAGGCACTCTCCGCCCGGGCGAAAAACTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGc                                                                                    <  1:2635866/93‑1 (MQ=255)
      ttttGATCCTCGAAGGCACTCTCCGCCCGGGCGAAAAACTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGc                                                                                    <  2:761030/93‑1 (MQ=255)
        ttGATCCTCGAAGGCACTCTCCGCCCGGGCGAAAAACTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGt                                                                                  >  1:4733327/1‑93 (MQ=255)
        ttGATCCTCGAAGGCACTCTCCGCCCGGGCGAAAAACTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGt                                                                                  >  2:1954921/1‑93 (MQ=255)
         tGATCCTCGAAGGCACTCTCCGCCCGGGCGAAAAACTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTg                                                                                 >  2:2387008/1‑93 (MQ=255)
          gATCCTCGAAGGCACTCTCCGCCCGGGCGAAAAACTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGa                                                                                >  2:1836687/1‑93 (MQ=255)
           aTCCTCGAAGGCCACCTCCGCCCGGGCGAAAAACTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAg                                                                               <  2:1542611/93‑1 (MQ=255)
            tCCTCGAAGGCACTCTCCGCCCGGGCGAAAAACTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAgg                                                                              >  1:2310261/1‑93 (MQ=255)
            tCCTCGAAGGCACTCTCCGCCCGGGCGAAAAACTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAgg                                                                              >  1:659080/1‑93 (MQ=255)
            tCCTCGAAGGCACTCTCCGCCCGGGCGAAAAACTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGCGAgg                                                                              >  2:80777/1‑93 (MQ=255)
              cTCGAAGGCACTCTCCGCCCGGGCGAAAAACTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCg                                                                            >  1:1690134/1‑93 (MQ=255)
               tCGAAGGCACTCTCCGCCCGGGCGAAAAACTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGa                                                                           >  2:3690055/1‑93 (MQ=255)
               tCGAAGGCACTCTCCGCCCGGGCGAAAAACTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGa                                                                           >  2:3903345/1‑93 (MQ=255)
                cGAAGGCTCTCTCCGCCCGGGCGAAAAACTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGAt                                                                          <  1:4680124/93‑1 (MQ=255)
                 gAAGGCACTCTCCGCCCGGGCGAAAAACTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGAtt                                                                         >  2:3812064/1‑93 (MQ=255)
                 gAAGGCACTCTCCGCCCGGGCGAAAAACTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGAtt                                                                         >  1:1738361/1‑93 (MQ=255)
                   agccACTCTCCGCCCGGGCGAAAACCTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCa                                                                       <  2:3484298/90‑1 (MQ=255)
                    ggCACTCTCCGCCCGGGCGAAAAACTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCaa                                                                      <  2:2353644/93‑1 (MQ=255)
                        ctctCCGCCCGGGCGAAAAACTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGtc                                                                  >  1:158551/1‑93 (MQ=255)
                        ctctCCGCCCGGGCGAAAAACTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGtc                                                                  >  1:3355277/1‑93 (MQ=255)
                        ctctCCGCCCGGGCGAAAAACTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGtc                                                                  >  2:356162/1‑93 (MQ=255)
                        ctctCCGCCCGGGCGAAAAACTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGtc                                                                  >  1:2209613/1‑93 (MQ=255)
                        ctctCCGCCCGGGCGAAAAACTCCCACCGGAACGCGAACTGCCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGtc                                                                  <  1:3924156/93‑1 (MQ=255)
                         tctcCGCCCGGGCGAAAAACTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGtct                                                                 <  1:3081985/93‑1 (MQ=255)
                         tctcCGCCCGGGCGAAAAACTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGtct                                                                 <  1:1992201/93‑1 (MQ=255)
                         cctcCGCCCGGGCGAAAAACTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGtct                                                                 <  1:3468720/92‑1 (MQ=255)
                          ctcCGCCCGGGCGAAAAACTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGtctc                                                                >  1:1249871/1‑93 (MQ=255)
                          ctcCGCCCGGGCGAAAAACTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGtctc                                                                <  1:252637/93‑1 (MQ=255)
                          ctcCGCCCGGGCGAAAAACTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGtctc                                                                >  2:1557267/1‑93 (MQ=255)
                           tcCGCCCGGGCGAAAAACTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTcg                                                               <  2:3017467/93‑1 (MQ=255)
                           tcCGCCCGGGCGAAAAACTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTcg                                                               >  1:4696885/1‑93 (MQ=255)
                            ccGCCCGGGCGAAAAACTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTcga                                                              <  2:1152270/93‑1 (MQ=255)
                             cGCCCGGGCGAAAATCTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTcgaa                                                             <  2:2303583/93‑1 (MQ=255)
                             cGCCCGGGCGAAAAACTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTcgaa                                                             >  1:200722/1‑93 (MQ=255)
                             cGCCCGGGCGAAAAACTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTcgaa                                                             <  2:580084/93‑1 (MQ=255)
                              gCCCGGGCGAAAAACTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTcgaag                                                            >  2:1758548/1‑93 (MQ=255)
                               cccGGGCGAAAAACTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTcgaagc                                                           <  2:782720/93‑1 (MQ=255)
                               cccGGGCGAAAAACTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTcgaagc                                                           <  1:3186183/93‑1 (MQ=255)
                                ccGGGCGAAAAACTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTcgaagcg                                                          >  1:1068876/1‑93 (MQ=255)
                                   ggCGAAAAACTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTcgaagcgaag                                                       >  2:414950/1‑93 (MQ=255)
                                    gCGAAAAACTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAgg                                                      <  1:4448167/93‑1 (MQ=255)
                                     cGAAAACCTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAggg                                                     <  1:189351/93‑1 (MQ=255)
                                       aaaaaCTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCt                                                   <  1:4619582/93‑1 (MQ=255)
                                       aaaaaCTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCt                                                   >  2:2611183/1‑93 (MQ=255)
                                         aaaCTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCttg                                                 <  1:2796369/93‑1 (MQ=255)
                                          aaCTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCttgt                                                <  2:36257/93‑1 (MQ=255)
                                            cTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCttgttg                                              >  1:2258841/1‑93 (MQ=255)
                                             ccccACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCTTGTTGc                                             <  1:1618330/92‑1 (MQ=255)
                                              cccACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCTTGTTGCt                                            >  1:4158914/1‑93 (MQ=255)
                                              cccACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCTTGTTGCt                                            <  2:2706162/93‑1 (MQ=255)
                                               ccACGGGAACGCGAACGGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCTTGTTGCtt                                           <  2:2194100/93‑1 (MQ=255)
                                               ccACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCTTGTTGCtt                                           <  2:917609/93‑1 (MQ=255)
                                                cACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCTTGTTGCTtc                                          <  1:1382622/93‑1 (MQ=255)
                                                 aCCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCTTGTTGCTtcg                                         >  1:3248303/1‑93 (MQ=255)
                                                 aCCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCTTGTTGCTtcg                                         >  2:1487659/1‑93 (MQ=255)
                                                  cgggAACGCGAACTGGCAAAACAGTTTGATGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCTTGTTGCTtcgt                                        <  2:1916451/91‑1 (MQ=255)
                                                   cGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCTTGTTGCTtcgtc                                       >  2:3176586/1‑93 (MQ=255)
                                                    ggAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCTTGTTGCTtcgtcg                                      >  2:987116/1‑93 (MQ=255)
                                                       aCGCGAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCTTGTTGCTTCGTCGCCa                                   >  1:1055902/1‑93 (MQ=255)
                                                        cgcgAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCTTGTTGCTTCGTCGCCAg                                  <  1:2377447/93‑1 (MQ=255)
                                                        cgcgAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCTTGTTGCTTCGTCGCCAg                                  <  1:2974562/93‑1 (MQ=255)
                                                         gcgAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAGGGGCTTGTTGCTTCGTCGCCAgg                                 <  2:1046775/93‑1 (MQ=255)
                                                          cgAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCTTGTTGCTTCGTCGCCAggg                                <  2:795318/93‑1 (MQ=255)
                                                          cgAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAAGGCTTGTTGCTTCGTCGCCAggg                                >  1:3003026/1‑93 (MQ=255)
                                                           gAACTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCTTGTTGCTTCGTCGCCAGGGt                               <  2:1857747/93‑1 (MQ=255)
                                                            aaCTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCTTGTTGCTTCGTCGCCAGGGTg                              <  2:2141332/93‑1 (MQ=255)
                                                            aaCTGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCATGTTGCTTCGTCGCCACGGTg                              >  2:1241666/1‑93 (MQ=255)
                                                               tGGCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCTTGTTGCTTCGTCGCCAGGGTggcg                           >  2:3213445/1‑93 (MQ=255)
                                                                ggCAAAACAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCTTGTTGCTTCGTCGCCAGGGTggcgg                          <  1:3979293/93‑1 (MQ=255)
                                                                   aaaaCAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCTTGTTGCTTCGTCGCCAGGGTggcggccc                       >  2:3253932/1‑91 (MQ=255)
                                                                   aaaaCAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCTTGTTGCTTCGTCGCCAGGGTGGCGGCAc                       >  2:4581760/1‑93 (MQ=255)
                                                                   aaaaCAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCTTGTTGCTTCGTCGCCAGGGTGGCGGCAc                       >  1:565009/1‑93 (MQ=255)
                                                                    aaaCAGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCTTGTTGCTTCGTCGCCAGGGTGGCGGCACt                      >  1:2421931/1‑93 (MQ=255)
                                                                        aGTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCTTGTTGCTTCGTCGCCAGGGTGGCGGCACTTTTg                  <  2:566622/93‑1 (MQ=255)
                                                                         gTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCTTGTTGCTTCGTCGCCAGGGTGGCGGCACTTTTGt                 >  2:2872436/1‑93 (MQ=255)
                                                                         gTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCTTGTTGCTTCGTCGCCAGGGTGGCGGCACTTTTGt                 >  2:3673768/1‑93 (MQ=255)
                                                                         gTTTGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCTTGTTGCTTCGTCGCCAGGGGGGGGGCACTTTTGt                 >  1:903338/1‑93 (MQ=255)
                                                                           ttGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCTTGTTGCTTCGTCGCCAGGGTGGCGGCACTTTTGTcc               >  1:3890598/1‑93 (MQ=255)
                                                                           ttGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCTTGTTGCTTCGTCGCCAGGGTGGCGGCACTTTTGTcc               <  1:1480905/93‑1 (MQ=255)
                                                                           ttGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCTTGTTGCTTCGTCGCCAGGGTGGCGGCACTTTTGTcc               <  2:1527728/93‑1 (MQ=255)
                                                                           ttGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCTTGTTGCTTCGTCGCCAGGGTGGCGGCACTTTTGTcc               >  2:4009742/1‑93 (MQ=255)
                                                                            tGACGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCTTGTTGCTTCGTCGCCAGGGTGGCGGCACTTTTGTCCa              <  2:4728146/93‑1 (MQ=255)
                                                                              cggTCCCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCTTGTTGCTTCGTCGCCAGGGTGGCGGCACTTTTGTCCaga            <  2:3396034/91‑1 (MQ=255)
                                                                              aCGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCTTGTTGCTTCGTCGCCAGGGTGGCGGCACTTTTGTCCaga            <  2:1214123/93‑1 (MQ=255)
                                                                               cGTCTCCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCTTGTTGCTTCGTCGCCAGGGTGGCGGCACTTTTGTCCagag           <  1:250185/93‑1 (MQ=255)
                                                                                 tctcCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCTTGTTGCTTCGTCGCCAGGGTGGCGGCACTTTTGTCCAGagca         <  1:995960/93‑1 (MQ=255)
                                                                                 tctcCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCTTGTTGCTTCGGCGCCAGGGTGGCGGCCCTTTTGTCCCGcgca         >  2:1774606/1‑93 (MQ=255)
                                                                                  ctcCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCTTGTTGCTTCGTCGCCAGGGTGGCGGCACTTTTGTCCAGagcag        <  2:4661741/93‑1 (MQ=255)
                                                                                   tcCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCTTGTTGCTTCGTCGCCAGGGTGGCGGCACTTTTGTCCAGagcagc       >  1:1429790/1‑93 (MQ=255)
                                                                                   tcCCGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCTTGTTGCTTCGTCGCCAGGGTGGCGGCACTTTTGTCCAGagcagc       >  2:2738569/1‑93 (MQ=255)
                                                                                    cccGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCTTGTTGCTTCGTCGCCAGGGTGGCGGCACTTTTGTCCAGAGCAGcc      <  1:3034546/93‑1 (MQ=255)
                                                                                     ccGCCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCTTGTTGCTTCGTCGCCAGGGTGGCAGCACTTTTGTCCAGAGCAGCCt     <  1:2870734/93‑1 (MQ=255)
                                                                                        tCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCTTGTTGCTTCGTCGCCAGGGTGGCGGCACTTTTGTCCAGAGCAGCCTATg  >  1:4033000/1‑93 (MQ=255)
                                                                                        |                                                                                            
TGGAGTTTTTGATCCTCGAAGGCACTCTCCGCCCGGGCGAAAAACTCCCACCGGAACGCGAACTGGCAAAACAGTTTGACGTCTCCCGTCCCTCCTTGCGTGAGGCGATTCAACGTCTCGAAGCGAAGGGCTTGTTGCTTCGTCGCCAGGGTGGCGGCACTTTTGTCCAGAGCAGCCTATG  >  NC_000913/122147‑122327

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