Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,262,255 G→C 18.9% 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,2550GC18.9% 470.2 / 105.2 191S196T (AGC→ACC) dusAtRNA‑dihydrouridine synthase A
Reads supporting (aligned to +/- strand):  ref base G (76/79);  new base C (21/15);  total (97/94)
Fisher's exact test for biased strand distribution p-value = 3.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.39e-01

Only 100 of 197 total aligned reads displayed.

GCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTAT  >  NC_000913/4262156‑4262354
                                                                                                   |                                                                                                   
gCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAg                                                                                                     <  1:903796/100‑1 (MQ=255)
gCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGGGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAg                                                                                                     >  2:710010/1‑100 (MQ=255)
  tATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGcc                                                                                                   <  2:15590/100‑1 (MQ=255)
  tATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGcc                                                                                                   <  1:2503181/100‑1 (MQ=255)
      aaTTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGccc                                                                                                  >  2:1957120/1‑97 (MQ=255)
       aTTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCgaa                                                                                               <  1:932157/99‑1 (MQ=255)
        tttCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCgaaag                                                                                             <  1:1409549/100‑1 (MQ=255)
        tttCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCgaaa                                                                                              <  1:3503306/99‑1 (MQ=255)
        tttCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaag                                                                                             <  1:3710874/100‑1 (MQ=255)
           ctctGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaagcca                                                                                          <  2:3050657/100‑4 (MQ=255)
           ctctGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACACACGTAAAGCCTGGCTTTCGGGGTTAAGCCCgaaagaaa                                                                                          >  1:3324726/1‑100 (MQ=255)
            tctGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaagccag                                                                                         <  1:1145946/100‑5 (MQ=39)
            tctGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaagccag                                                                                         <  2:2857003/100‑5 (MQ=39)
             ctGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAAc                                                                                        <  1:3586426/100‑1 (MQ=255)
              tGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAAcc                                                                                       >  1:1037810/1‑100 (MQ=255)
              tGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaagccaggc                                                                                       <  1:108293/100‑7 (MQ=38)
              tGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaagccaggc                                                                                       <  1:671994/100‑7 (MQ=38)
               gCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaagccaggct                                                                                      <  1:1458221/100‑8 (MQ=38)
               gCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaagccaggct                                                                                      <  2:1851877/100‑8 (MQ=38)
                   tttCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAACCCGTGa                                                                                   <  2:1389983/99‑1 (MQ=255)
                   tttCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGag                                                                                  <  2:58798/100‑2 (MQ=255)
                      catcaATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGACATc                                                                               >  2:1692382/1‑100 (MQ=255)
                       agcaATAAGGTTTCCGGCAAAGTCGAGTGTGAGATGTTCATCAACCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGaaa                                                                                 <  2:2579880/95‑1 (MQ=255)
                         caATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCccg                                                                            >  1:990752/1‑100 (MQ=255)
                         caATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCccg                                                                            <  1:1780870/100‑1 (MQ=255)
                          aaTACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCccgc                                                                           <  1:1943527/100‑1 (MQ=255)
                              ccGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCt                                                                       >  1:1669318/1‑100 (MQ=255)
                              ccGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCt                                                                       >  1:3808425/1‑100 (MQ=255)
                              ccGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCt                                                                       >  1:299218/1‑100 (MQ=255)
                              ccGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCt                                                                       >  1:1811402/1‑100 (MQ=255)
                               cGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTc                                                                      <  2:2496834/100‑1 (MQ=255)
                               cGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATACACCCTCGTTCAGCCTGGCTTTCGGGGTACAGCCCGAAAGAAAACAGTGAAAACacgcag                                                                         >  2:3371379/1‑95 (MQ=255)
                                gTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCg                                                                     >  1:1338395/1‑100 (MQ=255)
                                gTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCg                                                                     >  1:3241292/1‑100 (MQ=255)
                                 ttgccGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGa                                                                    <  2:2692021/97‑1 (MQ=255)
                                  ttCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGAt                                                                   <  1:2749688/100‑1 (MQ=255)
                                  ttCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGAt                                                                   <  1:1585868/100‑1 (MQ=255)
                                     cGGCAAAGGCGCGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTAt                                                                >  2:762240/1‑100 (MQ=255)
                                     cGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGACAGAAAACCGTGAAATCCCGCCGATCGCTTAt                                                                >  2:847764/1‑100 (MQ=255)
                                     cGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTAt                                                                >  1:1593396/1‑100 (MQ=255)
                                      ggCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATc                                                               >  2:3738351/1‑100 (MQ=255)
                                          aaGGCGGGTGTGAGAGGTTCATCATCCACGCATGTAAATCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCcgg                                                            <  2:3510314/99‑2 (MQ=255)
                                          aaGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCcgcg                                                           >  2:964229/1‑100 (MQ=255)
                                            ggCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCgtg                                                         <  1:975385/100‑1 (MQ=255)
                                              cGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGGTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCgtgtg                                                       <  1:341622/100‑1 (MQ=255)
                                                   gtgAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCa                                                  >  1:2962480/1‑100 (MQ=255)
                                                     gagaTGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCgtgggagtgaac                                                >  2:3005693/1‑91 (MQ=255)
                                                        aTGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGa                                             <  2:2445261/100‑1 (MQ=255)
                                                         tGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGaa                                            >  2:1538542/1‑100 (MQ=255)
                                                         tGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGaa                                            <  1:2913080/100‑1 (MQ=255)
                                                          gTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAg                                           <  1:2731261/100‑1 (MQ=255)
                                                            tcatcaTCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCg                                         >  2:3756105/1‑100 (MQ=255)
                                                              atcatcCACGCACGTTAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTg                                       <  1:3758651/100‑1 (MQ=255)
                                                              atcatcCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTg                                       <  1:450520/100‑1 (MQ=255)
                                                               tcatcCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGa                                      <  1:102672/100‑1 (MQ=255)
                                                                catcCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGAc                                     >  2:2403511/1‑100 (MQ=255)
                                                                 atcCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGGGTGTATCAACTGAAGCGTGACt                                    >  1:1894521/1‑100 (MQ=255)
                                                                    cacgcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTc                                 >  2:1279346/1‑100 (MQ=255)
                                                                    cacgcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTc                                 >  1:2766388/1‑100 (MQ=255)
                                                                     acgcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTcc                                <  2:2537763/100‑1 (MQ=255)
                                                                      cgcacgTAAAGCCTCGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCg                               <  2:2606246/100‑1 (MQ=255)
                                                                        cacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCa                             >  2:329795/1‑100 (MQ=255)
                                                                        cacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCAATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCa                             <  2:2021136/100‑1 (MQ=255)
                                                                           gTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCt                          >  2:463903/1‑100 (MQ=255)
                                                                            tcacggttttcTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTg                         >  2:3106946/11‑100 (MQ=37)
                                                                               cggttttcTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACa                      >  1:2803924/8‑100 (MQ=38)
                                                                                gCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGCTTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACaa                     <  2:972462/100‑1 (MQ=255)
                                                                                  cTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATg                   <  1:3629188/100‑1 (MQ=255)
                                                                                  cTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATg                   >  1:701970/1‑100 (MQ=255)
                                                                                   tGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGt                  >  1:2219743/1‑100 (MQ=255)
                                                                                    ttcTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTc                 >  1:55527/3‑100 (MQ=255)
                                                                                    ttcTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTc                 >  1:789323/3‑100 (MQ=255)
                                                                                    ttcTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTc                 >  1:1502842/3‑100 (MQ=255)
                                                                                    ggCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTc                 >  2:142760/1‑100 (MQ=255)
                                                                                     tcTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCg                >  2:1778014/2‑100 (MQ=255)
                                                                                     tcTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCg                >  1:1273216/2‑100 (MQ=255)
                                                                                     tcTTTCGGGCTAAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCg                >  2:1994267/2‑100 (MQ=255)
                                                                                     gCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTcc                <  2:2228219/100‑2 (MQ=255)
                                                                                      cTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCTCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGAGTTTCCGCATCTGACAATGGCGa               >  2:2275301/1‑100 (MQ=255)
                                                                                      cTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGa               <  2:2160669/100‑1 (MQ=255)
                                                                                      cTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGa               >  2:2745694/1‑100 (MQ=255)
                                                                                      cTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGa               >  1:2105320/1‑100 (MQ=255)
                                                                                      cTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGa               >  1:2650440/1‑100 (MQ=255)
                                                                                      cTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGa               >  2:3639718/1‑100 (MQ=255)
                                                                                       tttCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGGGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGAt              >  2:517080/1‑100 (MQ=255)
                                                                                       tttCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGAt              >  1:1736599/1‑100 (MQ=255)
                                                                                        ttCGGGGCTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCGACTGAAGCGTGACTTTCCGCATCTGACAATGTCGAtn             <  2:2046250/100‑2 (MQ=255)
                                                                                         tCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTa            <  2:1037801/100‑1 (MQ=255)
                                                                                          cGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTaa           <  2:1240288/100‑1 (MQ=255)
                                                                                           ggggTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAAc          >  2:1467678/1‑100 (MQ=255)
                                                                                            gggTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACg         >  2:949597/1‑100 (MQ=255)
                                                                                             ggTTAAGCCCGAAAGAAAACCGTGAACTCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACgg        >  1:2320612/1‑100 (MQ=255)
                                                                                             ggTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACgg        >  1:393352/1‑100 (MQ=255)
                                                                                             ggTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACgg        <  2:2352923/100‑1 (MQ=255)
                                                                                             ggTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACg         <  2:3241266/99‑1 (MQ=255)
                                                                                              gTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACggt       <  1:1935183/100‑1 (MQ=255)
                                                                                               ttAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACggtg      <  1:484426/100‑1 (MQ=255)
                                                                                                  aGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTa   <  1:255192/100‑1 (MQ=255)
                                                                                                   gCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTa   >  1:3162748/1‑99 (MQ=255)
                                                                                                   gCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTAt  >  1:1389315/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 < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: