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.

CTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAAC  >  NC_000913/4261250‑4261446
                                                                                                    |                                                                                                 
cTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAgg                                                                                                   >  1:448200/1‑100 (MQ=255)
  tACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGAc                                                                                                 >  1:2857859/1‑100 (MQ=255)
     aGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCa                                                                                              <  1:3632054/100‑1 (MQ=255)
      gACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCaa                                                                                             >  1:3199479/1‑100 (MQ=255)
      gACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCaa                                                                                             >  1:2928590/1‑100 (MQ=255)
          cTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAAtatt                                                                                         >  1:2943938/1‑100 (MQ=255)
          cTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAAtatt                                                                                         >  2:2830786/1‑100 (MQ=255)
          cTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAAtatt                                                                                         >  1:2027801/1‑100 (MQ=255)
            tCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAAtattta                                                                                       <  1:3368214/100‑1 (MQ=255)
             cTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAAtatttat                                                                                      >  1:2745496/1‑100 (MQ=255)
              tAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAAtatttatt                                                                                     <  2:1162168/100‑1 (MQ=255)
              tAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAAtatttatt                                                                                     <  1:2213074/100‑1 (MQ=255)
               aaaTTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTg                                                                                    >  2:2589567/1‑100 (MQ=255)
               aaaTTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTg                                                                                    <  1:3378864/100‑1 (MQ=255)
                aattaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTgg                                                                                   >  2:3202285/1‑100 (MQ=255)
                aattaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTgg                                                                                   <  1:2505251/100‑1 (MQ=255)
                   taaatTATTATTTGTTCGGGCTTAAATCGGAGAAATAATATGTGAAGTTATACGATAATATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGt                                                                                 <  2:1774321/99‑1 (MQ=255)
                    aattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTc                                                                               >  1:2271817/1‑100 (MQ=255)
                    aattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATAAATATTGGGtc                                                                                <  1:1800242/99‑2 (MQ=255)
                     attTATTATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGGAGGCATTATAT‑AAAAGAACCCAATAAATATTGGGTCCt                                                                              <  2:3593181/100‑1 (MQ=38)
                     attTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCt                                                                              >  1:1506636/1‑100 (MQ=255)
                      tttATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGAACCCAATAAATATTGGGTCCtt                                                                             <  1:2073321/100‑1 (MQ=38)
                       ttATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCttt                                                                            >  2:559479/1‑100 (MQ=255)
                       ttATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGAACCCAATAAATATTGGGTCCttt                                                                            <  1:2702355/100‑1 (MQ=38)
                        tATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCtttt                                                                           <  1:1343867/100‑1 (MQ=255)
                         aTGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGAACCCAATAAATATTGGGTCCtttta                                                                          <  2:2042417/100‑2 (MQ=38)
                           gatctGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAAGAACCCAATAAATATTGGGTCCttttat                                                                         <  2:490971/96‑3 (MQ=25)
                           gATTTGTTCGGGGTTAAATCGGAGAAATAATATTTTAAGTTTTTCGATAGTATGGAAGGCATTATAT‑AAAATGACCCAATATTTATTGGGTTCTTTTTtc                                                                        <  2:694334/100‑1 (MQ=255)
                           gATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTtc                                                                        <  1:2552133/100‑1 (MQ=255)
                           gATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTtc                                                                        <  2:3069095/100‑1 (MQ=255)
                           gATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGAACCCAATAAATATTGGGTCCttttata                                                                        <  1:1948462/100‑4 (MQ=37)
                           gATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGAACCCAATAAATATTGGGTCCttttata                                                                        <  2:3195260/100‑4 (MQ=37)
                            aTTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTtct                                                                       <  1:2315119/100‑1 (MQ=255)
                            aTTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTtct                                                                       >  2:2356571/1‑100 (MQ=255)
                             tttGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTtctc                                                                      <  1:2362039/100‑1 (MQ=255)
                              ttGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTctct                                                                     <  1:428028/100‑1 (MQ=255)
                              ttGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTctct                                                                     <  1:710422/100‑1 (MQ=255)
                               tGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTa                                                                    <  1:2087361/100‑1 (MQ=255)
                                gTTCGGGGTTAAATCGGAGAAATAATATGTGAAGGTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAAGATTAATTGGGTTCTTTTTTCTCTAt                                                                   >  2:887559/1‑100 (MQ=255)
                                 ttCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATc                                                                  <  1:2061491/100‑1 (MQ=255)
                                    ggggTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAAt                                                               >  1:2416142/1‑100 (MQ=255)
                                    ggggTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAAt                                                               >  1:1514257/1‑100 (MQ=255)
                                     gggTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATa                                                              <  2:1399641/100‑1 (MQ=255)
                                       gTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACt                                                            >  2:1241016/1‑100 (MQ=255)
                                           aaTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTa                                                        >  1:2701840/1‑100 (MQ=255)
                                              cGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCa                                                     >  1:3098652/1‑100 (MQ=255)
                                              cGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCa                                                     >  1:500144/1‑100 (MQ=255)
                                              cGGAGAAATAATATGTGAAGCTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCa                                                     >  2:3346019/1‑100 (MQ=255)
                                               ggAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAg                                                    <  1:323740/100‑1 (MQ=255)
                                                gagaAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAgg                                                   <  2:1498317/100‑1 (MQ=255)
                                                 agaAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAggg                                                  <  2:1819083/100‑1 (MQ=255)
                                                  gaAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAgggn                                                 <  2:2426593/100‑2 (MQ=255)
                                                    aataatATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGgaga                                               <  2:3335545/100‑1 (MQ=255)
                                                    aataatATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGgaga                                               <  2:1296558/100‑1 (MQ=255)
                                                     ataataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGat                                              >  2:2147484/1‑100 (MQ=255)
                                                      taataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAATGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGata                                             <  2:2874654/100‑1 (MQ=255)
                                                      taataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGata                                             <  1:55060/100‑1 (MQ=255)
                                                          atGTGAAGTTCTTCGATGGTATGGAAAGCATTATAT‑AAAAGGACCCAAGATTTATTGGGTTCTTTTTTGTCTATCACTACTATTAACattgctgatt                                            >  2:3354358/1‑88 (MQ=255)
                                                          atGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATc                                         >  1:3048497/1‑100 (MQ=255)
                                                             tgAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCAcc                                      <  1:1270835/100‑1 (MQ=255)
                                                             tgAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCAcc                                      >  2:2753280/1‑100 (MQ=255)
                                                             tgAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCAcc                                      <  1:2997690/100‑1 (MQ=255)
                                                              gAAGTTCTTCGATAGTCTGGAAGGCATTATAT‑AAAAGGGCCGAGTAGTTATTGGGTTCTTTTTTCTCTATCAATCCTATTAGCAGGGAGATATATCAcc                                      >  2:2249457/1‑99 (MQ=255)
                                                              gAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCa                                     >  1:94870/1‑100 (MQ=255)
                                                                   tcttcGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGtt                                <  1:1126202/100‑1 (MQ=255)
                                                                   tcttcGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGtt                                >  1:3154762/1‑100 (MQ=255)
                                                                    cttcGATAGTATGTAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGtt                                <  2:2491549/99‑1 (MQ=255)
                                                                    cttcGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGttt                               >  1:1328588/1‑100 (MQ=255)
                                                                      tcGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTaa                             >  2:3365833/1‑100 (MQ=255)
                                                                       cGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATTTATATTATTAGCAGGGAGATATATCACCATTGTTTact                            >  2:36514/1‑98 (MQ=255)
                                                                       cGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAt                            >  2:1434977/1‑100 (MQ=255)
                                                                       cGATAGTATGGAAGGCATTATAT‑AAAAGAACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAt                            <  2:3089786/100‑1 (MQ=255)
                                                                          taatATGGAAGGCATTATAT‑AAAAAGGCCCAATATTTATTGGGGTCTTTTTTCCCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtg                         <  2:2935461/97‑1 (MQ=255)
                                                                          tAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtg                         >  1:788108/1‑100 (MQ=255)
                                                                          tAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtg                         <  1:3365858/100‑1 (MQ=255)
                                                                            gTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtgtg                       <  2:1743261/100‑1 (MQ=255)
                                                                               tGGAAGCCATTATAT‑AAAAGGACCCAATATTTATTGTGTTCTTTTTTCTCTATCAATACTATTAGCAGGTAGATATATCACCAGAGTTTAAtgtgtg                       <  2:186865/97‑1 (MQ=255)
                                                                                   aGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAtttttt                <  1:1374206/100‑1 (MQ=255)
                                                                                     gCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttat              <  1:802232/100‑1 (MQ=255)
                                                                                     gCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTttta               <  1:1850650/99‑1 (MQ=255)
                                                                                      cATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttatt             <  1:3454966/100‑1 (MQ=255)
                                                                                      cATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttatt             <  1:1523179/100‑1 (MQ=255)
                                                                                       aTTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttattt            >  2:1264450/1‑100 (MQ=255)
                                                                                       aTTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATAAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttattt            >  2:614729/1‑100 (MQ=255)
                                                                                         tatat‑aAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTttatttat          <  1:1195782/100‑1 (MQ=255)
                                                                                         tatat‑aAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTttatttat          >  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)
                                                                                           tat‑aAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtc         >  2:3558662/1‑99 (MQ=255)
                                                                                            at‑aAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgt       <  1:591224/100‑1 (MQ=255)
                                                                                             t‑aAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATAGCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgtc      >  1:1689227/1‑100 (MQ=255)
                                                                                                 aaGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGaa   >  1:3033125/1‑100 (MQ=255)
                                                                                                  aGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAAc  <  1:1428805/100‑1 (MQ=255)
                                                                                                    |                                                                                                 
CTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAAC  >  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: ‑