Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,261,349 G→A 16.9% 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.9% 417.6 / 64.6 177intergenic (+43/‑320)yjbM/dusAuncharacterized protein/tRNA‑dihydrouridine synthase A
Reads supporting (aligned to +/- strand):  ref base G (68/79);  new base A (14/16);  total (82/95)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.61e-01

Only 100 of 179 total aligned reads displayed.

CTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAAC  >  NC_000913/4261250‑4261446
                                                                                                    |                                                                                                 
cTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAgg                                                                                                   >  1:448200/1‑100 (MQ=255)
  tACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGAc                                                                                                 >  1:2857859/1‑100 (MQ=255)
     aGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCa                                                                                              <  1:3632054/100‑1 (MQ=255)
      gACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCaa                                                                                             >  1:3199479/1‑100 (MQ=255)
      gACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCaa                                                                                             >  1:2928590/1‑100 (MQ=255)
          cTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtatt                                                                                         >  1:2943938/1‑100 (MQ=255)
          cTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtatt                                                                                         >  2:2830786/1‑100 (MQ=255)
          cTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtatt                                                                                         >  1:2027801/1‑100 (MQ=255)
            tCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtattta                                                                                       <  1:3368214/100‑1 (MQ=255)
             cTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtatttat                                                                                      >  1:2745496/1‑100 (MQ=255)
              tAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtatttatt                                                                                     <  2:1162168/100‑1 (MQ=255)
              tAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtatttatt                                                                                     <  1:2213074/100‑1 (MQ=255)
               aaaTTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTg                                                                                    >  2:2589567/1‑100 (MQ=255)
               aaaTTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTg                                                                                    <  1:3378864/100‑1 (MQ=255)
                aattaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTgg                                                                                   >  2:3202285/1‑100 (MQ=255)
                aattaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTgg                                                                                   <  1:2505251/100‑1 (MQ=255)
                   taaatTATTATTTGTTCGGGCTTAAATCGGAGAAATAATATGTGAAGTTATACGATAATATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGt                                                                                 <  2:1774321/99‑1 (MQ=255)
                    aattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTc                                                                               >  1:2271817/1‑100 (MQ=255)
                    aattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATAAATATTGGGtc                                                                                <  1:1800242/99‑2 (MQ=255)
                     attTATTATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGGAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCt                                                                              <  2:3593181/100‑1 (MQ=38)
                     attTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCt                                                                              >  1:1506636/1‑100 (MQ=255)
                      tttATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCtt                                                                             <  1:2073321/100‑1 (MQ=38)
                       ttATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCttt                                                                            >  2:559479/1‑100 (MQ=255)
                       ttATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCttt                                                                            <  1:2702355/100‑1 (MQ=38)
                        tATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCtttt                                                                           <  1:1343867/100‑1 (MQ=255)
                         aTGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCtttta                                                                          <  2:2042417/100‑2 (MQ=38)
                           gatctGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAAGAACCCAATAAATATTGGGTCCttttat                                                                         <  2:490971/96‑3 (MQ=25)
                           gATTTGTTCGGGGTTAAATCGGAGAAATAATATTTTAAGTTTTTCGATAGTATGGAAGGCATTATATAAAATGACCCAATATTTATTGGGTTCTTTTTtc                                                                        <  2:694334/100‑1 (MQ=255)
                           gATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtc                                                                        <  1:2552133/100‑1 (MQ=255)
                           gATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtc                                                                        <  2:3069095/100‑1 (MQ=255)
                           gATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCttttata                                                                        <  1:1948462/100‑4 (MQ=37)
                           gATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCttttata                                                                        <  2:3195260/100‑4 (MQ=37)
                            aTTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtct                                                                       <  1:2315119/100‑1 (MQ=255)
                            aTTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtct                                                                       >  2:2356571/1‑100 (MQ=255)
                             tttGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtctc                                                                      <  1:2362039/100‑1 (MQ=255)
                              ttGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTctct                                                                     <  1:428028/100‑1 (MQ=255)
                              ttGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTctct                                                                     <  1:710422/100‑1 (MQ=255)
                               tGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTa                                                                    <  1:2087361/100‑1 (MQ=255)
                                gTTCGGGGTTAAATCGGAGAAATAATATGTGAAGGTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAGATTAATTGGGTTCTTTTTTCTCTAt                                                                   >  2:887559/1‑100 (MQ=255)
                                 ttCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATc                                                                  <  1:2061491/100‑1 (MQ=255)
                                    ggggTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAAt                                                               >  1:2416142/1‑100 (MQ=255)
                                    ggggTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAAt                                                               >  1:1514257/1‑100 (MQ=255)
                                     gggTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATa                                                              <  2:1399641/100‑1 (MQ=255)
                                       gTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACt                                                            >  2:1241016/1‑100 (MQ=255)
                                           aaTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTa                                                        >  1:2701840/1‑100 (MQ=255)
                                              cGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCa                                                     >  1:3098652/1‑100 (MQ=255)
                                              cGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCa                                                     >  1:500144/1‑100 (MQ=255)
                                              cGGAGAAATAATATGTGAAGCTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCa                                                     >  2:3346019/1‑100 (MQ=255)
                                               ggAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAg                                                    <  1:323740/100‑1 (MQ=255)
                                                gagaAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAgg                                                   <  2:1498317/100‑1 (MQ=255)
                                                 agaAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAggg                                                  <  2:1819083/100‑1 (MQ=255)
                                                  gaAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAgggn                                                 <  2:2426593/100‑2 (MQ=255)
                                                    aataatATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGgaga                                               <  2:3335545/100‑1 (MQ=255)
                                                    aataatATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGgaga                                               <  2:1296558/100‑1 (MQ=255)
                                                     ataataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGat                                              >  2:2147484/1‑100 (MQ=255)
                                                      taataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAATGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGata                                             <  2:2874654/100‑1 (MQ=255)
                                                      taataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGata                                             <  1:55060/100‑1 (MQ=255)
                                                          atGTGAAGTTCTTCGATGGTATGGAAAGCATTATATAAAAGGACCCAAGATTTATTGGGTTCTTTTTTGTCTATCACTACTATTAACattgctgatt                                            >  2:3354358/1‑88 (MQ=255)
                                                          atGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATc                                         >  1:3048497/1‑100 (MQ=255)
                                                             tgAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCAcc                                      <  1:1270835/100‑1 (MQ=255)
                                                             tgAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCAcc                                      >  2:2753280/1‑100 (MQ=255)
                                                             tgAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCAcc                                      <  1:2997690/100‑1 (MQ=255)
                                                              gAAGTTCTTCGATAGTCTGGAAGGCATTATATAAAAGGGCCGAGTAGTTATTGGGTTCTTTTTTCTCTATCAATCCTATTAGCAGGGAGATATATCAcc                                      >  2:2249457/1‑99 (MQ=255)
                                                              gAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCa                                     >  1:94870/1‑100 (MQ=255)
                                                                   tcttcGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGtt                                <  1:1126202/100‑1 (MQ=255)
                                                                   tcttcGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGtt                                >  1:3154762/1‑100 (MQ=255)
                                                                    cttcGATAGTATGTAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGtt                                <  2:2491549/99‑1 (MQ=255)
                                                                    cttcGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGttt                               >  1:1328588/1‑100 (MQ=255)
                                                                      tcGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTaa                             >  2:3365833/1‑100 (MQ=255)
                                                                       cGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATTTATATTATTAGCAGGGAGATATATCACCATTGTTTact                            >  2:36514/1‑98 (MQ=255)
                                                                       cGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAt                            >  2:1434977/1‑100 (MQ=255)
                                                                       cGATAGTATGGAAGGCATTATATAAAAGAACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAt                            <  2:3089786/100‑1 (MQ=255)
                                                                          taatATGGAAGGCATTATATAAAAAGGCCCAATATTTATTGGGGTCTTTTTTCCCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtg                         <  2:2935461/97‑1 (MQ=255)
                                                                          tAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtg                         >  1:788108/1‑100 (MQ=255)
                                                                          tAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtg                         <  1:3365858/100‑1 (MQ=255)
                                                                            gTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtgtg                       <  2:1743261/100‑1 (MQ=255)
                                                                               tGGAAGCCATTATATAAAAGGACCCAATATTTATTGTGTTCTTTTTTCTCTATCAATACTATTAGCAGGTAGATATATCACCAGAGTTTAAtgtgtg                       <  2:186865/97‑1 (MQ=255)
                                                                                   aGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAtttttt                <  1:1374206/100‑1 (MQ=255)
                                                                                     gCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttat              <  1:802232/100‑1 (MQ=255)
                                                                                     gCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTttta               <  1:1850650/99‑1 (MQ=255)
                                                                                      cATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttatt             <  1:3454966/100‑1 (MQ=255)
                                                                                      cATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttatt             <  1:1523179/100‑1 (MQ=255)
                                                                                       aTTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttattt            >  2:1264450/1‑100 (MQ=255)
                                                                                       aTTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATAAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttattt            >  2:614729/1‑100 (MQ=255)
                                                                                         tatataAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTttatttat          <  1:1195782/100‑1 (MQ=255)
                                                                                         tatataAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTttatttat          >  2:3215297/1‑100 (MQ=255)
                                                                                         agaggaaaaaGAACCCAATAAATATTGGGTCCTTTTGTCTCTATAAATACTATTCGCGGGGAGATATTTTACCAGAGTTTGATGTGTGATTTtttattt            >  2:2672849/7‑99 (MQ=16)
                                                                                         agagaaaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttattta           >  2:2816601/7‑100 (MQ=25)
                                                                                          gagaaaaaaGAACCCAATAATTATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATTACCAGAGTTTAATGTGTGATTTTttatttat          >  2:2227183/6‑100 (MQ=37)
                                                                                          gagaaaaaaGAACCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTttatttat          >  1:1502191/6‑100 (MQ=38)
                                                                                          gagaaaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttattta           >  2:2673367/6‑99 (MQ=25)
                                                                                          gagaaaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTttatttat          >  1:2785490/6‑100 (MQ=25)
                                                                                          gagaaaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTttatttat          >  1:1024399/6‑100 (MQ=25)
                                                                                          gagaaaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTttatttat          >  1:2084214/6‑100 (MQ=37)
                                                                                          gagaaaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAAATATATCACCAGAGTTTAATGTGTGATTTTgtatttat          >  2:1888443/6‑100 (MQ=25)
                                                                                           tataAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtc         >  2:3558662/1‑99 (MQ=255)
                                                                                            ataAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgt       <  1:591224/100‑1 (MQ=255)
                                                                                             taAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATAGCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgtc      >  1:1689227/1‑100 (MQ=255)
                                                                                                 aaGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGaa   >  1:3033125/1‑100 (MQ=255)
                                                                                                  aGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAAc  <  1:1428805/100‑1 (MQ=255)
                                                                                                    |                                                                                                 
CTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAAC  >  NC_000913/4261250‑4261446

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: