Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,261,349 G→A 16.1% intergenic (+43/‑320) 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,3490GA16.1% 395.2 / 62.4 168intergenic (+43/‑320)yjbM/dusAuncharacterized protein/tRNA‑dihydrouridine synthase A
Reads supporting (aligned to +/- strand):  ref base G (58/82);  new base A (11/16);  total (69/99)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-02

Only 100 of 171 total aligned reads displayed.

CTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACC  >  NC_000913/4261250‑4261447
                                                                                                   |                                                                                                  
  gctcttccgatCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGAc                                                                                                  <  1:483076/90‑1 (MQ=255)
    cAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGAccc                                                                                                <  2:3829613/100‑1 (MQ=255)
    cAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTATTCGATAGTATGGAAGGCATTATATAAAAGGAccc                                                                                                >  1:398919/1‑100 (MQ=255)
      gACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCaa                                                                                              >  1:1774556/1‑100 (MQ=255)
       acacTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAat                                                                                             <  1:1047313/100‑1 (MQ=255)
            tCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtattta                                                                                        <  2:442624/100‑1 (MQ=255)
             cTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtatttat                                                                                       <  2:419559/100‑1 (MQ=255)
              tAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtatttatt                                                                                      >  2:3026426/1‑100 (MQ=255)
               aaaTTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTg                                                                                     <  2:32242/100‑1 (MQ=255)
                aattaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTgg                                                                                    >  1:2743006/1‑100 (MQ=255)
                aattaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTgg                                                                                    >  1:1015660/1‑100 (MQ=255)
                  ttaattTATGATTTGTTCGGTATTAAATCGGAAAAATAATATGTGAAGTTCTTCGATAATATGCAAGGCATTATATAAAAGTACCCAATATTTATTGGGt                                                                                  <  1:2054174/100‑1 (MQ=255)
                    aattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTc                                                                                <  1:4174494/100‑1 (MQ=255)
                    aattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGtcc                                                                                <  1:59884/100‑3 (MQ=39)
                     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                                                                            <  2:725324/100‑1 (MQ=255)
                        tATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCtttt                                                                            >  1:3312692/1‑100 (MQ=255)
                         aTGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCttttt                                                                           >  1:4073871/1‑100 (MQ=255)
                         aTGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCttttt                                                                           <  2:970778/100‑1 (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:1875315/100‑3 (MQ=37)
                          tGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCttttat                                                                          <  1:1503696/100‑3 (MQ=37)
                           gATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtc                                                                         <  1:2396988/100‑1 (MQ=255)
                           gATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtc                                                                         <  1:3878494/100‑1 (MQ=255)
                           gATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtc                                                                         <  1:2029914/100‑1 (MQ=255)
                           gATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtc                                                                         <  1:2457787/100‑1 (MQ=255)
                           gATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCttttata                                                                         <  2:2387397/100‑4 (MQ=37)
                           gATTTGTTCAGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtc                                                                         <  2:2296949/100‑1 (MQ=255)
                            aTTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtct                                                                        <  2:224589/100‑1 (MQ=255)
                            aTTTGTTCGGGGTTAAATCGGAAAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCttttat                                                                          <  2:1994342/98‑3 (MQ=37)
                             tttGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtctc                                                                       >  2:31147/1‑100 (MQ=255)
                             tttGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtctc                                                                       >  1:2321129/1‑100 (MQ=255)
                              ttGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTctct                                                                      <  1:1348860/100‑1 (MQ=255)
                              ttGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTctct                                                                      <  1:3089398/100‑1 (MQ=255)
                                 ttCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATc                                                                   <  1:2174292/100‑1 (MQ=255)
                                 ttCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATc                                                                   >  2:2173277/1‑100 (MQ=255)
                                 ttCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATc                                                                   <  1:1666478/100‑1 (MQ=255)
                                  tCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCa                                                                  <  1:974987/100‑1 (MQ=255)
                                   cGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCaa                                                                 >  2:3440536/1‑100 (MQ=255)
                                    ggggTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAAt                                                                <  2:2865448/100‑1 (MQ=255)
                                     gggTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAAt                                                                >  2:2472966/1‑99 (MQ=255)
                                        ttAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTa                                                            <  1:370599/100‑1 (MQ=255)
                                           aaTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTa                                                         <  2:2159474/100‑1 (MQ=255)
                                           aaTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTa                                                         <  1:858015/100‑1 (MQ=255)
                                           aaTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTa                                                         >  2:1716148/1‑100 (MQ=255)
                                             tCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGc                                                       >  2:2627327/1‑100 (MQ=255)
                                              cGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGc                                                       <  2:3013640/99‑1 (MQ=255)
                                              cGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCa                                                      <  2:3737794/100‑1 (MQ=255)
                                              cGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCa                                                      >  1:3271405/1‑100 (MQ=255)
                                              cGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCa                                                      >  1:2743552/1‑100 (MQ=255)
                                                 agaAATAATATGTGAATTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAggg                                                   <  2:1760389/100‑1 (MQ=255)
                                                   aaaTAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGgag                                                 <  2:4133997/100‑1 (MQ=255)
                                                   aaaTAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGgag                                                 <  1:2627334/100‑1 (MQ=255)
                                                   aaaTAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGgag                                                 <  1:1161683/100‑1 (MQ=255)
                                                   aaaTAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGgag                                                 <  1:579875/100‑1 (MQ=255)
                                                   aaaTAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGgag                                                 >  1:3551074/1‑100 (MQ=255)
                                                    aatcctatGTGACGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGgag                                                 <  1:188697/94‑1 (MQ=255)
                                                    aataatATGTGAAGTTATTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTATTTTATATATAAATAATATTAGCAGGgaga                                                >  2:365890/1‑100 (MQ=255)
                                                    aataatATGCGAAGTTCTTCCATAGTATTTAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTCTTCTCTATCAATACTATTAGCAGGtag                                                 <  2:2829931/99‑1 (MQ=255)
                                                     ataataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGat                                               <  2:4148577/100‑1 (MQ=255)
                                                     ataataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGat                                               <  1:1164741/100‑1 (MQ=255)
                                                      taataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGata                                              >  2:235049/1‑100 (MQ=255)
                                                      taataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGata                                              >  1:2292305/1‑100 (MQ=255)
                                                      taataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGata                                              >  2:845803/1‑100 (MQ=255)
                                                      taataTGTGAAATTCTTCGATAGTAATGGAGGAATTATAAAAAAGGATCCAAGATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCa                                                      <  2:3531812/92‑1 (MQ=255)
                                                        atatGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGatata                                            <  2:1628577/100‑1 (MQ=255)
                                                         tatGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGatatat                                           <  2:4090062/100‑1 (MQ=255)
                                                          atGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATc                                          <  2:43751/100‑1 (MQ=255)
                                                          atGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAGATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATc                                          <  1:338516/100‑1 (MQ=255)
                                                            gtgAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCAc                                        >  1:4160593/1‑100 (MQ=255)
                                                             tgAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCAcc                                       >  2:2728102/1‑100 (MQ=255)
                                                                     ttcGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTa                               <  1:292276/100‑1 (MQ=255)
                                                                         aTAGTATGGAAGGCATTATATAAAACGACCCCATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtc                            <  2:831331/99‑2 (MQ=255)
                                                                          tAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtg                          <  1:3115564/100‑1 (MQ=255)
                                                                            gTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtgtg                        <  1:3423184/100‑1 (MQ=255)
                                                                              aTGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAt                      >  1:3704789/1‑100 (MQ=255)
                                                                                 gAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAtttt                   >  1:1147613/1‑100 (MQ=255)
                                                                                  aaGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAttttt                  >  1:734931/1‑100 (MQ=255)
                                                                                   aGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAtttttt                 <  2:3610090/100‑1 (MQ=255)
                                                                                      cATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttatt              <  2:3790286/100‑1 (MQ=255)
                                                                                      cATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttat               <  2:3162407/99‑1 (MQ=255)
                                                                                       aTTATATAAAAGGACCCAATATTTATTGTGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttattt             <  1:2494185/100‑1 (MQ=255)
                                                                                        ttATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttattta            <  2:4016533/100‑1 (MQ=255)
                                                                                        ttATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttattta            <  2:949875/100‑1 (MQ=255)
                                                                                         gagaaaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTttatttat           >  2:124514/6‑100 (MQ=25)
                                                                                             taAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAAAACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTGATTTAtcgtc       <  2:2950786/100‑1 (MQ=255)
                                                                                              aaaaGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgtcg      >  2:2567429/1‑100 (MQ=255)
                                                                                              aaaaGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgtcg      >  1:706319/1‑100 (MQ=255)
                                                                                              aaaaGAACCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgtc       >  2:2113880/1‑99 (MQ=39)
                                                                                               aaaGGACCCAATATTTATTGTGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATTTGTGATTTTTTATTTATCGTCGa     <  2:1742813/100‑1 (MQ=255)
                                                                                               aaaGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGa     >  1:543024/1‑100 (MQ=255)
                                                                                               aaaGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGa     >  1:3773765/1‑100 (MQ=255)
                                                                                               aaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgtc       >  2:1425255/1‑98 (MQ=38)
                                                                                                 aGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAAc   >  1:1422123/1‑100 (MQ=255)
                                                                                                  ggACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtagtcg      >  2:3084075/1‑96 (MQ=255)
                                                                                                   |                                                                                                  
CTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACC  >  NC_000913/4261250‑4261447

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: