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.

CTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTA  >  NC_000913/4261260‑4261458
                                                                                                    |                                                                                                   
cTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAAtatt                                                                                                     >  2:846834/1‑100 (MQ=255)
  tCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAAtattta                                                                                                   >  1:2917180/1‑100 (MQ=255)
  tCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAAtattta                                                                                                   >  2:3287850/1‑100 (MQ=255)
  tCTAAATAAATCTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGAATGGAAGGCATTATATAAAAGG‑ACCCAAtattta                                                                                                   >  1:2242058/1‑100 (MQ=255)
   cTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAAtatttat                                                                                                  >  2:3107271/1‑100 (MQ=255)
    tAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAAtatttatt                                                                                                 <  2:42187/100‑1 (MQ=255)
     aaaTTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTg                                                                                                >  2:1475096/1‑100 (MQ=255)
      aattaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTgg                                                                                               <  1:855940/100‑1 (MQ=255)
      aattaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTgg                                                                                               >  2:1045496/1‑100 (MQ=255)
      aattaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTgg                                                                                               <  1:2298765/100‑1 (MQ=255)
       attaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTggg                                                                                              >  2:465447/1‑100 (MQ=255)
        ttaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGt                                                                                             <  1:3588980/100‑1 (MQ=255)
        ttaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGt                                                                                             <  2:1810592/100‑1 (MQ=255)
         taattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATAAATATTGGGtc                                                                                            <  2:1881748/100‑2 (MQ=255)
         taattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGGTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTgg                                                                                               <  2:2172794/97‑1 (MQ=255)
          aattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGA‑ACCCAATAAATATTGGGtcc                                                                                           <  1:3175326/100‑3 (MQ=39)
          aattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGA‑ACCCAATAAATATTGGGtcc                                                                                           <  2:1549179/100‑3 (MQ=39)
           attTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGA‑ACCCAATAAATATTGGGTCCt                                                                                          <  1:3530497/100‑1 (MQ=39)
            tttATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCtt                                                                                         >  2:3156662/1‑100 (MQ=255)
             ttATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGA‑ACCCAATAAATATTGGGTCCttt                                                                                        <  2:3709589/100‑1 (MQ=38)
             ttATGATTTGGTCGTGGTTAAATCGGAGAACTAATATGTGCAGTTCTTCTTTAGGATGTAAGGTCATATTTTAAATG‑ACCCAATATTTATTGGGTTc                                                                                           <  1:1217459/97‑1 (MQ=255)
              tATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCtttt                                                                                       >  1:1095754/1‑100 (MQ=255)
              tATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCtttt                                                                                       <  2:3736664/100‑1 (MQ=255)
              tATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGA‑ACCCAATAAATATTGGGTCCtttt                                                                                       <  2:1208695/100‑1 (MQ=38)
              tATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGA‑ACCCAATAAATATTGGGTCCtttt                                                                                       <  2:3817139/100‑1 (MQ=38)
               aTGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCttttt                                                                                      >  2:2263464/1‑100 (MQ=255)
               aTGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCttttt                                                                                      >  2:278743/1‑100 (MQ=255)
                tGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGA‑ACCCAATAAATATTGGGTCCttttat                                                                                     <  1:396222/100‑3 (MQ=37)
                 gATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCTTTTTtc                                                                                    <  1:4228652/100‑1 (MQ=255)
                  aTTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCTTTTTtct                                                                                   <  1:1288747/100‑1 (MQ=255)
                   tttGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCTTTTTtctc                                                                                  <  1:3279107/100‑1 (MQ=255)
                      gTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTAt                                                                               <  1:1870150/100‑1 (MQ=255)
                        tCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAAGAACCCAATAAATATTGGGTCCttttatataatgc                                                                              <  2:1577796/100‑10 (MQ=21)
                         cGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCaa                                                                            <  1:33846/100‑1 (MQ=255)
                         cGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCaa                                                                            >  1:4108558/1‑100 (MQ=255)
                         cGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCaa                                                                            <  1:3053106/100‑1 (MQ=255)
                          ggggTTAAATCGGAGAAATAATATGTGAAGTTGTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCAAt                                                                           >  2:4053218/1‑100 (MQ=255)
                            ggTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATAc                                                                         >  1:1034485/1‑100 (MQ=255)
                             gTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACt                                                                        >  2:3804274/1‑100 (MQ=255)
                               tAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTAt                                                                      >  1:2839429/1‑100 (MQ=255)
                                aaaTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTAtt                                                                     <  2:2572327/100‑1 (MQ=255)
                                aaaTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTAtt                                                                     <  1:1132049/100‑1 (MQ=255)
                                 aaTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTa                                                                    >  2:3088795/1‑100 (MQ=255)
                                 aaTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTa                                                                    >  2:936666/1‑100 (MQ=255)
                                  aTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAg                                                                   >  2:2500597/1‑100 (MQ=255)
                                   tagcagaAATAATATGTGAAGTTCTTCTATAGTAGGGAAGGTGTTATCTATTCGG‑AACGTTTATTTATTGGGGTCTTTGTTCTCTATCAAt                                                                           <  1:538213/87‑1 (MQ=255)
                                   tCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGc                                                                  >  2:125215/1‑100 (MQ=255)
                                    cGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCa                                                                 <  1:4195515/100‑1 (MQ=255)
                                     ggAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAg                                                                >  1:2430187/1‑100 (MQ=255)
                                      gagaAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAgg                                                               >  2:1301428/1‑100 (MQ=255)
                                      gagaAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATGGGGTTCTTTTTTCTCTATCAATACTATTAGCAgn                                                               <  1:1597538/100‑2 (MQ=255)
                                         aaaTAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGgag                                                            <  2:4130160/100‑1 (MQ=255)
                                            taataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGata                                                         <  2:3191247/100‑1 (MQ=255)
                                             aataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGatat                                                        <  1:2230311/100‑1 (MQ=255)
                                              atatGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGatata                                                       <  1:2018073/100‑1 (MQ=255)
                                                 tgtgAAGTTCTTCGATAGTATGGAAGGCATTATATAGAAGG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCa                                                    >  2:3006462/1‑100 (MQ=255)
                                                 tgtgAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCa                                                    <  2:2147401/100‑1 (MQ=255)
                                                 tgtgAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCa                                                    >  2:14599/1‑100 (MQ=255)
                                                   tgAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTATCAGGGAGCTATATCa                                                    >  1:598801/1‑98 (MQ=255)
                                                   tgAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCAcc                                                  <  2:2327178/100‑1 (MQ=255)
                                                    gAAGTTCTTCGATAGTATGGAAGGCAGTATATAAAAGG‑ACCCAGTATTTATTGGGTTCTTTTTTCTCTATCAATAATATTAGCAGGAAGatatatga                                                    >  1:3849112/1‑95 (MQ=255)
                                                      aGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCaga                                               >  2:1901931/1‑100 (MQ=255)
                                                       gTTCTTCGATAGTATGGAAGGCATTATATAAAATG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCagag                                              >  2:1974011/1‑100 (MQ=255)
                                                       gTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCagag                                              <  1:3222778/100‑1 (MQ=255)
                                                          cttcGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGttt                                           >  1:218997/1‑100 (MQ=255)
                                                          cttcGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGttt                                           >  2:340834/1‑100 (MQ=255)
                                                             cGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAt                                        >  1:1463273/1‑100 (MQ=255)
                                                             cGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAt                                        >  1:460842/1‑100 (MQ=255)
                                                              gATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtg                                       <  1:2119252/100‑1 (MQ=255)
                                                              gATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtg                                       <  2:3041244/100‑1 (MQ=255)
                                                              gATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtg                                       <  1:4097339/100‑1 (MQ=255)
                                                                tAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtg                                     <  2:3652938/100‑1 (MQ=255)
                                                                 aGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtgt                                    >  2:2715265/1‑100 (MQ=255)
                                                                       gAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAtttt                              >  1:4174038/1‑100 (MQ=255)
                                                                          ggCATTATATAAAAGG‑ACCCAATATTTATTGGGCTCTTTTTTTTATATCAATACTATTAGCAGGGAGATATACTCCCAGAGTTTAATGTGTGGttt                               >  1:1203616/1‑96 (MQ=255)
                                                                           gCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttat                          <  1:2509510/100‑1 (MQ=255)
                                                                             aTTATATAAAAGG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttattt                        <  2:663628/100‑1 (MQ=255)
                                                                             aTTATATAAAAGG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttattt                        >  2:3669004/1‑100 (MQ=255)
                                                                             aTTATATAAAAGG‑AACCAATACTTATTGGGTTCTTTATTCTCTATAAATACTATTAGCAGGGCGCTATATCACCAGAGTTTAATGTTTGATTTtttatt                         >  1:929678/1‑99 (MQ=255)
                                                                              ttATATAAAAGG‑ATCTAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGTGATATATCACCAGAGTTTAATGTGTGATTTtttatt                         <  1:2702307/98‑1 (MQ=255)
                                                                              ttATATAAAAGG‑ACTCTATATTTATTGGGTTTTTTTTTCTCTGTCAATACTTTTAGCAGGGAGATATATCACCAGAGTTTGATGTGTGATTTtttattt                        <  1:2541346/99‑1 (MQ=255)
                                                                                agaaaaaaGA‑ACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtc                     >  1:3092403/5‑100 (MQ=37)
                                                                                 tataAAAGG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcg                    <  1:3117359/100‑1 (MQ=255)
                                                                                 gaaaaaaGA‑ACCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTttatttat                      >  1:141818/4‑98 (MQ=38)
                                                                                 gaaaaaaGA‑ACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTttatttat                      >  1:2137358/4‑98 (MQ=37)
                                                                                  aaaaaaGA‑ACCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgt                   >  1:4267983/3‑100 (MQ=38)
                                                                                   taAAAGG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgtc                  >  2:3273344/1‑100 (MQ=255)
                                                                                   taAAAGG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgtc                  >  1:2876683/1‑100 (MQ=255)
                                                                                   taAAAGG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgtc                  <  2:429516/100‑1 (MQ=255)
                                                                                    aaaaGA‑ACCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATGTATGACCAGAGTTTAATGTGTGATTTTTTATTTAtcgtcg                 >  2:3404834/1‑100 (MQ=38)
                                                                                        gg‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCCATACTATTAGCAGGGAAATCTATCACCACAGTTTAATCTCTGATTTTTTCTTTATCGTCGAccc             >  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)
                                                                                                    |                                                                                                   
CTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGG‑ACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTA  >  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: ‑