Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,262,250 G→C 17.9% 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.9% 503.8 / 102.3 203G194G (GGG→GGCdusAtRNA‑dihydrouridine synthase A
Reads supporting (aligned to +/- strand):  ref base G (80/84);  new base C (19/17);  total (100/103)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.10e-01

Only 100 of 210 total aligned reads displayed.

GGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGG  >  NC_000913/4262151‑4262348
                                                                                                   |                                                                                                  
   cAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTa                                                                                                 >  1:798797/1‑100 (MQ=255)
     gCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAg                                                                                               <  2:1836328/100‑1 (MQ=255)
       tATGAATTTCTCTGCGATTTCATCAATACCGTTTCCTGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCGCGTAAAGCCTGGCTTTCGGGGTTAAGcc                                                                                             <  2:4087108/100‑1 (MQ=255)
        aTGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGccc                                                                                            >  2:488356/1‑100 (MQ=255)
        aTGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGccc                                                                                            >  1:3634710/1‑100 (MQ=255)
        aTGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGccc                                                                                            <  1:3539484/100‑1 (MQ=255)
           aaTTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCgaa                                                                                         >  1:4015480/1‑100 (MQ=255)
           aaTTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCg                                                                                           <  2:2348066/98‑1 (MQ=255)
           aaTTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaa                                                                                         <  1:2380643/100‑1 (MQ=255)
             tttCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCgaaag                                                                                       <  1:1020901/100‑1 (MQ=255)
             tttCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCgaaag                                                                                       <  1:2654995/100‑1 (MQ=255)
             tttCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCgaaa                                                                                        <  2:1287942/99‑1 (MQ=255)
             tttCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaag                                                                                       <  2:93197/100‑1 (MQ=255)
             tttCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaag                                                                                       <  1:3618708/100‑1 (MQ=255)
              ttCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaagc                                                                                      <  2:1616913/100‑2 (MQ=255)
              ttCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaagc                                                                                      <  2:753770/100‑2 (MQ=255)
               tctcTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaagcc                                                                                     <  1:3649692/100‑3 (MQ=255)
                ctctGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCTACGCATGTAAAGCCTGGCTTTCGGGGTTAAGCCCgaaagaaa                                                                                    <  2:705031/100‑1 (MQ=255)
                 tctGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaagccag                                                                                   <  2:1693351/100‑5 (MQ=39)
                   tcccgtttCCTCAATACCGGTTCCGGCAAAGGCGAGTGTGCGATGTGCCTCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaagcca                                                                                    <  2:106994/92‑4 (MQ=25)
                   tGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaagccaggc                                                                                 <  2:3218409/100‑7 (MQ=38)
                    gCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCg                                                                                <  1:3411436/100‑1 (MQ=255)
                    gCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCg                                                                                <  1:438578/100‑1 (MQ=255)
                    gCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaagccaggct                                                                                <  1:1310165/100‑8 (MQ=38)
                       atgtcatcaATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGa                                                                             <  1:3309266/97‑1 (MQ=255)
                       aTTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGa                                                                             <  2:2926134/100‑1 (MQ=255)
                         ttCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGaaa                                                                           >  2:2649735/1‑100 (MQ=255)
                          tcatcaATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAAt                                                                          >  1:4056233/1‑100 (MQ=255)
                          tcatcaATACCGTTTCCCGCAAAGTCGAGTGTGATATGTTCATCATTCACGCACGTAAAGCCTGGATTTCGGGGTTAATCTCTCAAGATAAACGTGaact                                                                          >  2:1355351/1‑98 (MQ=255)
                           catcaATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATc                                                                         >  1:3187622/1‑100 (MQ=255)
                           catcaATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATc                                                                         >  1:655750/1‑100 (MQ=255)
                           catcaATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATc                                                                         >  1:2071780/1‑100 (MQ=255)
                            atcaATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGAGGTTAAGCCCGAAAGAAAACCGTGAAATc                                                                         >  1:118520/1‑99 (MQ=255)
                             tcaATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATccc                                                                       >  2:206258/1‑100 (MQ=255)
                              caATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCccg                                                                      >  1:461695/1‑100 (MQ=255)
                               aaTACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCCTCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCccgc                                                                     <  1:2950569/100‑1 (MQ=255)
                               aaTACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCccgc                                                                     <  2:425723/100‑1 (MQ=255)
                               aaTACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCccgc                                                                     <  2:3789301/100‑1 (MQ=255)
                               aaTACCGTTGCCGGCAAAGGCGAGTGTGATATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCccgc                                                                     <  1:2048763/100‑1 (MQ=255)
                                     gTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCg                                                               >  2:310835/1‑100 (MQ=255)
                                     gTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCg                                                               >  1:3587948/1‑100 (MQ=255)
                                      tttCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCg                                                               <  2:3926149/99‑1 (MQ=255)
                                       ttCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGAt                                                             >  1:2459209/1‑100 (MQ=255)
                                         ccGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTa                                                           <  2:2963970/100‑1 (MQ=255)
                                          cGGCTAAGGCGAGTGTGAGATTTTCATCATCCACGCGCGTAAAGCCTGGCTTTTGGGATTAAGCCTGCAAGAAAACAGTGAAATCCCGCCGCTCg                                                               >  2:863140/1‑95 (MQ=255)
                                          cGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTAt                                                          >  1:2873145/1‑100 (MQ=255)
                                           ggCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATc                                                         >  2:2257908/1‑100 (MQ=255)
                                           ggCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAACCCCGCCGCTCGATTATc                                                         <  2:3008154/100‑1 (MQ=255)
                                            gCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATcc                                                        <  2:2156167/100‑1 (MQ=255)
                                                   cGAGTGTGGGATGTTCATAATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCgtgtg                                                 <  2:3017841/100‑1 (MQ=255)
                                                     aGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTa                                               >  1:3130122/1‑100 (MQ=255)
                                                       tgtgAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATc                                             <  1:2017746/100‑1 (MQ=255)
                                                         tgAGATGTTCACCATCCACGCACGGAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGAGCGTATCa                                            >  2:3622531/1‑99 (MQ=255)
                                                          gagaTGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGACATCCCGCCGCTCGATTATCCGCGTGTGTATCAAc                                          >  2:3331868/1‑100 (MQ=255)
                                                          gagaTGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAAc                                          <  2:517405/100‑1 (MQ=255)
                                                          gagaTGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAAc                                          <  1:2532546/100‑1 (MQ=255)
                                                            gaTGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTg                                        >  2:1457317/1‑100 (MQ=255)
                                                             aTGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACt                                         <  1:1750260/98‑1 (MQ=255)
                                                              tGTTCATCTTCCACGCACGTAAAGCCTGGCTTTCGGGTTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACt                                         <  2:1986955/97‑1 (MQ=255)
                                                               gTTCGTCATCGACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAGAACCGTGAACTCCCGCCGCTCGATTATCCGCGTGTGTATCAACTCAAg                                     >  2:1317175/1‑100 (MQ=255)
                                                               gTTCATCAGCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAg                                     <  1:3266636/100‑1 (MQ=255)
                                                                  catcatCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGt                                  >  1:1426152/1‑100 (MQ=255)
                                                                   atcatcCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTg                                 <  2:2080613/100‑1 (MQ=255)
                                                                    tcatcCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGa                                <  1:30746/100‑1 (MQ=255)
                                                                     catcCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGAc                               >  2:1899650/1‑100 (MQ=255)
                                                                      atcCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACt                              >  1:1829463/1‑100 (MQ=255)
                                                                          acgcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTcc                          <  1:3140717/100‑1 (MQ=255)
                                                                           cgcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCg                         <  1:2705730/100‑1 (MQ=255)
                                                                           cgcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCg                         <  1:3534980/100‑1 (MQ=255)
                                                                            gcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGc                        >  2:3740182/1‑100 (MQ=255)
                                                                             cacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCa                       >  2:1931737/1‑100 (MQ=255)
                                                                               cgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATc                     <  1:1980739/100‑1 (MQ=255)
                                                                               cgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATc                     >  1:3954082/1‑100 (MQ=255)
                                                                               cgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATc                     >  1:215063/1‑100 (MQ=255)
                                                                               cgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATc                     <  2:2572361/100‑1 (MQ=255)
                                                                               cgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATc                     >  1:617201/1‑100 (MQ=255)
                                                                                  cacggttttcTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTg                   >  1:2077737/10‑99 (MQ=38)
                                                                                    aGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACa                <  2:1200216/100‑1 (MQ=255)
                                                                                    aGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACa                <  2:355167/100‑1 (MQ=255)
                                                                                     gCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACa                >  2:3874652/1‑99 (MQ=255)
                                                                                      gttttcTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAAt              >  2:1819125/6‑100 (MQ=39)
                                                                                        tttcTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGt            >  2:3648094/4‑100 (MQ=255)
                                                                                         ttcTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGGCTATCAACCTAAGCGTGACTTTCCGCATCCGACAATg             >  1:917753/3‑98 (MQ=255)
                                                                                         ggCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTc           <  1:1751536/100‑1 (MQ=255)
                                                                                         ggCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCACATCTGACAATGTc           >  1:2059567/1‑100 (MQ=255)
                                                                                          tcTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCg          >  2:382227/2‑100 (MQ=255)
                                                                                          tcTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCg          >  1:1724711/2‑100 (MQ=255)
                                                                                          tcTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCg          >  2:3155369/2‑100 (MQ=255)
                                                                                          tcTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCg          >  2:2797762/2‑100 (MQ=255)
                                                                                           cTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGa         >  2:4120791/1‑100 (MQ=255)
                                                                                           cTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGa         >  1:687048/1‑100 (MQ=255)
                                                                                           cTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGa         >  1:2573643/1‑100 (MQ=255)
                                                                                           cTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGa         >  1:490431/1‑100 (MQ=255)
                                                                                            tttCGGGGTTAAGCCCGAAATAAAACCGTGAAACCCCGCCGCTCCATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCg          <  1:451703/98‑1 (MQ=255)
                                                                                            tttCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGAt        >  1:3775774/1‑100 (MQ=255)
                                                                                            tttCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGAt        >  2:2469013/1‑100 (MQ=255)
                                                                                             ttCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCTGCATCTGACAATGTCGAtt       >  1:1537463/1‑100 (MQ=255)
                                                                                              tCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTa      <  1:2251632/100‑1 (MQ=255)
                                                                                                 gggTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACg   >  2:579751/1‑100 (MQ=255)
                                                                                                  ggTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACgg  <  1:3153351/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
GGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGG  >  NC_000913/4262151‑4262348

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