Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,262,255 G→C 17.7% S196T (AGC→ACC)  dusA → tRNA‑dihydrouridine synthase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,262,2550GC17.7% 537.6 / 105.1 216S196T (AGC→ACC) dusAtRNA‑dihydrouridine synthase A
Reads supporting (aligned to +/- strand):  ref base G (95/82);  new base C (15/23);  total (110/106)
Fisher's exact test for biased strand distribution p-value = 1.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01

Only 100 of 224 total aligned reads displayed.

GCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTAT  >  NC_000913/4262156‑4262354
                                                                                                   |                                                                                                   
gCTATGAATTTCTCTGCGATTTCATCAATCCCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAg                                                                                                     <  1:1526237/100‑1 (MQ=255)
gCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAg                                                                                                     >  1:3109365/1‑100 (MQ=255)
   aTGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGccc                                                                                                  >  1:1121425/1‑100 (MQ=255)
      aaTTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaa                                                                                               <  1:2918207/100‑1 (MQ=255)
       aTTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCgaaa                                                                                              <  1:2076492/100‑1 (MQ=255)
       aTTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCgaaa                                                                                              >  1:4233182/1‑100 (MQ=255)
        tttCTCTGCGATTTCATCAATACCGTGTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaag                                                                                             <  1:3811077/100‑1 (MQ=255)
         ttCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGGGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaagc                                                                                            <  1:2261137/100‑2 (MQ=255)
         ttCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCgaaaga                                                                                            <  1:2840165/100‑1 (MQ=255)
         ttCTCTGCGATTTCATCAATACCGTTGCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaagc                                                                                            <  1:1809805/100‑2 (MQ=255)
         ttCTCTGCGATTTCAGCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaagc                                                                                            <  1:3027284/100‑2 (MQ=255)
           ctctGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaagcca                                                                                          <  1:1864475/100‑4 (MQ=255)
           ctctGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaagcca                                                                                          <  2:1960501/100‑4 (MQ=255)
           ctctGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaagcca                                                                                          <  2:2706389/100‑4 (MQ=255)
            tctGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaagccag                                                                                         <  1:336529/100‑5 (MQ=39)
            tctGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaagccag                                                                                         <  1:2142500/100‑5 (MQ=39)
            tctGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaagccag                                                                                         <  2:3104383/100‑5 (MQ=39)
             ctGCGATTTCATCAATACCGTTTCCGGCAAAGGTGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAAc                                                                                        >  2:4009067/1‑100 (MQ=255)
             ctGCGATTTCATCAATAACGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAAc                                                                                        <  1:2187218/100‑1 (MQ=255)
              tGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAAcc                                                                                       >  1:948489/1‑100 (MQ=255)
               gCGATTTCATCAATACCGTTTACGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCCCGTAAAGCCGGGCTTTCGGGGTTAAGCCCGAAAGAAAACCg                                                                                      >  1:3468279/1‑100 (MQ=255)
                 gATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTg                                                                                    <  2:265866/100‑1 (MQ=255)
                 gATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTg                                                                                    <  1:1173501/100‑1 (MQ=255)
                  aTTTCATCAATACCGTTTCCGTCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGTGTTAAGCCCGAAAGAAAACCGTGa                                                                                   >  1:4267525/1‑100 (MQ=255)
                     tcatcaATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAAt                                                                                >  2:552402/1‑100 (MQ=255)
                      catcaATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATc                                                                               >  2:153746/1‑100 (MQ=255)
                      catcaATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATc                                                                               >  2:621350/1‑100 (MQ=255)
                      catcaATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATc                                                                               >  1:1012689/1‑100 (MQ=255)
                       atcaATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATcc                                                                              >  2:3595547/1‑100 (MQ=255)
                        tcaATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATc                                                                               >  2:2700717/1‑98 (MQ=255)
                         caATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTTATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATccc                                                                             >  1:2585738/1‑99 (MQ=255)
                          aaTACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCccgc                                                                           >  1:3257690/1‑100 (MQ=255)
                           aTACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCccgcc                                                                          <  2:3896312/100‑1 (MQ=255)
                             aCCGTTTCCGGCAAAGGCGAGTGTGAGTTGTTCATCATCGACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCccg                                                                            <  1:1885716/96‑1 (MQ=255)
                             aCCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCcgccgc                                                                        >  2:1114113/1‑100 (MQ=255)
                              ccGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCt                                                                       >  1:835851/1‑100 (MQ=255)
                               cGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTc                                                                      <  2:4338869/100‑1 (MQ=255)
                                gTTTCCGGCAAAGGCGAGTTTGAGATGTTCATTATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAATTCCGTGAAATCCCGCCGCTc                                                                      <  1:1252573/99‑1 (MQ=255)
                                gTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCg                                                                     <  2:183507/100‑1 (MQ=255)
                                    ccGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTa                                                                 >  2:569865/1‑100 (MQ=255)
                                     cGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTAt                                                                >  2:2889443/1‑100 (MQ=255)
                                     cGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCCCg                                                                     >  1:3385752/1‑95 (MQ=255)
                                      ggCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATc                                                               >  1:4012545/1‑100 (MQ=255)
                                      ggCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATc                                                               >  2:1046179/1‑100 (MQ=255)
                                       gCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATcc                                                              >  2:241996/1‑100 (MQ=255)
                                       gCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATcc                                                              >  2:3812260/1‑100 (MQ=255)
                                       gCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATcc                                                              <  1:2448570/100‑1 (MQ=255)
                                       gCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATcc                                                              >  1:3952608/1‑100 (MQ=255)
                                        cAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCcg                                                             <  2:2610406/100‑1 (MQ=255)
                                         aaaGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCcgc                                                            <  1:135385/100‑1 (MQ=255)
                                          aaGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCcgcg                                                           <  2:2467024/100‑1 (MQ=255)
                                            ggCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCgtg                                                         >  2:3451158/1‑100 (MQ=255)
                                              cGAGTGTGAGATGTTCATCAGCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCgtgtg                                                       <  2:3305228/100‑1 (MQ=255)
                                                aGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTa                                                     <  2:3095098/100‑1 (MQ=255)
                                                   gtgAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCa                                                  >  1:3569251/1‑100 (MQ=255)
                                                    tgAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCaa                                                 >  1:4348604/1‑100 (MQ=255)
                                                     gggaTGTTCAGCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAAc                                                <  2:78861/98‑1 (MQ=255)
                                                     gagaTGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGCTTATCCGCGTGTGTATCAAc                                                >  1:2536751/1‑100 (MQ=255)
                                                     gagaTGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAAc                                                >  2:4143722/1‑100 (MQ=255)
                                                     gagaTGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAAc                                                >  2:3138357/1‑100 (MQ=255)
                                                      agaTGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACt                                               >  2:1124409/1‑100 (MQ=255)
                                                      agaTGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACt                                               <  2:2672843/100‑1 (MQ=255)
                                                       gaTGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACt                                               <  2:3499734/99‑1 (MQ=255)
                                                       gaTGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTg                                              >  1:3984472/1‑100 (MQ=255)
                                                            tcatcaTCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAg                                           <  2:1667316/98‑1 (MQ=255)
                                                             catcatCCACGCACGTAAAGCCTGGCTTTCGGAGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGt                                        >  1:3880133/1‑100 (MQ=255)
                                                              atcatcCACGCACGTAAAGCGTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGCAGCg                                         >  1:2691296/1‑98 (MQ=255)
                                                               tcatcCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCTAAAGAAAACCGTGAAATCCCGCCGCTCAATTATCCGCGTGTGTATCAACTGAAGCGTGa                                      <  1:817788/100‑1 (MQ=255)
                                                                catcCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGAc                                     <  2:1812591/100‑1 (MQ=255)
                                                                catcCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGAc                                     >  1:3715355/1‑100 (MQ=255)
                                                                    cacgcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACttt                                  >  2:2550260/1‑99 (MQ=255)
                                                                    cacgcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTc                                 <  2:3923786/100‑1 (MQ=255)
                                                                     acgcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTcc                                >  1:1738125/1‑100 (MQ=255)
                                                                     acgcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTcc                                <  2:1478152/100‑1 (MQ=255)
                                                                            tcacggttttcTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTg                         >  2:1421362/11‑100 (MQ=37)
                                                                            tAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAATCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTTAAGCGTGACTGTCCGCATCTg                         <  1:45775/100‑1 (MQ=255)
                                                                            tAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTg                         <  2:4348824/100‑1 (MQ=255)
                                                                               cggttttcTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTGTCCGCGAGTGTATCAACTGAAGCGTGACATTCCGCATCTggc                       >  1:2815172/8‑97 (MQ=37)
                                                                               cggttttcTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACa                      >  2:336274/8‑100 (MQ=38)
                                                                               aGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACa                      <  2:3661671/100‑1 (MQ=255)
                                                                                 gttttcTTTCGGGCTTAACCCCGGAAGAAAACCGTGAAATTCCGCCGCTCCATTATCCGCGAATGTATCAACTCAAGCGTGACTTTCCGCCTCTGa                        >  1:3046019/6‑96 (MQ=255)
                                                                                  cTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATg                   >  2:2200497/1‑100 (MQ=255)
                                                                                    ggCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTc                 >  2:2839585/1‑100 (MQ=255)
                                                                                    ggCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTc                 <  2:298404/100‑1 (MQ=255)
                                                                                      cTTTCGGGGTTAAGAGCGAAAGCAAAGCGTGAGATCCCGGCGGTTGAATATGCGGGTGTGTATCAACTGAAGAGCGACTGTCTGCAGCTGaacatctc                 >  1:3638287/1‑91 (MQ=255)
                                                                                      cTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGa               >  2:1203132/1‑100 (MQ=255)
                                                                                      cTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGa               >  2:123891/1‑100 (MQ=255)
                                                                                      cTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGa               >  2:4223970/1‑100 (MQ=255)
                                                                                      cTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATACGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGa               >  1:776542/1‑100 (MQ=255)
                                                                                       tttCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGAt              <  1:1421305/100‑1 (MQ=255)
                                                                                        ttCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGAtt             <  2:3569270/100‑1 (MQ=255)
                                                                                           ggggTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAAc          <  2:2478047/100‑1 (MQ=255)
                                                                                            ggtttAAGCCCGCAAGAAAACCGTGAAATCCCGCCTCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGTCTTTTCGCATCTGACAATGTCGATTAACg         <  1:1216981/97‑1 (MQ=255)
                                                                                             ggTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACgg        >  2:3606812/1‑100 (MQ=255)
                                                                                             ggTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACgg        <  1:689651/100‑1 (MQ=255)
                                                                                              gTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACggt       <  1:2398942/100‑1 (MQ=255)
                                                                                                tAAGCCCTAAAAAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACgg        >  1:2299480/1‑97 (MQ=255)
                                                                                                tAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGCTTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGCCAATGTCGATTAACggtgg     >  1:3459050/1‑100 (MQ=255)
                                                                                                tAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACggtgg     >  2:467618/1‑100 (MQ=255)
                                                                                                   gCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTAt  >  1:255288/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
GCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTAT  >  NC_000913/4262156‑4262354

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