Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,261,359 T→A 18.7% 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,3590TA18.7% 498.0 / 96.8 208intergenic (+53/‑310)yjbM/dusAuncharacterized protein/tRNA‑dihydrouridine synthase A
Reads supporting (aligned to +/- strand):  ref base T (76/92);  new base A (19/20);  total (95/113)
Fisher's exact test for biased strand distribution p-value = 7.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 218 total aligned reads displayed.

CTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTA  >  NC_000913/4261260‑4261458
                                                                                                    |                                                                                                   
cTTCTAAATTAATTTATGATTTGTTCGGGGTAAAATCGGAGAAATAATATGTGAAGTTCTACGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAAtatt                                                                                                     <  1:1387365/100‑1 (MQ=255)
  tCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAAtattta                                                                                                   <  1:4181063/100‑1 (MQ=255)
     aaaTTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTg                                                                                                <  1:1038170/100‑1 (MQ=255)
      aattaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTgg                                                                                               <  1:516957/100‑1 (MQ=255)
       attcattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAAGCATTATAT‑AAAACGACTCAAtatttcatgg                                                                                               >  1:3432122/1‑94 (MQ=255)
       attaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTggg                                                                                              >  2:1758350/1‑100 (MQ=255)
          aattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGAACCCAATAAATATTGGGtcc                                                                                           <  2:687256/100‑3 (MQ=39)
              tATGGTTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCtttt                                                                                       <  1:240231/100‑1 (MQ=255)
               aTGATTTGTTCGGGGTTAAATCGGAGAGATAATATGTGAAGTTCTTCGATAGTATGGAAGGCTTTATAT‑AAAAGAACCCAATAAATATTGGGTCCtttta                                                                                      <  1:2351147/100‑2 (MQ=37)
                tGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCtttttt                                                                                     <  2:3799504/100‑1 (MQ=255)
                 gATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCTATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTtc                                                                                    <  1:1144095/100‑1 (MQ=255)
                 gATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTtc                                                                                    <  1:660810/100‑1 (MQ=255)
                 gATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTtc                                                                                    <  2:2138617/100‑1 (MQ=255)
                   tttGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTtctc                                                                                  <  2:840665/100‑1 (MQ=255)
                   tttGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTtctc                                                                                  <  2:1687325/100‑1 (MQ=255)
                    ttGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGAACCCAATAAATATTGGGTCCttttatataa                                                                                 <  2:3533432/100‑7 (MQ=25)
                     tGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTa                                                                                <  1:3894819/100‑1 (MQ=255)
                     tGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTa                                                                                <  2:3844237/100‑1 (MQ=255)
                     tGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTa                                                                                <  1:1223349/100‑1 (MQ=255)
                      gTTCGGGGTTAAATCGGAGTAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTAt                                                                               >  2:1981689/1‑100 (MQ=255)
                      gTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTAt                                                                               >  2:2829015/1‑100 (MQ=255)
                      gTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTAt                                                                               >  2:459722/1‑100 (MQ=255)
                       ttCGGTGTTAAATCGGAGAAATAATATGTAAAGGTATTCGATAGTATGGAAGTCATTATAT‑AAAAGAACCCAATAAATATTGGGTCCttttatataatg                                                                               <  1:1861033/99‑9 (MQ=18)
                       ttCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATc                                                                              <  2:2434754/100‑1 (MQ=255)
                        tCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAAGAACCCAATAAATATTGGGTCCttttatataatgc                                                                              <  1:3452187/100‑10 (MQ=18)
                        tCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGAACCCAATAAATATTGGGTCCttttatataatgcc                                                                             <  1:3809902/100‑11 (MQ=21)
                           gggTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATa                                                                          >  1:3057428/1‑100 (MQ=255)
                           gggTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATa                                                                          <  2:3210406/100‑1 (MQ=255)
                            ggTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATAc                                                                         <  1:80365/100‑1 (MQ=255)
                            ggTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATAc                                                                         >  2:3932481/1‑100 (MQ=255)
                                 aaTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTa                                                                    >  1:3740934/1‑100 (MQ=255)
                                  aTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAg                                                                   >  2:3067505/1‑100 (MQ=255)
                                   tCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAg                                                                   >  2:881247/1‑99 (MQ=255)
                                   tCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGc                                                                  <  1:4015327/100‑1 (MQ=255)
                                   tCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGc                                                                  >  2:2218667/1‑100 (MQ=255)
                                    cGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTAtcgaag                                                                 >  2:277656/1‑95 (MQ=255)
                                    cGGAGAAATAATATGTGAAGTTATTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCa                                                                 >  1:2449603/1‑100 (MQ=255)
                                     ggAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAg                                                                <  2:4070836/100‑1 (MQ=255)
                                      gagaAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAgg                                                               <  2:1778312/100‑1 (MQ=255)
                                          aataatATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGgaga                                                           <  1:4043245/100‑1 (MQ=255)
                                           ataataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGat                                                          >  1:284280/1‑100 (MQ=255)
                                           ataataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGat                                                          >  2:3184598/1‑100 (MQ=255)
                                            taataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGata                                                         >  2:2511996/1‑100 (MQ=255)
                                            taataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGata                                                         <  2:2024437/100‑1 (MQ=255)
                                            taataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGata                                                         <  2:793704/100‑1 (MQ=255)
                                             aataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGatat                                                        <  2:3965149/100‑1 (MQ=255)
                                              atatGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGatata                                                       <  1:168043/100‑1 (MQ=255)
                                              atatGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGatata                                                       <  2:2550986/100‑1 (MQ=255)
                                                 tttgAAGTTCTTCAATAGTATGGAAGGCATTATAT‑AAAAGGAGCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCa                                                    >  2:2983304/3‑100 (MQ=255)
                                                  gtgAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCAc                                                   >  2:2267460/1‑100 (MQ=255)
                                                   tgAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCAcc                                                  <  2:3379784/100‑1 (MQ=255)
                                                    gAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCa                                                 >  1:4237871/1‑100 (MQ=255)
                                                       gTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCagag                                              >  2:4221390/1‑100 (MQ=255)
                                                          cttcGATAGTATGAAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGttt                                           <  2:3483225/100‑1 (MQ=255)
                                                               aTAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATCCTATTAGCCGTGAGATATATCACCAGAGTTTAAtgt                                      >  2:56475/1‑100 (MQ=255)
                                                                tAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtg                                     <  2:2782944/100‑1 (MQ=255)
                                                                 aGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtgt                                    <  2:1954336/100‑1 (MQ=255)
                                                                 aGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtgt                                    <  1:1804309/100‑1 (MQ=255)
                                                                 aGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtgt                                    <  1:4058218/100‑1 (MQ=255)
                                                                  gTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtgtg                                   >  2:2187681/1‑100 (MQ=255)
                                                                  gTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtgtg                                   >  1:3261667/1‑100 (MQ=255)
                                                                    aTGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAt                                 >  2:1003329/1‑100 (MQ=255)
                                                                    aTGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAt                                 <  2:1501025/100‑1 (MQ=255)
                                                                      ggAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAttt                               >  1:457665/1‑100 (MQ=255)
                                                                        aaGGCATTATTT‑ATCAGGACCCAATATGTATTGGGTTCTTTTTTCTCTATCACTACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAttttt                             <  1:320650/100‑1 (MQ=255)
                                                                         aGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAtttttt                            >  2:768181/1‑100 (MQ=255)
                                                                          ggCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGTAGGGAGATATATCACCAGAGTTTAATGTGTGAtttttt                            <  1:2702476/99‑1 (MQ=255)
                                                                          ggCATTATAT‑AAAAGGACCCAAGATTTAGTGGGTTCTTTTTCCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTttta                           <  2:1435896/100‑1 (MQ=255)
                                                                           tgatagagaaaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTttta                           >  1:2499827/11‑100 (MQ=21)
                                                                           gCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttat                          <  2:128844/100‑1 (MQ=255)
                                                                             aTTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttattt                        <  2:3587008/100‑1 (MQ=255)
                                                                             aTTATAT‑AAAAGGACCCAATATATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCATCAGAGTTTAATGTGTGATGTttgattt                        >  1:1533882/1‑100 (MQ=255)
                                                                              tagagaaaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttattt                        >  1:192596/8‑100 (MQ=25)
                                                                               tatat‑aAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTttatttat                      >  2:1920768/1‑100 (MQ=255)
                                                                                gagaaaaaaGAACCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGATTTTAATGTGTGATTTTttatttat                      >  1:1094815/6‑100 (MQ=37)
                                                                                atat‑aAAAGGACCCAATATTTATTGGGCTCTTTCTTCTCTATAAAAACTATTAACAGGGCGATATATCACCCGCCTTCAACGTGTGATTTTTTCTTtctc                     >  1:1495080/1‑100 (MQ=255)
                                                                                  gaaaaaaGAACCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcg                    >  2:1002464/4‑100 (MQ=38)
                                                                                   t‑aAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgtc                  <  2:3934822/100‑1 (MQ=255)
                                                                                    aaaaaGAACCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCAGCAGAGTTTAATGTGTGATTTTTTGTTTAtcgtc                  >  2:2594037/2‑100 (MQ=38)
                                                                                    aaaaaGAACCCAATAAAAATTGGGTTGTTTTTTCTCTATCAATGCTTTTAACTGGGAGATATATTAGTAGAGTTGGATGTGAGATTTgttaattagttta                  >  1:4118727/2‑95 (MQ=255)
                                                                                        agaaCCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAAc              >  2:1633751/4‑100 (MQ=255)
                                                                                         ggACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAAcc             >  1:221961/1‑100 (MQ=255)
                                                                                         gaaCCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAAcc             >  2:3354621/3‑100 (MQ=255)
                                                                                         gaaCCCAATAAATATTGGGTTC‑TTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCt            >  2:1603767/3‑100 (MQ=255)
                                                                                           aCCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTg           >  2:4153577/1‑100 (MQ=255)
                                                                                           aCCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTg           <  1:2792268/100‑1 (MQ=255)
                                                                                           aCCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTg           <  1:2113552/100‑1 (MQ=255)
                                                                                           aCCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTg           <  2:1576634/100‑1 (MQ=255)
                                                                                            cccAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTgg          <  2:2885465/100‑1 (MQ=255)
                                                                                             ccAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGa         >  1:591733/1‑100 (MQ=255)
                                                                                             ccAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGa         <  1:2478593/100‑1 (MQ=255)
                                                                                             ccAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGa         >  2:541320/1‑100 (MQ=255)
                                                                                             ccAATAATTATTGAGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGCTATATCACCAGAGATTAATGTGTGATTTTTTATTTATCGTCGAATCt            >  1:1728418/1‑97 (MQ=255)
                                                                                              cAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTg           >  1:1249463/1‑97 (MQ=255)
                                                                                              cAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGAt        <  1:1556131/100‑1 (MQ=255)
                                                                                                atatTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGAtt       >  1:4222363/1‑99 (MQ=255)
                                                                                                  atttattGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGtt    <  1:3029050/100‑1 (MQ=255)
                                                                                                  atttattGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGtt    <  1:1002491/100‑1 (MQ=255)
                                                                                                   tttattGGGTTCTTTTTTCTCTATCAATACTATTAGCTGTGATATATATCACGAGAGTTTAATGTTTGATTTTTTATTTATCGTCTAACCTGGATTGttt   <  1:970174/100‑1 (MQ=255)
                                                                                                   tttattGGGTTCTTTTTTCTCTATCAATACTATGAGCAGGGAGATATATCTCCAGAGTTTAATGTGTTATTTGTTATTTATCGTCGAACCTTTATTGttt   <  1:361572/100‑1 (MQ=255)
                                                                                                    |                                                                                                   
CTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTA  >  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: ‑