Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,262,255 G→C 16.7% S196T (AGC→ACC)  dusA → tRNA‑dihydrouridine synthase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,262,2550GC16.7% 282.0 / 49.3 108S196T (AGC→ACC) dusAtRNA‑dihydrouridine synthase A
Reads supporting (aligned to +/- strand):  ref base G (45/45);  new base C (12/6);  total (57/51)
Fisher's exact test for biased strand distribution p-value = 3.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.89e-01

Only 100 of 109 total aligned reads displayed.

GCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTA  >  NC_000913/4262156‑4262353
                                                                                                   |                                                                                                  
gCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAg                                                                                                    >  1:197423/1‑100 (MQ=255)
 cTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGc                                                                                                   >  1:2487981/1‑100 (MQ=255)
  tATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGcc                                                                                                  >  1:79265/1‑100 (MQ=255)
  tATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGcc                                                                                                  >  1:1249430/1‑100 (MQ=255)
   aTGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGccc                                                                                                 <  1:1684772/100‑1 (MQ=255)
    tGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCg                                                                                                <  1:1450409/100‑1 (MQ=255)
      aaTTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCgaa                                                                                              >  2:572347/1‑100 (MQ=255)
      aaTTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCgaa                                                                                              <  1:122791/100‑1 (MQ=255)
       aTTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCgaaa                                                                                             <  1:2521654/100‑1 (MQ=255)
        tttCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaag                                                                                            <  2:1907132/100‑1 (MQ=255)
        tttCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaag                                                                                            <  2:2043265/100‑1 (MQ=255)
        tttCTCTGCAATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGCGAGATGTTCATCATCCACGCACGAAAAGCCTGGCTTTCGGGGTTAAGCCCgaaag                                                                                            <  2:2301559/100‑1 (MQ=255)
          tctcTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaagcc                                                                                          <  1:1719278/100‑3 (MQ=255)
           ctctGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCgaaagaaa                                                                                         >  2:353713/1‑100 (MQ=255)
            tctGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaagccag                                                                                        <  1:321743/100‑5 (MQ=39)
            tctGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGGGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaagccag                                                                                        <  2:1028299/100‑5 (MQ=39)
             ctGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaagccagg                                                                                       <  1:1086156/100‑6 (MQ=39)
              tGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAAcc                                                                                      <  1:1537549/100‑1 (MQ=255)
               gCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCg                                                                                     >  2:1143697/1‑100 (MQ=255)
                  aTTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGa                                                                                  <  2:2235600/100‑1 (MQ=255)
                   tttCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGaa                                                                                 <  1:143280/100‑1 (MQ=255)
                   tttCATCAATACCGTTTCCGACAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGaa                                                                                 <  2:1698461/100‑1 (MQ=255)
                      catcaATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATc                                                                              >  1:569759/1‑100 (MQ=255)
                      catcaATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATc                                                                              >  2:70960/1‑100 (MQ=255)
                        tcaATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATccc                                                                            >  1:1298149/1‑100 (MQ=255)
                            caCCGTTTCCGGCAAAGGCGAGTGTGAGATGTTTATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGTCCGAAAGAAAACCGTGAAATCccgccg                                                                        <  1:941667/99‑1 (MQ=255)
                             aCCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCcgccgc                                                                       >  1:671255/1‑100 (MQ=255)
                              ccGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCt                                                                      >  2:1325947/1‑100 (MQ=255)
                               cGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTc                                                                     <  2:79264/100‑1 (MQ=255)
                                gTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCg                                                                    >  2:501018/1‑100 (MQ=255)
                                  ttCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGa                                                                   >  2:217217/1‑99 (MQ=255)
                                  ttCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGAt                                                                  <  2:2406650/100‑1 (MQ=255)
                                  ttCCAGCAAAGGCTAGTGTGAAATGTCCATCATCCACGCACGTAAAGCCGGGCTTTTGGCTTAAAGCCCGAAAGAAAACCGTGAAATCTCGCCGCTCg                                                                    <  2:1041248/98‑1 (MQ=255)
                                   tCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGAtt                                                                 >  2:1547942/1‑100 (MQ=255)
                                    cgggCAAGGGGGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTGCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGAt                                                                  <  2:1040197/96‑1 (MQ=255)
                                    ccGGCAAAGGCGAGTGTGAGATGTTCCTCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTa                                                                <  2:812576/100‑1 (MQ=255)
                                      ggCAAAGGCGAGTGTGAGATGTTCATCAGCCACGCACGTAGATCCTGGTTTCCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCGCGATTa                                                                <  2:1006940/98‑1 (MQ=255)
                                         aaaGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCcgc                                                           <  1:2255189/100‑1 (MQ=255)
                                          aaGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCcgcg                                                          <  2:320541/100‑1 (MQ=255)
                                           aGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCgt                                                         <  1:1051202/100‑1 (MQ=255)
                                            ggCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCgtg                                                        >  1:2442522/1‑100 (MQ=255)
                                            ggCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCgtg                                                        >  1:403140/1‑100 (MQ=255)
                                            ggCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCgtg                                                        <  1:1155471/100‑1 (MQ=255)
                                            ggCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCgtg                                                        <  2:2343266/100‑1 (MQ=255)
                                            ggCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCgtg                                                        <  2:1320847/100‑1 (MQ=255)
                                              cGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCgtgtg                                                      <  2:2156015/100‑1 (MQ=255)
                                              cGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCgtgtg                                                      <  2:1003130/100‑1 (MQ=255)
                                              cGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCgtgtg                                                      <  2:1377858/100‑1 (MQ=255)
                                                   gtgAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCa                                                 >  2:1712035/1‑100 (MQ=255)
                                                     gagaTGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAAc                                               >  1:553171/1‑100 (MQ=255)
                                                      agaTGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAAc                                               <  2:1283528/99‑1 (MQ=255)
                                                      agaTGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACt                                              <  1:1962183/100‑1 (MQ=255)
                                                          gTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAg                                          >  1:1691457/1‑100 (MQ=255)
                                                            tcatcaTCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCg                                        >  1:35511/1‑100 (MQ=255)
                                                            tcatcaTCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCg                                        <  1:1953449/100‑1 (MQ=255)
                                                             catcatCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGt                                       <  2:1604088/100‑1 (MQ=255)
                                                               tcatcCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCg                                        >  2:2114264/1‑97 (MQ=255)
                                                                catcCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGAc                                    <  2:2227654/100‑1 (MQ=255)
                                                                     acgcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTcc                               >  2:1453342/1‑100 (MQ=255)
                                                                     acgcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTcc                               >  1:1326000/1‑100 (MQ=255)
                                                                      cgcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCg                              <  1:1830044/100‑1 (MQ=255)
                                                                        cacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCa                            >  2:1661716/1‑100 (MQ=255)
                                                                        cacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCa                            >  1:164468/1‑100 (MQ=255)
                                                                            tAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTg                        <  2:1483788/100‑1 (MQ=255)
                                                                             aaaGCCTGGTTTTCGGGGTTAAGCCCTAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTg                        <  2:1647055/99‑1 (MQ=255)
                                                                               aGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACa                     <  1:1044647/100‑1 (MQ=255)
                                                                                ggttttcTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACaa                    >  1:1377724/7‑100 (MQ=38)
                                                                                gCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACaa                    >  1:2247230/1‑100 (MQ=255)
                                                                                 gttttcTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAAt                   >  1:1757233/6‑100 (MQ=39)
                                                                                 gttttcTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAAt                   >  1:1171012/6‑100 (MQ=39)
                                                                                 ccTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAACTCCAGCCGCTCGATTATCCGCGTGTGTATCAACTGGAGCGTGACTTTCCGCATCTGACAAt                   <  2:2394996/100‑1 (MQ=255)
                                                                                   tttcTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGt                 >  2:1675595/4‑100 (MQ=255)
                                                                                   tttcTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATAAACTAACGCGTGACTTTCCGCATc                          >  2:1245504/4‑91 (MQ=37)
                                                                                    ttcTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTc                >  2:183252/3‑100 (MQ=255)
                                                                                    ttcTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTc                >  2:763680/3‑100 (MQ=255)
                                                                                    ggCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTc                <  1:1096900/100‑1 (MQ=255)
                                                                                    ggCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTc                >  1:2447542/1‑100 (MQ=255)
                                                                                     gCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCg               >  1:208222/1‑100 (MQ=255)
                                                                                      cTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGa              <  1:2030962/100‑1 (MQ=255)
                                                                                      cTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGa              >  1:1553347/1‑100 (MQ=255)
                                                                                      cTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGa              >  1:235337/1‑100 (MQ=255)
                                                                                      cTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGa              >  1:421973/1‑100 (MQ=255)
                                                                                      cTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGa              >  1:781121/1‑100 (MQ=255)
                                                                                      cTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGa              >  1:1449316/1‑100 (MQ=255)
                                                                                       tttCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGAt             >  2:328685/1‑100 (MQ=255)
                                                                                       tttCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGAt             >  1:597619/1‑100 (MQ=255)
                                                                                       tttCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGAt             >  1:2338284/1‑100 (MQ=255)
                                                                                         tCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTa           >  2:1246352/1‑100 (MQ=255)
                                                                                         tCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTa           >  2:438808/1‑100 (MQ=255)
                                                                                          cGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTaa          >  1:1060792/1‑100 (MQ=255)
                                                                                           ggggTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAAATGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAAc         >  2:16562/1‑100 (MQ=255)
                                                                                             ggTTAAGCCCGAAAGAAATCCGTGAAATCCCTCCGCTCTATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACg        <  2:1700409/99‑1 (MQ=255)
                                                                                              gTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACggt      >  2:1360733/1‑100 (MQ=255)
                                                                                              gTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACggt      >  1:1290876/1‑100 (MQ=255)
                                                                                               ttAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACggtg     >  2:1674288/1‑100 (MQ=255)
                                                                                               ttAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACggtg     <  2:403136/100‑1 (MQ=255)
                                                                                                tAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACggtgg    <  1:784451/100‑1 (MQ=255)
                                                                                                tAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACggtgg    >  1:809383/1‑100 (MQ=255)
                                                                                                  aGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTa  <  1:2305118/100‑1 (MQ=255)
                                                                                                  aGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTa  >  2:2014253/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
GCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTA  >  NC_000913/4262156‑4262353

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