Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,262,250 G→C 17.8% G194G (GGG→GGC dusA → tRNA‑dihydrouridine synthase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,262,2500GC17.8% 252.5 / 51.1 102G194G (GGG→GGCdusAtRNA‑dihydrouridine synthase A
Reads supporting (aligned to +/- strand):  ref base G (41/42);  new base C (12/6);  total (53/49)
Fisher's exact test for biased strand distribution p-value = 2.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.33e-01

Only 100 of 109 total aligned reads displayed.

GGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGT  >  NC_000913/4262151‑4262349
                                                                                                   |                                                                                                   
ggACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCgggg                                                                                                     <  2:363440/100‑1 (MQ=255)
ggACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCgggg                                                                                                     <  1:1582690/100‑1 (MQ=255)
ggACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCgggg                                                                                                     >  2:1953435/1‑100 (MQ=255)
 gACAGCTATGAATTTCTCTGCGATTGCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGt                                                                                                    <  2:1836141/100‑1 (MQ=255)
     gCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAg                                                                                                >  1:197423/1‑100 (MQ=255)
     gCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCAAGCACGTAAAGCCTGGCTTTCGGGGTTaa                                                                                                 >  1:1585182/1‑99 (MQ=255)
     gCTATGAATTGTTCAGCGATTTCATCAATACCTTTTCAGGAAAATACGAGTTTGAAATGTTCATCATCCACCCATGTAAATTGTGGCTTTCGGGGTTa                                                                                                  <  2:184988/98‑1 (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)
                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)
                       aTTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGa                                                                              <  2:2235600/100‑1 (MQ=255)
                        tttCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGaa                                                                             <  1:143280/100‑1 (MQ=255)
                        tttCATCAATACCGTTTCCGACAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGaa                                                                             <  2:1698461/100‑1 (MQ=255)
                           catcaATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATc                                                                          <  1:1763844/100‑1 (MQ=255)
                           catcaATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATc                                                                          >  2:70960/1‑100 (MQ=255)
                           catcaATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATc                                                                          >  1:569759/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)
                                   ccGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCt                                                                  >  1:2406664/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:1155471/100‑1 (MQ=255)
                                                 ggCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCgtg                                                    >  1:403140/1‑100 (MQ=255)
                                                 ggCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCgtg                                                    <  2:1320847/100‑1 (MQ=255)
                                                 ggCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCgtg                                                    <  2:2343266/100‑1 (MQ=255)
                                                 ggCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCgtg                                                    >  1:2442522/1‑100 (MQ=255)
                                                   cGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCgtgtg                                                  <  2:2156015/100‑1 (MQ=255)
                                                   cGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCgtgtg                                                  <  2:1377858/100‑1 (MQ=255)
                                                   cGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCgtgtg                                                  <  2:1003130/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                                      >  2:1044639/1‑100 (MQ=255)
                                                               gTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAg                                      >  1:1691457/1‑100 (MQ=255)
                                                                 tcatcaTCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCg                                    <  1:1953449/100‑1 (MQ=255)
                                                                 tcatcaTCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCg                                    >  1:35511/1‑100 (MQ=255)
                                                                 tcatcaTCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCg                                    >  2:1575722/1‑100 (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)
                                                                                   aaGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGAc                  <  1:630780/100‑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)
                                                                                           cTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGa          >  1:1449316/1‑100 (MQ=255)
                                                                                           cTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGa          >  1:421973/1‑100 (MQ=255)
                                                                                           cTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGa          >  1:781121/1‑100 (MQ=255)
                                                                                            tttCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGAt         >  2:328685/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  >  1:1290876/1‑100 (MQ=255)
                                                                                                   gTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACggt  >  2:1360733/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
GGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGT  >  NC_000913/4262151‑4262349

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: