Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 964,150 T→A 15.8% intergenic (+38/‑170) ihfB → / → ycaI integration host factor (IHF), DNA‑binding protein, beta subunit/ComEC family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913964,1500TA15.8% 393.8 / 61.5 159intergenic (+38/‑170)ihfB/ycaIintegration host factor (IHF), DNA‑binding protein, beta subunit/ComEC family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (62/70);  new base A (12/13);  total (74/85)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.29e-01

Only 100 of 159 total aligned reads displayed.

TACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCT  >  NC_000913/964053‑964247
                                                                                                 |                                                                                                 
  cGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCggg                                                                                               <  1:3229830/100‑1 (MQ=255)
   gTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGt                                                                                              >  1:335282/1‑100 (MQ=255)
     tCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGc                                                                                            <  2:2342914/100‑1 (MQ=255)
      ccTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGc                                                                                            <  1:3203695/99‑1 (MQ=255)
      ccTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCt                                                                                           >  1:2724694/1‑100 (MQ=255)
      ccTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCt                                                                                           >  1:3507916/1‑100 (MQ=255)
      ccTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCt                                                                                           >  2:852982/1‑100 (MQ=255)
       cTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCtt                                                                                          >  1:336133/1‑100 (MQ=255)
       cTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCtt                                                                                          <  2:1434376/100‑1 (MQ=255)
       cTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCt                                                                                           >  2:3946005/1‑99 (MQ=255)
       cTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCCGACAGGTGCtt                                                                                          <  1:4138957/100‑1 (MQ=255)
       cTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAAGATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCCGACAGGTGCtt                                                                                          <  1:2568394/100‑1 (MQ=255)
         cagtttAAAGCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTGTTTTACTTTAACTTGAACGAGAGAAAAGCACCTGTGGTGGGCtttt                                                                                        >  2:2481619/4‑100 (MQ=255)
         cACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCtttt                                                                                        >  1:2755704/1‑100 (MQ=255)
         cACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCCGACAGGTGCttag                                                                                        >  2:4141100/1‑98 (MQ=255)
         cACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCCGACAGGTGCttag                                                                                        >  2:4138928/1‑98 (MQ=255)
          aCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCttttt                                                                                       >  2:1764188/1‑100 (MQ=255)
          aCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCttttt                                                                                       <  2:2346716/100‑1 (MQ=255)
          aCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCCGACAGGTGCttttc                                                                                       <  2:3993284/100‑2 (MQ=255)
          aCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCCGACAGGTGCttttc                                                                                       <  1:4073303/100‑2 (MQ=255)
           cTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCtttttt                                                                                      >  2:3682385/1‑100 (MQ=255)
           cTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCtttttt                                                                                      <  1:4026151/100‑1 (MQ=255)
            tttAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTc                                                                                     <  1:3306337/100‑1 (MQ=255)
             ttAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCCGACAGGTGCTTTtctct                                                                                    <  1:406378/100‑2 (MQ=38)
             ttAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCCGACAGGTGCTTTtctct                                                                                    <  2:879848/100‑2 (MQ=38)
             ttAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCCGACAGGTGCTTTtctct                                                                                    <  2:3065473/100‑2 (MQ=38)
             ttAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCCGACAGGTGCTTTtctct                                                                                    <  2:1475185/100‑2 (MQ=38)
              tAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCAt                                                                                   >  1:3489367/1‑100 (MQ=255)
               aaaCCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCAtt                                                                                  <  1:2308147/100‑1 (MQ=255)
                aaCCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCAttt                                                                                 >  1:3778251/1‑100 (MQ=255)
                 aCCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTtc                                                                                <  2:3794046/100‑1 (MQ=255)
                 aCCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTtc                                                                                <  2:757789/100‑1 (MQ=255)
                 aCCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTtc                                                                                <  1:3655609/100‑1 (MQ=255)
                 aCCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCCGACAGGTGCttttctctcgt                                                                                 <  2:723266/99‑8 (MQ=37)
                   cTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTtctc                                                                              <  2:3417673/100‑1 (MQ=255)
                    tGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTctct                                                                             <  2:2634606/100‑1 (MQ=255)
                    tGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTctct                                                                             >  1:3335522/1‑100 (MQ=255)
                    tGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTctct                                                                             >  1:1627244/1‑100 (MQ=255)
                     ggTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTa                                                                            <  1:1635475/100‑1 (MQ=255)
                     ggTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTa                                                                            >  1:630097/1‑100 (MQ=255)
                       tAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAAt                                                                          >  1:2150984/1‑100 (MQ=255)
                         aaGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATct                                                                        <  2:1191010/100‑1 (MQ=255)
                         aaGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATct                                                                        <  2:2810218/100‑1 (MQ=255)
                           gAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATctgg                                                                      >  1:3852776/1‑100 (MQ=255)
                           gAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATctgg                                                                      <  1:1752557/100‑1 (MQ=255)
                            attggcgcGATCGCGCCAATATTTACAGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCAACTGTCGGGTGCTTTTTTCATTTCTCTAATct                                                                        >  2:928407/5‑97 (MQ=255)
                            aaCTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATctgga                                                                     >  2:825655/1‑100 (MQ=255)
                            aaCTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATctgg                                                                      <  2:3397343/99‑1 (MQ=255)
                              cTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATctggaac                                                                   <  1:825655/100‑1 (MQ=255)
                              cTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATctggaac                                                                   <  1:1692913/100‑1 (MQ=255)
                               tGCGCGATCGTGACAATAGTTACGGTCAAGTTCTTTAGTCAAACTTGAACGAGAGAAAAGCAGCTTTCGGGTGCACTTTTCATTTCTCTAATctgcaa                                                                    <  2:167326/98‑1 (MQ=255)
                               tGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATctggaact                                                                  >  1:2460140/1‑100 (MQ=255)
                                gcgcGATCGCGCCCATACTGACGGTTAAGTTTTTTACTCACACGTGAACGAGAGAAGAGCACCTCTCGGGTGCTTTTTTTATTTCTCTAATcgtga                                                                     >  2:1565584/1‑92 (MQ=255)
                                gcgcGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATctggaactg                                                                 >  2:52488/1‑100 (MQ=255)
                                 cgcgATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATctggaactgg                                                                >  1:1863856/1‑100 (MQ=255)
                                  gcgATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATctggaactgga                                                               <  1:3819135/100‑1 (MQ=255)
                                         cgcCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGc                                                         <  2:1923415/99‑1 (MQ=255)
                                            cAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCg                                                     <  1:2443166/100‑1 (MQ=255)
                                            cAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCg                                                     >  2:288644/1‑100 (MQ=255)
                                            cAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATATGGAAATGGAAGCTGCCTc                                                      >  2:66926/1‑99 (MQ=255)
                                             aaTATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGc                                                    <  1:791016/100‑1 (MQ=255)
                                                atTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCagag                                                >  1:3838876/1‑100 (MQ=255)
                                                   tACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCGGCCTCGCagagct                                              <  1:2825398/100‑3 (MQ=255)
                                                    aCGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGttt                                             <  2:3860660/100‑1 (MQ=255)
                                                     cGGTTAAGTTTTTTTCTTTATTTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGtttt                                            >  2:121210/1‑100 (MQ=255)
                                                     cGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGGCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGtttt                                            <  2:157016/100‑1 (MQ=255)
                                                      ggTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTg                                           <  2:2798677/100‑1 (MQ=255)
                                                      ggTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTg                                           <  1:2626485/100‑1 (MQ=255)
                                                      ggTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTg                                           <  1:1232806/100‑1 (MQ=255)
                                                       gTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGa                                          >  1:780808/1‑100 (MQ=255)
                                                         tAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAAc                                        <  1:3126796/100‑1 (MQ=255)
                                                         tAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAAc                                        <  2:3138587/100‑1 (MQ=255)
                                                            gTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGt                                     >  1:3718987/1‑100 (MQ=255)
                                                                   aCTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCAcc                              <  2:22459/100‑1 (MQ=255)
                                                                   aCTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCAcc                              >  1:4189947/1‑100 (MQ=255)
                                                                   aCTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCAcc                              <  1:3387380/100‑1 (MQ=255)
                                                                    cTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCAccc                             >  2:2460965/1‑100 (MQ=255)
                                                                        aaCTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTc                         <  1:2428034/100‑1 (MQ=255)
                                                                        aaCTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTc                         >  1:3022280/1‑100 (MQ=255)
                                                                         aCTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCg                        >  1:598080/1‑100 (MQ=255)
                                                                          cTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGt                       <  2:2306693/100‑1 (MQ=255)
                                                                            tGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTa                     <  1:1434145/100‑1 (MQ=255)
                                                                             gAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTaa                    >  2:3709126/1‑100 (MQ=255)
                                                                                cGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAAtt                 >  2:3211213/1‑100 (MQ=255)
                                                                                cGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAAtt                 >  2:1268533/1‑100 (MQ=255)
                                                                                aatgaaaaaaGCACCCGACAGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAAtt                 >  2:3806813/7‑100 (MQ=38)
                                                                                 gtgagaAACGCACCTGTCGGGTGCTTTTTTCATTTCTCTAGTCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCCCCCTTCCGTTAAAt                  <  2:2686173/96‑1 (MQ=255)
                                                                                 gagagaAAAGCACCTGGCGGGTGCTTGTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAAttc                <  2:101114/100‑1 (MQ=255)
                                                                                  tgaaaaaaGCACCCGACAGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGGAGAGTGTTGAACAGGTTTCACCCTTTCGTTAAAttct               >  2:2835034/5‑100 (MQ=37)
                                                                                   gagaAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAAttctt              >  2:2346055/1‑100 (MQ=255)
                                                                                   gagaAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAAttctt              >  1:1361599/1‑100 (MQ=255)
                                                                                   gagaAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGGTTTGAACAGTTTTCACCCTTTCGTTAAAttctt              >  1:539121/1‑100 (MQ=255)
                                                                                   gaaaaaaGCACCCGACAGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAAttctt              >  1:4154775/4‑100 (MQ=38)
                                                                                    aaaaaaGCACCCGACAGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAAttcttc             >  2:1131576/3‑100 (MQ=255)
                                                                                     gaAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATtcttct            <  2:3038630/100‑1 (MQ=255)
                                                                                     aaaaaGCACCCGACAGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATtcttct            >  1:3811344/2‑100 (MQ=255)
                                                                                       aaaGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTTGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTg           <  2:3674717/99‑1 (MQ=255)
                                                                                          gCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAata       >  1:4139784/1‑100 (MQ=255)
                                                                                           cACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAatat      >  1:1006387/1‑100 (MQ=255)
                                                                                               tGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCt  >  1:2569855/1‑100 (MQ=255)
                                                                                                 |                                                                                                 
TACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCT  >  NC_000913/964053‑964247

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: