Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,261,359 T→A 15.6% intergenic (+53/‑310) 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,3590TA15.6% 490.7 / 75.6 195intergenic (+53/‑310)yjbM/dusAuncharacterized protein/tRNA‑dihydrouridine synthase A
Reads supporting (aligned to +/- strand):  ref base T (82/81);  new base A (12/18);  total (94/101)
Fisher's exact test for biased strand distribution p-value = 3.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.76e-01

Only 100 of 211 total aligned reads displayed.

CTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTA  >  NC_000913/4261260‑4261458
                                                                                                    |                                                                                                   
cTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtatt                                                                                                     >  2:846834/1‑100 (MQ=255)
  tCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtattta                                                                                                   >  1:2917180/1‑100 (MQ=255)
  tCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtattta                                                                                                   >  2:3287850/1‑100 (MQ=255)
  tCTAAATAAATCTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGAATGGAAGGCATTATATAAAAGGACCCAAtattta                                                                                                   >  1:2242058/1‑100 (MQ=255)
   cTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtatttat                                                                                                  >  2:3107271/1‑100 (MQ=255)
    tAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtatttatt                                                                                                 <  2:42187/100‑1 (MQ=255)
     aaaTTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTg                                                                                                >  2:1475096/1‑100 (MQ=255)
      aattaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTgg                                                                                               <  1:855940/100‑1 (MQ=255)
      aattaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTgg                                                                                               >  2:1045496/1‑100 (MQ=255)
      aattaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTgg                                                                                               <  1:2298765/100‑1 (MQ=255)
       attaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTggg                                                                                              >  2:465447/1‑100 (MQ=255)
        ttaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGt                                                                                             <  1:3588980/100‑1 (MQ=255)
        ttaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGt                                                                                             <  2:1810592/100‑1 (MQ=255)
         taattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATAAATATTGGGtc                                                                                            <  2:1881748/100‑2 (MQ=255)
         taattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGGTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTgg                                                                                               <  2:2172794/97‑1 (MQ=255)
          aattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGtcc                                                                                           <  1:3175326/100‑3 (MQ=39)
          aattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGtcc                                                                                           <  2:1549179/100‑3 (MQ=39)
           attTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCt                                                                                          <  1:3530497/100‑1 (MQ=39)
            tttATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCtt                                                                                         >  2:3156662/1‑100 (MQ=255)
             ttATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCttt                                                                                        <  2:3709589/100‑1 (MQ=38)
             ttATGATTTGGTCGTGGTTAAATCGGAGAACTAATATGTGCAGTTCTTCTTTAGGATGTAAGGTCATATTTTAAATGACCCAATATTTATTGGGTTc                                                                                           <  1:1217459/97‑1 (MQ=255)
              tATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCtttt                                                                                       >  1:1095754/1‑100 (MQ=255)
              tATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCtttt                                                                                       <  2:3736664/100‑1 (MQ=255)
              tATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCtttt                                                                                       <  2:1208695/100‑1 (MQ=38)
              tATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCtttt                                                                                       <  2:3817139/100‑1 (MQ=38)
               aTGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCttttt                                                                                      >  2:2263464/1‑100 (MQ=255)
               aTGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCttttt                                                                                      >  2:278743/1‑100 (MQ=255)
                tGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCttttat                                                                                     <  1:396222/100‑3 (MQ=37)
                 gATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtc                                                                                    <  1:4228652/100‑1 (MQ=255)
                  aTTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtct                                                                                   <  1:1288747/100‑1 (MQ=255)
                   tttGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtctc                                                                                  <  1:3279107/100‑1 (MQ=255)
                      gTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTAt                                                                               <  1:1870150/100‑1 (MQ=255)
                        tCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAAGAACCCAATAAATATTGGGTCCttttatataatgc                                                                              <  2:1577796/100‑10 (MQ=21)
                         cGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCaa                                                                            <  1:33846/100‑1 (MQ=255)
                         cGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCaa                                                                            >  1:4108558/1‑100 (MQ=255)
                         cGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCaa                                                                            <  1:3053106/100‑1 (MQ=255)
                          ggggTTAAATCGGAGAAATAATATGTGAAGTTGTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAAt                                                                           >  2:4053218/1‑100 (MQ=255)
                            ggTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATAc                                                                         >  1:1034485/1‑100 (MQ=255)
                             gTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACt                                                                        >  2:3804274/1‑100 (MQ=255)
                               tAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTAt                                                                      >  1:2839429/1‑100 (MQ=255)
                                aaaTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTAtt                                                                     <  2:2572327/100‑1 (MQ=255)
                                aaaTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTAtt                                                                     <  1:1132049/100‑1 (MQ=255)
                                 aaTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTa                                                                    >  2:3088795/1‑100 (MQ=255)
                                 aaTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTa                                                                    >  2:936666/1‑100 (MQ=255)
                                  aTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAg                                                                   >  2:2500597/1‑100 (MQ=255)
                                   tagcagaAATAATATGTGAAGTTCTTCTATAGTAGGGAAGGTGTTATCTATTCGGAACGTTTATTTATTGGGGTCTTTGTTCTCTATCAAt                                                                           <  1:538213/87‑1 (MQ=255)
                                   tCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGc                                                                  >  2:125215/1‑100 (MQ=255)
                                    cGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCa                                                                 <  1:4195515/100‑1 (MQ=255)
                                     ggAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAg                                                                >  1:2430187/1‑100 (MQ=255)
                                      gagaAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAgg                                                               >  2:1301428/1‑100 (MQ=255)
                                      gagaAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATGGGGTTCTTTTTTCTCTATCAATACTATTAGCAgn                                                               <  1:1597538/100‑2 (MQ=255)
                                         aaaTAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGgag                                                            <  2:4130160/100‑1 (MQ=255)
                                            taataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGata                                                         <  2:3191247/100‑1 (MQ=255)
                                             aataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGatat                                                        <  1:2230311/100‑1 (MQ=255)
                                              atatGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGatata                                                       <  1:2018073/100‑1 (MQ=255)
                                                 tgtgAAGTTCTTCGATAGTATGGAAGGCATTATATAGAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCa                                                    >  2:3006462/1‑100 (MQ=255)
                                                 tgtgAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCa                                                    <  2:2147401/100‑1 (MQ=255)
                                                 tgtgAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCa                                                    >  2:14599/1‑100 (MQ=255)
                                                   tgAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTATCAGGGAGCTATATCa                                                    >  1:598801/1‑98 (MQ=255)
                                                   tgAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCAcc                                                  <  2:2327178/100‑1 (MQ=255)
                                                    gAAGTTCTTCGATAGTATGGAAGGCAGTATATAAAAGGACCCAGTATTTATTGGGTTCTTTTTTCTCTATCAATAATATTAGCAGGAAGatatatga                                                    >  1:3849112/1‑95 (MQ=255)
                                                      aGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCaga                                               >  2:1901931/1‑100 (MQ=255)
                                                       gTTCTTCGATAGTATGGAAGGCATTATATAAAATGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCagag                                              >  2:1974011/1‑100 (MQ=255)
                                                       gTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCagag                                              <  1:3222778/100‑1 (MQ=255)
                                                          cttcGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGttt                                           >  1:218997/1‑100 (MQ=255)
                                                          cttcGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGttt                                           >  2:340834/1‑100 (MQ=255)
                                                             cGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAt                                        >  1:1463273/1‑100 (MQ=255)
                                                             cGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAt                                        >  1:460842/1‑100 (MQ=255)
                                                              gATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtg                                       <  1:2119252/100‑1 (MQ=255)
                                                              gATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtg                                       <  2:3041244/100‑1 (MQ=255)
                                                              gATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtg                                       <  1:4097339/100‑1 (MQ=255)
                                                                tAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtg                                     <  2:3652938/100‑1 (MQ=255)
                                                                 aGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtgt                                    >  2:2715265/1‑100 (MQ=255)
                                                                       gAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAtttt                              >  1:4174038/1‑100 (MQ=255)
                                                                          ggCATTATATAAAAGGACCCAATATTTATTGGGCTCTTTTTTTTATATCAATACTATTAGCAGGGAGATATACTCCCAGAGTTTAATGTGTGGttt                               >  1:1203616/1‑96 (MQ=255)
                                                                           gCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttat                          <  1:2509510/100‑1 (MQ=255)
                                                                             aTTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttattt                        <  2:663628/100‑1 (MQ=255)
                                                                             aTTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttattt                        >  2:3669004/1‑100 (MQ=255)
                                                                             aTTATATAAAAGGAACCAATACTTATTGGGTTCTTTATTCTCTATAAATACTATTAGCAGGGCGCTATATCACCAGAGTTTAATGTTTGATTTtttatt                         >  1:929678/1‑99 (MQ=255)
                                                                              ttATATAAAAGGATCTAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGTGATATATCACCAGAGTTTAATGTGTGATTTtttatt                         <  1:2702307/98‑1 (MQ=255)
                                                                              ttATATAAAAGGACTCTATATTTATTGGGTTTTTTTTTCTCTGTCAATACTTTTAGCAGGGAGATATATCACCAGAGTTTGATGTGTGATTTtttattt                        <  1:2541346/99‑1 (MQ=255)
                                                                                agaaaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtc                     >  1:3092403/5‑100 (MQ=37)
                                                                                 tataAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcg                    <  1:3117359/100‑1 (MQ=255)
                                                                                 gaaaaaaGAACCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTttatttat                      >  1:141818/4‑98 (MQ=38)
                                                                                 gaaaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTttatttat                      >  1:2137358/4‑98 (MQ=37)
                                                                                  aaaaaaGAACCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgt                   >  1:4267983/3‑100 (MQ=38)
                                                                                   taAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgtc                  >  2:3273344/1‑100 (MQ=255)
                                                                                   taAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgtc                  >  1:2876683/1‑100 (MQ=255)
                                                                                   taAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgtc                  <  2:429516/100‑1 (MQ=255)
                                                                                    aaaaGAACCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATGTATGACCAGAGTTTAATGTGTGATTTTTTATTTAtcgtcg                 >  2:3404834/1‑100 (MQ=38)
                                                                                        ggACCCAATATTTATTGGGTTCTTTTTTCTCTATCCATACTATTAGCAGGGAAATCTATCACCACAGTTTAATCTCTGATTTTTTCTTTATCGTCGAccc             >  1:2810291/1‑100 (MQ=255)
                                                                                           aCCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTg           <  2:797002/100‑1 (MQ=255)
                                                                                            cccAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTgg          >  1:2267296/1‑100 (MQ=255)
                                                                                            cccAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTgg          >  2:1042021/1‑100 (MQ=255)
                                                                                            cccAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTgg          <  2:3206714/100‑1 (MQ=255)
                                                                                                atatTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTg      <  2:3538479/100‑1 (MQ=255)
                                                                                                 tatttattGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGt     <  2:2855590/100‑1 (MQ=255)
                                                                                                   tttattGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGttt   <  1:664720/100‑1 (MQ=255)
                                                                                                   tttattGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGttt   >  2:3327844/1‑100 (MQ=255)
                                                                                                    ttactGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGtt    <  1:1993461/98‑1 (MQ=255)
                                                                                                    |                                                                                                   
CTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTA  >  NC_000913/4261260‑4261458

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: