Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,261,359 T→A 17.5% intergenic (+53/‑310) yjbM → / → dusA uncharacterized protein/tRNA‑dihydrouridine synthase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,261,3590TA17.5% 418.4 / 77.1 172intergenic (+53/‑310)yjbM/dusAuncharacterized protein/tRNA‑dihydrouridine synthase A
Reads supporting (aligned to +/- strand):  ref base T (67/75);  new base A (13/17);  total (80/92)
Fisher's exact test for biased strand distribution p-value = 8.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.49e-01

Only 100 of 184 total aligned reads displayed.

CTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTA  >  NC_000913/4261260‑4261458
                                                                                                    |                                                                                                   
cTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAAtatt                                                                                                     >  1:2943938/1‑100 (MQ=255)
   cTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAAtatttat                                                                                                  >  1:2745496/1‑100 (MQ=255)
    tAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAAtatttatt                                                                                                 <  1:2213074/100‑1 (MQ=255)
     aaaTTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTg                                                                                                <  1:1291110/100‑1 (MQ=255)
     aaaTTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTg                                                                                                >  2:2589567/1‑100 (MQ=255)
     aaaTTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTg                                                                                                <  1:555841/100‑1 (MQ=255)
      aattaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTgg                                                                                               >  2:1803393/1‑100 (MQ=255)
         taattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATAAATATTGGGtc                                                                                            <  1:3376344/100‑2 (MQ=255)
          aattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTc                                                                                           >  1:2271817/1‑100 (MQ=255)
          aattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGAACCCAATAAATATTGGGtcc                                                                                           <  1:3607094/100‑3 (MQ=39)
          aattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGAACCCAATAAATATTGGGtcc                                                                                           <  1:79650/100‑3 (MQ=39)
           attTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCt                                                                                          >  1:1420514/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)
                 gATTTGTTCGGGGTTAAATCGGAGAAATAATATTTTAAGTTTTTCGATAGTATGGAAGGCATTATAT‑AAAATGACCCAATATTTATTGGGTTCTTTTTtc                                                                                    <  2:694334/100‑1 (MQ=255)
                 gATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTtc                                                                                    <  1:3635197/100‑1 (MQ=255)
                 gATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGAACCCAATAAATATTGGGTCCttttata                                                                                    <  2:3195260/100‑4 (MQ=37)
                  aattGTTCGGGGTTAAATCGGATAAATAATATAACAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGAACCCAATAAATATTGGGTCCttttatat                                                                                   <  2:3214605/98‑5 (MQ=21)
                  aTTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTtct                                                                                   >  2:2356571/1‑100 (MQ=255)
                  aTTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTtct                                                                                   <  1:2315119/100‑1 (MQ=255)
                   tttGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTtctc                                                                                  <  1:945652/100‑1 (MQ=255)
                   tttGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTtctc                                                                                  <  1:2364028/100‑1 (MQ=255)
                   tttGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTtctc                                                                                  <  1:1138949/100‑1 (MQ=255)
                     tGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTa                                                                                <  1:2087361/100‑1 (MQ=255)
                      gTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAAGAACCCAATAAATATTGGGTCCttttatagatc                                                                                >  1:490973/1‑93 (MQ=25)
                      gTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTAt                                                                               >  1:1502567/1‑100 (MQ=255)
                      gTTCGGGGTTAAATCGGAGAAATAATATGTGAAGGTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAAGATTAATTGGGTTCTTTTTTCTCTAt                                                                               >  2:887559/1‑100 (MQ=255)
                        tCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGAACCCAATAAATATTGGGTCCttttatataatgcc                                                                             <  2:3047272/100‑11 (MQ=21)
                          ggggTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAAt                                                                           >  1:1514257/1‑100 (MQ=255)
                           gggTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATa                                                                          <  2:1399641/100‑1 (MQ=255)
                            ggTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATTAATAc                                                                         >  2:1833864/1‑100 (MQ=255)
                            ggTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATAc                                                                         <  1:2549898/100‑1 (MQ=255)
                             gTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACt                                                                        >  2:1241016/1‑100 (MQ=255)
                             gTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACt                                                                        >  1:1999704/1‑100 (MQ=255)
                                 aaTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTa                                                                    >  1:2064349/1‑100 (MQ=255)
                                   tCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGc                                                                  <  1:2489600/100‑1 (MQ=255)
                                    cGGAGAAATAATATGTGAAGTTCTTCGATATTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCa                                                                 <  2:1271151/100‑1 (MQ=255)
                                    cGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCa                                                                 >  1:3098652/1‑100 (MQ=255)
                                    cGGAGAAATAATATGTGAAGCTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCa                                                                 >  2:3346019/1‑100 (MQ=255)
                                     ggAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAg                                                                <  1:323740/100‑1 (MQ=255)
                                     ggAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAg                                                                >  2:3454941/1‑100 (MQ=255)
                                     ggAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAg                                                                >  2:3568335/1‑100 (MQ=255)
                                      gagaAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAgg                                                               <  2:1498317/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)
                                                atGTGAAGTTCTTCGATGGTATGGAAAGCATTATAT‑AAAAGGACCCAAGATTTATTGGGTTCTTTTTTGTCTATCACTACTATTAACattgctgatt                                                        >  2:3354358/1‑88 (MQ=255)
                                                atGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATc                                                     >  1:1472058/1‑100 (MQ=255)
                                                atGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATc                                                     >  1:3048497/1‑100 (MQ=255)
                                                 tgtgAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCa                                                    <  1:530218/100‑1 (MQ=255)
                                                  gtgAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCAc                                                   <  1:1827416/100‑1 (MQ=255)
                                                   tgAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCAcc                                                  >  2:2753280/1‑100 (MQ=255)
                                                   tgAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCAcc                                                  <  1:2997690/100‑1 (MQ=255)
                                                   tgAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCAcc                                                  >  2:1750214/1‑100 (MQ=255)
                                                   tgAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCAc                                                   >  1:2852145/1‑99 (MQ=255)
                                                         tcttcGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGtt                                            >  1:3154762/1‑100 (MQ=255)
                                                          cttcGATAGTATGTAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGtt                                            <  2:2491549/99‑1 (MQ=255)
                                                            tcGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTaa                                         <  2:308333/100‑1 (MQ=255)
                                                            tcGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTa                                          <  2:2822208/99‑1 (MQ=255)
                                                             cGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATTTATATTATTAGCAGGGAGATATATCACCATTGTTTact                                        >  2:36514/1‑98 (MQ=255)
                                                             cGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAt                                        >  2:1434977/1‑100 (MQ=255)
                                                                taatATGGAAGGCATTATAT‑AAAAAGGCCCAATATTTATTGGGGTCTTTTTTCCCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtg                                     <  2:2935461/97‑1 (MQ=255)
                                                                tAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtg                                     <  1:3365858/100‑1 (MQ=255)
                                                                  gTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtgtg                                   >  1:3390678/1‑100 (MQ=255)
                                                                  gTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtgtg                                   <  2:1743261/100‑1 (MQ=255)
                                                                     tGGAAGGCATTATAT‑AAAAGGACCTAATATTTATTGGGTTCTTTTTTCTCTGTTAATAATATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAt                                 <  2:3308765/99‑1 (MQ=255)
                                                                        aaGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAttttt                             >  1:701001/1‑100 (MQ=255)
                                                                         aGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAtttttt                            <  1:1374206/100‑1 (MQ=255)
                                                                          ggCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTttta                           <  1:2033209/100‑1 (MQ=255)
                                                                           gCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttat                          <  2:3619981/100‑1 (MQ=255)
                                                                           gCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttat                          <  1:802232/100‑1 (MQ=255)
                                                                            cATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttatt                         <  1:1523179/100‑1 (MQ=255)
                                                                            cATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttatt                         >  1:2591701/1‑100 (MQ=255)
                                                                             aTTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATAAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttattt                        >  2:614729/1‑100 (MQ=255)
                                                                              ttATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttattta                       <  1:1750222/100‑1 (MQ=255)
                                                                               agaggaaaaaGAACCCAATAAATATTGGGTCCTTTTGTCTCTATAAATACTATTCGCGGGGAGATATTTTACCAGAGTTTGATGTGTGATTTtttattt                        >  2:2672849/7‑99 (MQ=16)
                                                                                gagaaaaaaGAACCCAATAATTATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATTACCAGAGTTTAATGTGTGATTTTttatttat                      >  2:2227183/6‑100 (MQ=37)
                                                                                gagaaaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttattta                       >  2:2673367/6‑99 (MQ=25)
                                                                                gagaaaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTttatttat                      >  1:1024399/6‑100 (MQ=25)
                                                                                gagaaaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTttatttat                      >  1:126548/6‑100 (MQ=25)
                                                                                gagaaaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTttatttat                      >  1:2785490/6‑100 (MQ=25)
                                                                                gagaaaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAAATATATCACCAGAGTTTAATGTGTGATTTTgtatttat                      >  2:1888443/6‑100 (MQ=25)
                                                                                   t‑aAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATAGCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgtc                  >  1:1689227/1‑100 (MQ=255)
                                                                                       aaGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGaa               >  1:3033125/1‑100 (MQ=255)
                                                                                          aaCCCAATAAATATTGGGTTCTTTTTTCTCTATTTATACTATTAGCAGGGAGATATATCACCATATTTTAATGTGTGATTTTTTATTTATAGTCTaa               >  2:311774/2‑97 (MQ=255)
                                                                                            cccAATATTTATTGGGTTCTTTTTTTTTTTTTTTTAATATTAGAAAGGAGATATATCACCATATTTTTATGTGTGATTTTTTTTTTAtctt                   >  2:1506924/1‑89 (MQ=255)
                                                                                            cccAATATTTATTGGGTTCTTTTTTTTCTTTCTATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCt            >  2:1058635/1‑98 (MQ=255)
                                                                                            cccAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTgg          >  1:2487628/1‑100 (MQ=255)
                                                                                            cccAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTg           <  2:1230862/99‑1 (MQ=255)
                                                                                            cccAATACTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTgg          >  1:1911668/1‑100 (MQ=255)
                                                                                              cAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGAt        <  2:448198/100‑1 (MQ=255)
                                                                                              cAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGAt        >  1:762679/1‑100 (MQ=255)
                                                                                                atatTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTg      >  1:101366/1‑100 (MQ=255)
                                                                                                   tttattGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGttt   <  2:224514/100‑1 (MQ=255)
                                                                                                    ttattGGGTTCTTTTTTCTCTATCAATACTATTTGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTa  <  2:278255/100‑1 (MQ=255)
                                                                                                    |                                                                                                   
CTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTA  >  NC_000913/4261260‑4261458

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