Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,261,358 T→A 18.3% 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,3580TA18.3% 394.0 / 83.2 170intergenic (+52/‑311)yjbM/dusAuncharacterized protein/tRNA‑dihydrouridine synthase A
Reads supporting (aligned to +/- strand):  ref base T (59/78);  new base A (14/17);  total (74/96)
Fisher's exact test for biased strand distribution p-value = 8.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.72e-01

Only 100 of 186 total aligned reads displayed.

CTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTT  >  NC_000913/4261260‑4261457
                                                                                                  |                                                                                                   
  tCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtattta                                                                                                  <  2:442624/100‑1 (MQ=255)
   cTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtatttat                                                                                                 <  2:419669/100‑1 (MQ=255)
   cTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtatttat                                                                                                 <  2:419559/100‑1 (MQ=255)
    tAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTCTGGTAGGCATTATATAAAAGGACCCAAtatttatt                                                                                                >  2:1737922/1‑100 (MQ=255)
    tAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtatttatt                                                                                                <  2:756271/100‑1 (MQ=255)
    tAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtatttatt                                                                                                >  2:3026426/1‑100 (MQ=255)
    tAAATTAATTTATGATTTGTTCGGAGTTAAATCGGAGAAATAATATGTGAAGTGCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtatttatt                                                                                                <  2:2015493/100‑1 (MQ=255)
     aaaTTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTg                                                                                               <  2:32242/100‑1 (MQ=255)
      aattaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTgg                                                                                              >  1:1015660/1‑100 (MQ=255)
          aattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTc                                                                                          <  1:4174494/100‑1 (MQ=255)
          aattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGtcc                                                                                          <  1:59884/100‑3 (MQ=39)
           attTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCt                                                                                         <  1:1949320/100‑1 (MQ=38)
           attTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCt                                                                                         <  1:2005294/100‑1 (MQ=38)
             ttATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCttt                                                                                       <  2:1476755/100‑1 (MQ=255)
             ttATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCttt                                                                                       <  1:2073891/100‑1 (MQ=38)
             ttATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCttt                                                                                       <  1:4156941/100‑1 (MQ=38)
              tATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCtttt                                                                                      >  1:3312692/1‑100 (MQ=255)
               aTGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCttttt                                                                                     <  2:970778/100‑1 (MQ=255)
               aTGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCttttt                                                                                     >  1:4073871/1‑100 (MQ=255)
               aTGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCttttt                                                                                     <  1:2447337/100‑1 (MQ=255)
               aTGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCtttta                                                                                     <  1:3794341/100‑2 (MQ=38)
                tGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCtttttt                                                                                    >  2:1550629/1‑100 (MQ=255)
                tGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCttttat                                                                                    <  1:1503696/100‑3 (MQ=37)
                tGATTTGTTCGGAGTTAAATCTGAGAAATAATATGTGAAGTCCTTCTATAGTATGGAAGTCATTATATAAAAGAACCCAATAAATATTGGGTCCttttat                                                                                    <  2:1992967/100‑3 (MQ=21)
                 gATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtc                                                                                   >  2:3435813/1‑100 (MQ=255)
                 gATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtc                                                                                   <  1:2029914/100‑1 (MQ=255)
                 gATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtc                                                                                   <  1:3878494/100‑1 (MQ=255)
                 gATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtc                                                                                   >  1:3231229/1‑100 (MQ=255)
                 gATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtc                                                                                   <  2:2140008/100‑1 (MQ=255)
                 gATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCttttata                                                                                   <  2:2387397/100‑4 (MQ=37)
                 gATTTGTTCAGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtc                                                                                   <  2:2296949/100‑1 (MQ=255)
                  aTTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtct                                                                                  >  2:719632/1‑100 (MQ=255)
                  aTTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtct                                                                                  <  2:224589/100‑1 (MQ=255)
                   tttGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtctc                                                                                 >  2:31147/1‑100 (MQ=255)
                    ttGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTctct                                                                                <  1:3089398/100‑1 (MQ=255)
                    ttGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTctct                                                                                <  1:1348860/100‑1 (MQ=255)
                       ttCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATc                                                                             >  2:2173277/1‑100 (MQ=255)
                       ttCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATc                                                                             <  1:1666478/100‑1 (MQ=255)
                       ttCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATc                                                                             <  1:2174292/100‑1 (MQ=255)
                        tCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCa                                                                            <  1:974987/100‑1 (MQ=255)
                                 aaTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTa                                                                   <  1:1832345/100‑1 (MQ=255)
                                 aaTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTa                                                                   <  1:858015/100‑1 (MQ=255)
                                  aTCGGATAAATAATATGTCAAGTTCTCCGATAGTATGGAAGCCATTATATAAAAGGACCTAATATTTATTCGGTTCTTTTTTCTCTATCAATACTATTAg                                                                  <  2:3133056/100‑1 (MQ=255)
                                   tCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGc                                                                 >  2:2627327/1‑100 (MQ=255)
                                    cGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCa                                                                <  2:3737794/100‑1 (MQ=255)
                                    cGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCa                                                                >  1:2743552/1‑100 (MQ=255)
                                    cGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCa                                                                >  1:3271405/1‑100 (MQ=255)
                                       agaAATAATATGTGAATTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAggg                                                             <  2:1760389/100‑1 (MQ=255)
                                       agaAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAggg                                                             <  1:355690/100‑1 (MQ=255)
                                         aaaTAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGgag                                                           <  1:579875/100‑1 (MQ=255)
                                         aaaTAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGgag                                                           <  1:2627334/100‑1 (MQ=255)
                                         aaaTAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGgag                                                           <  2:4133997/100‑1 (MQ=255)
                                          aataatATGTGAAGTTATTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTATTTTATATATAAATAATATTAGCAGGgaga                                                          >  2:365890/1‑100 (MQ=255)
                                          aataatATGCGAAGTTCTTCCATAGTATTTAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTCTTCTCTATCAATACTATTAGCAGGtag                                                           <  2:2829931/99‑1 (MQ=255)
                                           ataataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGat                                                         <  1:1164741/100‑1 (MQ=255)
                                           ataataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGat                                                         <  2:4148577/100‑1 (MQ=255)
                                            taataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGata                                                        >  2:235049/1‑100 (MQ=255)
                                            taataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGata                                                        >  2:845803/1‑100 (MQ=255)
                                            taataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGata                                                        >  1:2292305/1‑100 (MQ=255)
                                             aagatatgtaGTTCTTCTATATTTTTGAGGTCATTATATAAAAGGACACAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGat                                                         <  2:3784769/89‑1 (MQ=255)
                                               tatGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGatatat                                                     <  2:4090062/100‑1 (MQ=255)
                                                atGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATc                                                    <  2:1678199/100‑1 (MQ=255)
                                                atGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAGATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATc                                                    <  1:338516/100‑1 (MQ=255)
                                                 tgtgAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCa                                                   >  1:2325055/1‑100 (MQ=255)
                                                   tgAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCAcc                                                 >  2:2728102/1‑100 (MQ=255)
                                                          cttgtaTAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCagag                                             <  2:3473798/92‑1 (MQ=255)
                                                           ttcGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTa                                         <  1:292276/100‑1 (MQ=255)
                                                            tcGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTaa                                        <  1:3738065/100‑1 (MQ=255)
                                                                tAGTATGGAAGTCATTATATAAAAAGAACCAAGAGTGAGTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGGGCTTAAtg                                      <  2:2905952/98‑1 (MQ=255)
                                                                  gTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtgtg                                  >  1:351874/1‑100 (MQ=255)
                                                                    aTGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAt                                >  1:3704789/1‑100 (MQ=255)
                                                                    aTGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAt                                >  1:1150447/1‑100 (MQ=255)
                                                                    aTGGAAGGCATTAGATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAt                                <  1:167632/100‑1 (MQ=255)
                                                                         aGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAtttttt                           <  2:3610090/100‑1 (MQ=255)
                                                                         aGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAtttttt                           <  1:4173577/100‑1 (MQ=255)
                                                                             tagagaaaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGCTTAATGTGTGATTTtttattt                       >  1:2412879/8‑100 (MQ=25)
                                                                             aTTATATAAAAGGACCCAATATTTATTGTGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttattt                       <  1:2494185/100‑1 (MQ=255)
                                                                             aTTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttattt                       <  1:110970/100‑1 (MQ=255)
                                                                              ttATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttattta                      <  2:4016533/100‑1 (MQ=255)
                                                                              ttATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttattta                      <  2:949875/100‑1 (MQ=255)
                                                                               gagaaaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCCGAGTTTAATGTGTGATTTTctatttat                     >  2:280276/6‑100 (MQ=25)
                                                                               gagaaaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTttatttat                     >  2:124514/6‑100 (MQ=25)
                                                                                    aaaaGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgtcg                >  1:706319/1‑100 (MQ=255)
                                                                                    aaaaGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgtcg                >  2:2567429/1‑100 (MQ=255)
                                                                                    aaaaGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgtcg                <  1:2056035/100‑1 (MQ=255)
                                                                                    aaaaGAACCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgtcg                >  1:2467552/1‑100 (MQ=39)
                                                                                     aaaGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGa               >  1:543024/1‑100 (MQ=255)
                                                                                     aaaGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGa               >  1:3773765/1‑100 (MQ=255)
                                                                                     aaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgtc                 >  2:1425255/1‑98 (MQ=38)
                                                                                       agaaCCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAAc             >  2:2993897/4‑100 (MQ=255)
                                                                                        ggACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtagtcg                >  2:3084075/1‑96 (MQ=255)
                                                                                          aCCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCATAGTTTAATGTGTGATTTTTTATTTATCGTCGAACTTg          >  2:2956184/1‑100 (MQ=255)
                                                                                           cccAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTgg         >  1:1882486/1‑100 (MQ=255)
                                                                                             cAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGAt       >  2:2177547/1‑100 (MQ=255)
                                                                                              aaTATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGAtt      >  2:2451736/1‑100 (MQ=255)
                                                                                                 atttattGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGtt   >  1:444732/1‑100 (MQ=255)
                                                                                                  tttattGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGttt  <  2:3958994/100‑1 (MQ=255)
                                                                                                  tttattGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGttt  >  1:3959021/1‑100 (MQ=255)
                                                                                                  tttattGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGttt  <  1:2728109/100‑1 (MQ=255)
                                                                                                  tttattGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGttt  <  1:1014056/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
CTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTT  >  NC_000913/4261260‑4261457

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: