Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,262,250 G→C 15.1% G194G (GGG→GGC dusA → tRNA‑dihydrouridine synthase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,262,2500GC15.1% 506.1 / 77.3 186G194G (GGG→GGCdusAtRNA‑dihydrouridine synthase A
Reads supporting (aligned to +/- strand):  ref base G (91/65);  new base C (14/14);  total (106/80)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.65e-01

Only 100 of 196 total aligned reads displayed.

GGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGG  >  NC_000913/4262151‑4262348
                                                                                                   |                                                                                                  
ggACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCgggg                                                                                                    <  1:4399477/100‑1 (MQ=255)
ggACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCgggg                                                                                                    >  1:3445632/1‑100 (MQ=255)
  aCAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGtt                                                                                                  <  2:2525736/100‑1 (MQ=255)
   cAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTa                                                                                                 >  1:4358636/1‑100 (MQ=255)
     gCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAg                                                                                               <  2:772905/100‑1 (MQ=255)
      cTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGAAACAAGCGCGTCTGAGATATTCATCATAAACGCACATAAAGACTGGCTTTCGGGCTTAAg                                                                                               >  1:3734539/1‑99 (MQ=255)
       tATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGcc                                                                                             >  2:3608286/1‑100 (MQ=255)
        aTGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGccc                                                                                            >  2:3251222/1‑100 (MQ=255)
        aTGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGccc                                                                                            <  1:737034/100‑1 (MQ=255)
         tGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCg                                                                                           >  1:3572473/1‑100 (MQ=255)
          gAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCga                                                                                          >  1:3614917/1‑100 (MQ=255)
           aaTTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCgaa                                                                                         >  2:1443929/1‑100 (MQ=255)
            aTTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCgaaa                                                                                        <  2:487987/100‑1 (MQ=255)
            aTTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCgaaa                                                                                        <  1:3394709/100‑1 (MQ=255)
            aTTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCgaaa                                                                                        <  1:3604675/100‑1 (MQ=255)
            aTTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaa                                                                                        <  2:485998/100‑1 (MQ=255)
             tttCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaag                                                                                       <  2:1240520/100‑1 (MQ=255)
             tttCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaag                                                                                       <  2:3753897/100‑1 (MQ=255)
             tttCTCTGCGATTTCATAAATACCGTTTCCGGCAAAGGCGAGTATGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCgaaag                                                                                       <  1:2529744/100‑1 (MQ=255)
              ttCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaagc                                                                                      <  2:1368002/100‑2 (MQ=255)
                ctctGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCgaaagaaa                                                                                    >  1:620217/1‑100 (MQ=255)
                 tctGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaagccag                                                                                   <  1:2355596/100‑5 (MQ=39)
                 tctGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaagccag                                                                                   <  2:2936741/100‑5 (MQ=39)
                  ctGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTGAGCCCGAAAGAAAAc                                                                                  >  1:1002359/1‑100 (MQ=255)
                  ctGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAAc                                                                                  <  2:2883316/100‑1 (MQ=255)
                  ctGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAAc                                                                                  >  2:3046736/1‑100 (MQ=255)
                  ctGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaagccagg                                                                                  <  1:3228664/100‑6 (MQ=39)
                  ctGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaagccagg                                                                                  <  2:4274635/100‑6 (MQ=39)
                  ctGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGATTTCGGGGTTAAGACCGAAAGAAAAc                                                                                  >  1:1002136/1‑100 (MQ=255)
                    gCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaagccaggct                                                                                <  1:2776597/100‑8 (MQ=38)
                     cGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGt                                                                               >  2:2592360/1‑100 (MQ=255)
                      gATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTg                                                                              <  2:3862795/100‑1 (MQ=255)
                      gATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTg                                                                              <  1:3246738/100‑1 (MQ=255)
                      gATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTg                                                                              >  1:1709178/1‑100 (MQ=255)
                       aTTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGa                                                                             <  1:4411259/100‑1 (MQ=255)
                        tttCATCCATACCGTTTCCGGCAAAGGCGAGTGTGCGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGaa                                                                            <  2:2403236/100‑1 (MQ=255)
                        tttCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGaa                                                                            <  2:2736420/100‑1 (MQ=255)
                         ttCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGaaa                                                                           <  1:3148733/100‑1 (MQ=255)
                           catcaATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATc                                                                         >  2:899338/1‑100 (MQ=255)
                           catcaATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATc                                                                         >  1:263458/1‑100 (MQ=255)
                            atcaATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATcc                                                                        >  2:2482473/1‑100 (MQ=255)
                             tcaATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATccc                                                                       >  1:1663547/1‑100 (MQ=255)
                             tcaATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATccc                                                                       >  1:1865793/1‑100 (MQ=255)
                               aaTACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCccgc                                                                     >  2:3577458/1‑100 (MQ=255)
                               aaTACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCccgc                                                                     <  2:2911588/100‑1 (MQ=255)
                               aaTACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCccgc                                                                     <  2:3996617/100‑1 (MQ=255)
                               aaTACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCccgc                                                                     <  2:1761264/100‑1 (MQ=255)
                               aaTACCGTTCCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATccc                                                                       <  1:1106226/98‑1 (MQ=255)
                                aTACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCccgcc                                                                    >  1:3516321/1‑100 (MQ=255)
                                 tACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCccgccg                                                                   <  1:3592960/100‑1 (MQ=255)
                                       ttCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGAt                                                             >  2:2855438/1‑100 (MQ=255)
                                        tCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGAtt                                                            <  2:1326876/100‑1 (MQ=255)
                                        tCCGGCAAAGGCGAGAGTGAGATGTTCATCATCCCCGCACGTAAAGCCTGGCTTTAGGGTTTTAGCCCGAAGGAAAACCGTGAAATccgaccg                                                                   >  1:3675347/1‑88 (MQ=255)
                                         ccGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTa                                                           >  1:4226012/1‑100 (MQ=255)
                                          cGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTAt                                                          >  2:1951545/1‑100 (MQ=255)
                                           ggCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATc                                                         >  2:2178026/1‑100 (MQ=255)
                                           ggCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATc                                                         >  1:3841833/1‑100 (MQ=255)
                                           ggCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCACCGCTCGATTATc                                                         >  1:2243442/1‑100 (MQ=255)
                                             cAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCcg                                                       <  2:2054235/100‑1 (MQ=255)
                                                 ggCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCgtg                                                   >  2:1331861/1‑100 (MQ=255)
                                                  gCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCgtgt                                                  >  2:3520202/1‑100 (MQ=255)
                                                   cGAGTGTGGGATGTTCAGCATCCACGTACGTAAACCCTGGCCGTCGAGGTTAACCACGAAAGAAAACCGTGAAATCCCGCCGGTCGATTAGCCGCgggtg                                                 <  1:1469906/100‑1 (MQ=255)
                                                   cGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCgtgtg                                                 <  1:3046718/100‑1 (MQ=255)
                                                          gagaTGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAAc                                          <  2:1619650/100‑1 (MQ=255)
                                                           agaTGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACt                                         <  1:3192981/100‑1 (MQ=255)
                                                             aTGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGa                                       <  1:4059493/100‑1 (MQ=255)
                                                              tGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGaa                                      >  2:146830/1‑100 (MQ=255)
                                                                 tcatcaTCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCg                                   >  2:2643367/1‑100 (MQ=255)
                                                                 tcatcaTCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCg                                   >  1:1077987/1‑100 (MQ=255)
                                                                 tcatcaTCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCg                                   >  2:174249/1‑100 (MQ=255)
                                                                      atcCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACt                              >  2:695514/1‑100 (MQ=255)
                                                                      atcCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACt                              >  2:231214/1‑100 (MQ=255)
                                                                      atcCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACt                              >  2:2163291/1‑100 (MQ=255)
                                                                         cacgcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTc                           >  2:614715/1‑100 (MQ=255)
                                                                         cacgcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTc                           >  1:2918203/1‑100 (MQ=255)
                                                                            gcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGc                        <  1:1939917/100‑1 (MQ=255)
                                                                                 tAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTg                   <  2:105298/100‑1 (MQ=255)
                                                                                 tAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTg                   >  2:3378396/1‑100 (MQ=255)
                                                                                  aaaGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGa                  >  2:1503145/1‑100 (MQ=255)
                                                                                   acggttttcTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGAc                 >  1:167318/9‑100 (MQ=38)
                                                                                   aaGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGAc                 >  1:153974/1‑100 (MQ=255)
                                                                                    aGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACa                >  1:1337790/1‑100 (MQ=255)
                                                                                      ccTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAAt              <  2:355951/100‑1 (MQ=255)
                                                                                       cTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATg             <  2:3213272/100‑1 (MQ=255)
                                                                                        tGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGt            <  1:1503139/100‑1 (MQ=255)
                                                                                        tGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGt            <  2:2020220/100‑1 (MQ=255)
                                                                                         ggCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCTCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTc           <  1:2816254/100‑1 (MQ=255)
                                                                                         ggCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAAGCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTc           <  2:4208335/100‑1 (MQ=255)
                                                                                          tcTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCg          >  2:4351971/2‑100 (MQ=255)
                                                                                          gCTTTCGGGGTTAAGCCCGAAAGAAAACCCTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCACCTGAAGCGTGACTTGCCGCATCTGACAATGTCg          <  1:2649359/100‑1 (MQ=255)
                                                                                           cTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGa         >  1:138012/1‑100 (MQ=255)
                                                                                            tttCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGAt        <  2:164315/100‑1 (MQ=255)
                                                                                            tttCGGGGTTAAGCCCGAAAGAAAAACGTGAAATCACGCCGCTCGATTATCCGCGGGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGAt        >  1:3291878/1‑100 (MQ=255)
                                                                                             ttCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGAtt       <  2:1672741/100‑1 (MQ=255)
                                                                                                ggggTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAAc    >  1:891245/1‑100 (MQ=255)
                                                                                                ggggTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTACCAGCTGAAGAGTGACTTTCCGCATCTGACAATGTCGATTAAc    >  1:3093462/1‑100 (MQ=255)
                                                                                                 gggTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACg   >  2:4111025/1‑100 (MQ=255)
                                                                                                 gggTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACg   >  2:4316532/1‑100 (MQ=255)
                                                                                                 gggTTAAGCCCCAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAACGTGACTTTCCGCATCTGACAATGTCGATTAACg   >  1:2318246/1‑100 (MQ=255)
                                                                                                  ggTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACgg  <  2:2960709/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
GGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGG  >  NC_000913/4262151‑4262348

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