Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,261,358 T→A 15.6% 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,3580TA15.6% 505.4 / 72.4 194intergenic (+52/‑311)yjbM/dusAuncharacterized protein/tRNA‑dihydrouridine synthase A
Reads supporting (aligned to +/- strand):  ref base T (81/81);  new base A (12/18);  total (94/100)
Fisher's exact test for biased strand distribution p-value = 3.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 210 total aligned reads displayed.

ACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTT  >  NC_000913/4261259‑4261457
                                                                                                   |                                                                                                   
 cTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtatt                                                                                                    >  2:846834/1‑100 (MQ=255)
    cTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtatttat                                                                                                 >  2:3107271/1‑100 (MQ=255)
     tAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtatttatt                                                                                                <  2:42187/100‑1 (MQ=255)
     tAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtatttatt                                                                                                >  1:2291658/1‑100 (MQ=255)
      aaaTTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtatttatt                                                                                                <  1:2855395/99‑1 (MQ=255)
       aattaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTgg                                                                                              >  2:1045496/1‑100 (MQ=255)
        attaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTggg                                                                                             >  2:465447/1‑100 (MQ=255)
          taattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATAAATATTGGGtc                                                                                           <  2:1881748/100‑2 (MQ=255)
           aattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTc                                                                                          <  2:1423335/100‑1 (MQ=255)
            attTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATCTATAGGGTTCt                                                                                         <  1:2232864/100‑1 (MQ=255)
             tttATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCtt                                                                                        >  1:3057298/1‑100 (MQ=255)
             tttATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCtt                                                                                        >  2:3156662/1‑100 (MQ=255)
              ttATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCttt                                                                                       <  2:3709589/100‑1 (MQ=38)
              ttATGATTTGGTCGTGGTTAAATCGGAGAACTAATATGTGCAGTTCTTCTTTAGGATGTAAGGTCATATTTTAAATGACCCAATATTTATTGGGTTc                                                                                          <  1:1217459/97‑1 (MQ=255)
              ttATGATGTGTTCGGGGATAAGTCGGAGAAATAAGATGTGAAGGTCTGCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCttt                                                                                       <  1:4034419/100‑1 (MQ=25)
               tATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCtttt                                                                                      <  2:3736664/100‑1 (MQ=255)
               tATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCtttt                                                                                      >  1:1095754/1‑100 (MQ=255)
               tATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCtttt                                                                                      <  2:1208695/100‑1 (MQ=38)
                  gATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtc                                                                                   >  1:2123131/1‑100 (MQ=255)
                  gATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtc                                                                                   <  1:4228652/100‑1 (MQ=255)
                  gATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtc                                                                                   <  2:2263504/100‑1 (MQ=255)
                  gATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtc                                                                                   <  2:2213943/100‑1 (MQ=255)
                   aTTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtct                                                                                  >  2:3011029/1‑100 (MQ=255)
                   aTTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtct                                                                                  <  1:1288747/100‑1 (MQ=255)
                       gTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTAt                                                                              <  1:1870150/100‑1 (MQ=255)
                          cGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCa                                                                            >  2:1565008/1‑99 (MQ=255)
                           ggggTTAAATCGGAGAAATAATATGTGAAGTTGTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAAt                                                                          >  2:4053218/1‑100 (MQ=255)
                           ggggTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAAt                                                                          >  1:287839/1‑100 (MQ=255)
                             ggTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATAc                                                                        >  1:1034485/1‑100 (MQ=255)
                              gTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACt                                                                       >  2:3804274/1‑100 (MQ=255)
                                 aaaTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTAtt                                                                    <  1:1132049/100‑1 (MQ=255)
                                 aaaTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTAtt                                                                    <  2:2572327/100‑1 (MQ=255)
                                  aaTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTa                                                                   >  2:3088795/1‑100 (MQ=255)
                                  aaTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTa                                                                   >  2:398473/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)
                                    tCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGc                                                                 <  2:3289867/100‑1 (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                                                           <  1:4100258/100‑1 (MQ=255)
                                          aaaTAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGgag                                                           <  2:4130160/100‑1 (MQ=255)
                                            ataataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGat                                                         <  1:3202722/100‑1 (MQ=255)
                                             taataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTGTTTTTTCTCTATCAATACTATTAGCAGGGAGata                                                        >  1:2252180/1‑100 (MQ=255)
                                             taataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGata                                                        <  2:3191247/100‑1 (MQ=255)
                                              aataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGatat                                                       <  1:2787418/100‑1 (MQ=255)
                                               atatGTGAACGTCTTCGATAGCATGGAAGGCATTATATAAAAGGACCCAACAGCTATTGGGTTCGTTGTCCTCAATCGATACGATTAGCAGAGAGat                                                         <  1:2624478/97‑1 (MQ=255)
                                                tatGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGatatat                                                     <  1:57194/100‑1 (MQ=255)
                                                  tgtgAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCa                                                   <  2:2147401/100‑1 (MQ=255)
                                                    tgAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCAcc                                                 <  2:1534321/100‑1 (MQ=255)
                                                     gAAGTTGTTCGTTAGTATTGAAGGTATTATATCAAAGGGACTAATATTTATTGCGTTCGTTTTTGTCTATGAATACTATTAGCAGGGAGATATATCTcc                                                 <  1:2141064/99‑1 (MQ=255)
                                                     gAAGTTCTTCGATAGTATGGAAGGCAGTATATAAAAGGACCCAGTATTTATTGGGTTCTTTTTTCTCTATCAATAATATTAGCAGGAAGatatatga                                                   >  1:3849112/1‑95 (MQ=255)
                                                       aGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCaga                                              >  2:1901931/1‑100 (MQ=255)
                                                        gTTCTTCGATAGTATGGAAGGCATTATATAAAATGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCagag                                             >  2:1974011/1‑100 (MQ=255)
                                                           cttcGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGttt                                          >  2:340834/1‑100 (MQ=255)
                                                              cGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAt                                       >  1:460842/1‑100 (MQ=255)
                                                              cGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAt                                       >  1:1463273/1‑100 (MQ=255)
                                                              cGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAt                                       >  1:22180/1‑100 (MQ=255)
                                                               gATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtg                                      <  2:3041244/100‑1 (MQ=255)
                                                               gATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtg                                      <  1:87169/100‑1 (MQ=255)
                                                                 tAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtg                                    <  2:3652938/100‑1 (MQ=255)
                                                                  aGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtgt                                   >  2:2715265/1‑100 (MQ=255)
                                                                    tATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGa                                 >  2:1672860/1‑100 (MQ=255)
                                                                     aTGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAt                                <  2:2839440/100‑1 (MQ=255)
                                                                        gAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAtttt                             >  1:4174038/1‑100 (MQ=255)
                                                                        gAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAttt                              <  1:199700/99‑1 (MQ=255)
                                                                           ggCATTATATAAAAGGACCCAATATTTATTGGGCTCTTTTTTTTATATCAATACTATTAGCAGGGAGATATACTCCCAGAGTTTAATGTGTGGttt                              >  1:1203616/1‑96 (MQ=255)
                                                                              aTTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttattt                       <  2:663628/100‑1 (MQ=255)
                                                                              aTTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttattt                       >  2:3669004/1‑100 (MQ=255)
                                                                                gagaaaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTttatttat                     >  1:3567570/6‑100 (MQ=37)
                                                                                 agaaaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtc                    >  1:3092403/5‑100 (MQ=37)
                                                                                  tataAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcg                   <  1:3117359/100‑1 (MQ=255)
                                                                                  gaaaaaaGAACCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTttatttat                     >  1:141818/4‑98 (MQ=38)
                                                                                   ataAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgt                  >  1:322892/1‑100 (MQ=255)
                                                                                   aaaaaaGAACCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgt                  >  1:4267983/3‑100 (MQ=38)
                                                                                    taAAAGGACCCAATATTTATTGGGTTCTTTTTTTTTTTTTTATACTATTAGCAGGGAGATATATCACAAGAGTTTAATGTGTGATTTTTTATTTAtc                    >  1:2093798/1‑97 (MQ=255)
                                                                                    taAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgtc                 >  2:3273344/1‑100 (MQ=255)
                                                                                    taAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgtc                 >  2:335595/1‑100 (MQ=255)
                                                                                    taAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgtc                 <  2:429516/100‑1 (MQ=255)
                                                                                    taAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgtc                 <  1:737379/100‑1 (MQ=255)
                                                                                    taAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGACATATATCACCGGAGTTTAATGTGTGATTCTTTATTTAtcgtc                 <  2:2887932/100‑1 (MQ=255)
                                                                                    taAAAGGACCAAATATTTATTGGGTTCTTTTTTTTCTATCAATACTATTATCATGGAGATTTATCACCAGAGTTTAATGGTTGATTTTTTATTTAtcgtc                 <  1:350165/100‑1 (MQ=255)
                                                                                    naaaaGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGCTATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgtc                 >  1:3801741/2‑100 (MQ=255)
                                                                                    aaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGCGATATATCACCAGCGTTTAATGTGTGATTTTTTATTTAtcgtc                 >  1:2832759/2‑100 (MQ=37)
                                                                                     aaaaGAACCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATGTATGACCAGAGTTTAATGTGTGATTTTTTATTTAtcgtcg                >  2:3404834/1‑100 (MQ=38)
                                                                                        agaaCCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGaa              >  1:2521817/4‑99 (MQ=255)
                                                                                           aCCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTg          >  2:2172362/1‑100 (MQ=255)
                                                                                           aCCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTg          <  2:797002/100‑1 (MQ=255)
                                                                                            cccAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTgg         <  2:3206714/100‑1 (MQ=255)
                                                                                            cccAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTgg         >  2:1042021/1‑100 (MQ=255)
                                                                                            cccAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTgg         >  1:3288035/1‑100 (MQ=255)
                                                                                            cccAATATTTATTGGGTTCTTGTTTCTCTATCAATACTATTAGCAGGGAGATATATCACTAGAGTTTAATGTGTGATTTTTTATTTATCTTCGAActt           >  1:625932/1‑96 (MQ=255)
                                                                                            cccAATATTTAGTGGGTTCTTTTTTCTCTATCAATACTATTCGCAGGGATATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTgg         >  1:2325087/1‑100 (MQ=255)
                                                                                               aaTATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGAtt      <  2:1464497/100‑1 (MQ=255)
                                                                                               aaTATTTATTGGGTTCTTTTTTCTCTATAAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGAt       <  1:2926220/99‑1 (MQ=255)
                                                                                                atatTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTg     <  2:3538479/100‑1 (MQ=255)
                                                                                                  atttattGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGtt   <  1:3938364/100‑1 (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: