Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,261,358 T→A 16.9% intergenic (+52/‑311) yjbM → / → dusA uncharacterized protein/tRNA‑dihydrouridine synthase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,261,3580TA16.9% 478.5 / 76.6 187intergenic (+52/‑311)yjbM/dusAuncharacterized protein/tRNA‑dihydrouridine synthase A
Reads supporting (aligned to +/- strand):  ref base T (80/75);  new base A (15/17);  total (95/92)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 202 total aligned reads displayed.

ACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTT  >  NC_000913/4261259‑4261457
                                                                                                   |                                                                                                   
acTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAatat                                                                                                     <  1:2005812/100‑1 (MQ=255)
 cTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGGCATAATATTTGAAGTTCTTCGTTAGTATGGACGGCATTATATAAAAGGACCCAAtatt                                                                                                    >  2:2193258/1‑100 (MQ=255)
   tCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtattta                                                                                                  <  2:3491120/100‑1 (MQ=255)
   tCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtattta                                                                                                  <  2:3234872/100‑1 (MQ=255)
   tCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtattta                                                                                                  >  2:2091763/1‑100 (MQ=255)
   tCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtatt                                                                                                    >  1:1787156/1‑98 (MQ=255)
   tCTAAATTAATTTATGAGTTGTTCGGGGTTAAATCGGAGAAAGAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtattta                                                                                                  <  1:3105774/100‑1 (MQ=255)
     tAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtatttatt                                                                                                >  1:3935870/1‑100 (MQ=255)
      aaaTTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTg                                                                                               >  2:893086/1‑100 (MQ=255)
      aaaTTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTg                                                                                               <  2:3414302/100‑1 (MQ=255)
        attaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTggg                                                                                             <  1:1132398/100‑1 (MQ=255)
        attaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTggg                                                                                             >  2:4058517/1‑100 (MQ=255)
        attaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTggg                                                                                             <  1:3863709/100‑1 (MQ=255)
           aattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTc                                                                                          >  1:3664246/1‑100 (MQ=255)
           aattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGtcc                                                                                          <  2:1572840/100‑3 (MQ=39)
            attTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCt                                                                                         <  2:3951839/100‑1 (MQ=255)
            attTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCt                                                                                         >  1:2738609/1‑100 (MQ=255)
            attTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATAAATATTGGGTCCt                                                                                         <  2:1641889/100‑1 (MQ=255)
            attTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCt                                                                                         <  1:2647109/100‑1 (MQ=38)
            attTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTACGGTAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCt                                                                                         <  1:2078529/100‑1 (MQ=38)
             tttATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCtt                                                                                        >  1:3260427/1‑100 (MQ=255)
             tttATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCtt                                                                                        >  1:944664/1‑100 (MQ=255)
              ttATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCttt                                                                                       >  1:3603674/1‑100 (MQ=255)
               tATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATAAATATTGGGTCCtttt                                                                                      <  2:3088893/100‑1 (MQ=39)
                aTGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCttttt                                                                                     <  1:2778885/100‑1 (MQ=255)
                aTGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCTATAAATATTGGGTCCttttn                                                                                     <  1:2690265/100‑2 (MQ=38)
                aTGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCtttta                                                                                     <  1:1509324/100‑2 (MQ=38)
                 tGATTTGTTCGTGTTTAACTCGAATACATAATATATGAAGTTCTTCGATATTATTGAAGTCATTATATAACAGAACCCAATAAATAATGGGTCCttttat                                                                                    <  1:725790/100‑3 (MQ=16)
                 tGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCttttat                                                                                    <  2:2938480/100‑3 (MQ=37)
                  gATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtc                                                                                   <  2:1723858/100‑1 (MQ=255)
                  gATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtc                                                                                   <  2:3510659/100‑1 (MQ=255)
                  gATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtc                                                                                   <  2:77731/100‑1 (MQ=255)
                   aTTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtct                                                                                  >  1:1317931/1‑100 (MQ=255)
                   aTTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCttttatat                                                                                  <  2:4095553/100‑5 (MQ=37)
                      tGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTa                                                                               <  2:475719/100‑1 (MQ=255)
                       gTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTAt                                                                              >  2:1596550/1‑100 (MQ=255)
                        ttCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATc                                                                             <  2:535111/100‑1 (MQ=255)
                        gtCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATc                                                                             >  1:4086360/2‑100 (MQ=255)
                          cGGGGTTAAATCGGCGAAATAATCTGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAATACCCAATATTTATTGGATTCTTTTTTCTCTATCaa                                                                           >  1:813770/1‑100 (MQ=255)
                          cGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCaa                                                                           >  2:4248738/1‑100 (MQ=255)
                           ggggTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAAt                                                                          <  2:3802672/100‑1 (MQ=255)
                           ggggTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAAt                                                                          >  2:2125806/1‑100 (MQ=255)
                             ggTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATAc                                                                        >  2:1193677/1‑100 (MQ=255)
                                 aaaTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTAtt                                                                    >  1:2459134/1‑100 (MQ=255)
                                 aaaTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCAAATATTTATTGGGTTCTTTTTTCTCTATCAATACTAtt                                                                    >  1:349651/1‑100 (MQ=255)
                                   aTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAg                                                                  >  1:4311324/1‑100 (MQ=255)
                                    tCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGc                                                                 <  2:3888221/100‑1 (MQ=255)
                                    tCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGc                                                                 <  2:1173808/100‑1 (MQ=255)
                                     cGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCa                                                                >  2:231194/1‑100 (MQ=255)
                                     cGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTCTCTCTATCAATACTATTAGCa                                                                >  2:1772371/1‑100 (MQ=255)
                                      ggAGAAATAATATGAGGAGTTCCTCGATAGTATGGAAGGCATTATATAAAAGGACCCGATATTTATTGTGTTCTTTTTTCTCTATCAATACTATTAGCAg                                                               >  1:1435768/1‑100 (MQ=255)
                                       gagaAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAgg                                                              <  2:3050243/100‑1 (MQ=255)
                                        agaAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAggg                                                             <  2:3458587/100‑1 (MQ=255)
                                         gaAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGga                                                            >  1:1332235/1‑100 (MQ=255)
                                           aataatATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGgaga                                                          >  1:2316846/1‑100 (MQ=255)
                                               atatGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGatata                                                      <  1:106206/100‑1 (MQ=255)
                                                  tgtgAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCa                                                   >  2:2520246/1‑100 (MQ=255)
                                                    tgAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCAcc                                                 >  1:4370844/1‑100 (MQ=255)
                                                     gAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCa                                                <  2:2361356/100‑1 (MQ=255)
                                                     gAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCa                                                >  2:2424075/1‑100 (MQ=255)
                                                      aatttcttcGATATTATGGAATGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCa                                                <  1:3128972/96‑1 (MQ=255)
                                                        gTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCagag                                             >  2:732902/1‑100 (MQ=255)
                                                        gTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTATTTATCAATACTATTAGCAGGGAGATATATCACCag                                               >  1:1813559/1‑98 (MQ=255)
                                                         ttcttcGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGt                                            <  1:1948608/100‑1 (MQ=255)
                                                          tcttcGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGt                                            >  1:1703109/1‑99 (MQ=255)
                                                           cttcGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGttt                                          >  1:3307353/1‑100 (MQ=255)
                                                             tcGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTaa                                        >  2:1318489/1‑100 (MQ=255)
                                                               gATAGTCTGGTAGGCATTATATAAAAGCACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAACAGGCAGATCTTTCACCAGAGTTTACtg                                      <  1:3441822/100‑1 (MQ=255)
                                                                aTAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgt                                     >  2:3110680/1‑100 (MQ=255)
                                                                    tATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGa                                 <  2:2003244/100‑1 (MQ=255)
                                                                       ggAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAttt                              >  1:564658/1‑100 (MQ=255)
                                                                        gAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAtttt                             >  2:4194121/1‑100 (MQ=255)
                                                                         aaGGCATTATATAAAAGGACCCAATAATTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAttttt                            >  1:3645494/1‑100 (MQ=255)
                                                                             cATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttatt                        <  2:328565/100‑1 (MQ=255)
                                                                                tatataAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTttatttat                     <  2:1392915/100‑1 (MQ=255)
                                                                                gagaaaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATTTGTGATTTTttatttat                     >  1:2295482/6‑100 (MQ=25)
                                                                                gagaaaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTttatttat                     >  2:3982569/6‑100 (MQ=25)
                                                                                 agaaaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtc                    >  2:4404197/5‑100 (MQ=37)
                                                                                   ataAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTGTCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgt                  >  2:2678208/1‑100 (MQ=255)
                                                                                    taAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgtc                 >  2:2787946/1‑100 (MQ=255)
                                                                                     aaaaGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgtcg                <  2:3664271/100‑1 (MQ=255)
                                                                                     aaaaGAACCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgtcg                >  2:273945/1‑100 (MQ=39)
                                                                                      aaaGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGa               >  2:478977/1‑100 (MQ=255)
                                                                                      aaaGAACCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGa               >  2:299919/1‑100 (MQ=255)
                                                                                        agaaCCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTAATAGCAGGGCGATATATCACCAGAGGTTAATGTGTGATTTTTTATTTATCGTCGaa              >  1:1611265/4‑99 (MQ=255)
                                                                                         ggACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGa               >  1:292726/1‑97 (MQ=255)
                                                                                          gACACAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCt           <  1:401770/100‑1 (MQ=255)
                                                                                           aCCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTg          <  2:144531/100‑1 (MQ=255)
                                                                                           aCCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCAACAGAGTTTAATTTGTGATTTTTTATTTATCGTCGAAcc            >  1:4207605/1‑98 (MQ=255)
                                                                                            cccAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTgg         <  1:1866570/100‑1 (MQ=255)
                                                                                            cccAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTgg         >  2:3229084/1‑100 (MQ=255)
                                                                                            cccAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTgg         >  1:203217/1‑100 (MQ=255)
                                                                                             ccAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCt           >  1:3842154/1‑97 (MQ=255)
                                                                                             ccAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGa        <  2:3636637/100‑1 (MQ=255)
                                                                                             ccAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGa        >  1:466568/1‑100 (MQ=255)
                                                                                             ccAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGACTTTAATGTGTTATTTTTTATTTAtcat                  >  1:2867420/1‑88 (MQ=255)
                                                                                                atatTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTg     <  1:1172772/100‑1 (MQ=255)
                                                                                                 tatttattGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGt    >  2:2881542/1‑100 (MQ=255)
                                                                                                  atttattGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGtt   >  2:1453690/1‑100 (MQ=255)
                                                                                                   tttattGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGttt  >  2:3011440/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
ACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTT  >  NC_000913/4261259‑4261457

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: