Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 964,148 T→C 15.2% intergenic (+36/‑172) 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,1480TC15.2% 368.6 / 56.3 158intergenic (+36/‑172)ihfB/ycaIintegration host factor (IHF), DNA‑binding protein, beta subunit/ComEC family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (62/72);  new base C (11/13);  total (73/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 = 3.68e-01

Only 100 of 161 total aligned reads displayed.

AAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCT  >  NC_000913/964049‑964247
                                                                                                   |                                                                                                   
aaaaTACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCt                                                                                                     >  1:1430731/1‑100 (MQ=255)
    tACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTc                                                                                                  <  2:1174415/99‑1 (MQ=255)
      cGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCggg                                                                                               <  1:3229830/100‑1 (MQ=255)
       gTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGt                                                                                              >  1:335282/1‑100 (MQ=255)
          ccTCACTTTAAACCTGGTAAAGCACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAACGCACCTGTCAGGTGCt                                                                                           >  2:2491988/1‑100 (MQ=255)
          ccTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCt                                                                                           >  2:852982/1‑100 (MQ=255)
          ccTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCt                                                                                           >  1:3507916/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)
           cTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCCGACAGGTGCtt                                                                                          <  1:4141129/100‑1 (MQ=255)
           cTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAAGATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCCGACAGGTGCtt                                                                                          <  1:2568394/100‑1 (MQ=255)
            tCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCttt                                                                                         <  1:1676668/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)
              aCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCttttt                                                                                       >  2:1764188/1‑100 (MQ=255)
              aCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCCGACAGGTGCttttc                                                                                       <  2:3993284/100‑2 (MQ=255)
               cTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCtttttt                                                                                      <  1:4026151/100‑1 (MQ=255)
               cTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCtttttt                                                                                      >  2:3682385/1‑100 (MQ=255)
                tttAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTc                                                                                     <  1:3306337/100‑1 (MQ=255)
                tttAAACCTGGAAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTc                                                                                     <  2:2364946/100‑1 (MQ=255)
                 ttAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCCGACAGGTGCTTTtctct                                                                                    <  1:406378/100‑2 (MQ=38)
                 ttAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCCGACAGGTGCTTTtctct                                                                                    <  2:879848/100‑2 (MQ=38)
                  tAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCAt                                                                                   >  1:3489367/1‑100 (MQ=255)
                   aaaCCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCAtt                                                                                  <  1:2308147/100‑1 (MQ=255)
                   aaaCCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCCGACAGGTGCttttctctcg                                                                                  <  1:2038938/100‑7 (MQ=38)
                    aaCCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCAttt                                                                                 >  1:3778251/1‑100 (MQ=255)
                     aCCTGTTAAAGAACTGCGCGATCGCGCCAATATTTACTGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTtc                                                                                <  1:3722978/100‑1 (MQ=255)
                     aCCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTtc                                                                                <  1:702473/100‑1 (MQ=255)
                     aCCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTtc                                                                                <  1:3655609/100‑1 (MQ=255)
                     aCCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCCGACAGGTGCttttctctcgt                                                                                 <  2:723266/99‑8 (MQ=37)
                       cTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTtctc                                                                              >  2:1258617/1‑100 (MQ=255)
                       cTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTtctc                                                                              <  2:2077980/100‑1 (MQ=255)
                       cTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTtctc                                                                              <  2:3417673/100‑1 (MQ=255)
                        tGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTctct                                                                             >  1:3335522/1‑100 (MQ=255)
                        tGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTctct                                                                             <  2:2803022/100‑1 (MQ=255)
                         ggTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTa                                                                            <  1:1635475/100‑1 (MQ=255)
                          gTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTaa                                                                           >  1:827260/1‑100 (MQ=255)
                          gTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTaa                                                                           >  2:2313227/1‑100 (MQ=255)
                           tAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAAt                                                                          >  1:2150984/1‑100 (MQ=255)
                            aaaGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATc                                                                         <  2:179491/100‑1 (MQ=255)
                             aaGGACTGCGCGATCGCTCCAATATTGAAGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGGGCTTTTTTCATTTCTCTAAtn                                                                         <  2:1554439/99‑2 (MQ=255)
                             aaGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATct                                                                        <  2:1191010/100‑1 (MQ=255)
                               gAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATctgg                                                                      <  1:1752557/100‑1 (MQ=255)
                                attggcgcGATCGCGCCAATATTTACAGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCAACTGTCGGGTGCTTTTTTCATTTCTCTAATct                                                                        >  2:928407/5‑97 (MQ=255)
                                aaCTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATctgga                                                                     >  2:825655/1‑100 (MQ=255)
                                 aCTGCGCGATCGAGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATctgga                                                                     <  1:2717151/99‑1 (MQ=255)
                                   tGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATctggaact                                                                  >  1:2460140/1‑100 (MQ=255)
                                    gcgcGATCGCGCCCATACTGACGGTTAAGTTTTTTACTCACACGTGAACGAGAGAAGAGCACCTCTCGGGTGCTTTTTTTATTTCTCTAATcgtga                                                                     >  2:1565584/1‑92 (MQ=255)
                                    gcgcGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATctggaactg                                                                 >  2:52488/1‑100 (MQ=255)
                                    gcgcGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATctggaactg                                                                 >  2:3126782/1‑100 (MQ=255)
                                      gcgATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATctggaactgga                                                               <  1:3819135/100‑1 (MQ=255)
                                       cgATCGCGCCAATCTTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCAACTGTCGGGTGCTTTTTTCATTTCTCTAATctggaactgga                                                               <  2:3341186/99‑1 (MQ=255)
                                                cAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATATGGAAATGGAAGCTGCCTc                                                      >  2:66926/1‑99 (MQ=255)
                                                   tatTTACGGTTAAGTCTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCag                                                  <  1:4136108/100‑1 (MQ=255)
                                                    atTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCagag                                                >  1:3838876/1‑100 (MQ=255)
                                                      ttACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGGACTGGAAGCTGCCTCGCagag                                                >  1:1517061/1‑99 (MQ=255)
                                                       tACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCGGCCTCGCagagct                                              <  1:2825398/100‑3 (MQ=255)
                                                        aCGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGttt                                             <  2:3860660/100‑1 (MQ=255)
                                                          ggTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTg                                           <  1:1232806/100‑1 (MQ=255)
                                                          ggTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTg                                           <  1:2626485/100‑1 (MQ=255)
                                                          ggTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTg                                           <  1:56224/100‑1 (MQ=255)
                                                           gTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGa                                          >  1:780808/1‑100 (MQ=255)
                                                           gTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGACGGGTGCTTTTTTCATTTCTCTAATCTGGACCTGGAAGCTGCCTCGCAGAGTTTTGa                                          >  1:81200/1‑100 (MQ=255)
                                                            ttAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGaa                                         >  2:3787959/1‑100 (MQ=255)
                                                             tAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAAc                                        <  1:3126796/100‑1 (MQ=255)
                                                               aGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAg                                      >  2:2880614/1‑100 (MQ=255)
                                                                       aCTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCAcc                              <  1:3387380/100‑1 (MQ=255)
                                                                          cAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCtt                           >  1:3220086/1‑100 (MQ=255)
                                                                            aaCTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTc                         <  1:2428034/100‑1 (MQ=255)
                                                                            aaCTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTc                         <  1:2428262/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)
                                                                                 gAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTaa                    >  2:3709126/1‑100 (MQ=255)
                                                                                 gAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTaa                    <  2:3061904/100‑1 (MQ=255)
                                                                                   aCGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAAt                  <  2:2346282/100‑1 (MQ=255)
                                                                                   aCGAGAGAAAAGCACCTGTCGGGTGCTTTTTGCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAAt                  <  2:1576465/100‑1 (MQ=255)
                                                                                    cGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAAtt                 >  2:3211213/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)
                                                                                      tgaaaaaaGCACCCGACAGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAAttct               >  1:3961426/5‑100 (MQ=38)
                                                                                      tgaaaaaaGCACCCGACAGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGAAGAGTATTGAACAGTTTTCACCCTTTCGTTAAAt                  >  2:73551/5‑97 (MQ=37)
                                                                                       gagaAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAAttctt              >  1:1361599/1‑100 (MQ=255)
                                                                                       gagaAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAAttctt              >  2:2346055/1‑100 (MQ=255)
                                                                                       gaaaaaaGCACCCGACAGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAAttctt              >  1:4154775/4‑100 (MQ=38)
                                                                                        agaAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAAttcttc             <  1:50020/100‑1 (MQ=255)
                                                                                         gaAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATtcttct            <  2:3038630/100‑1 (MQ=255)
                                                                                         gaAAAGCACCTGTCGGGCGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGATGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATtcttct            >  2:1338552/1‑100 (MQ=255)
                                                                                         aaaaaGCACCCGACAGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATtcttct            >  1:3811344/2‑100 (MQ=255)
                                                                                           aaaGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTTGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTg           <  2:3674717/99‑1 (MQ=255)
                                                                                           aaaGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGa          >  2:4081646/1‑100 (MQ=255)
                                                                                              gCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAat        <  1:1740419/99‑1 (MQ=255)
                                                                                               cACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAatat      >  1:3367166/1‑100 (MQ=255)
                                                                                               cACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAatat      >  1:1006387/1‑100 (MQ=255)
                                                                                                 ccTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGc    <  2:2711190/100‑1 (MQ=255)
                                                                                                   tGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCt  >  1:2569855/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
AAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCT  >  NC_000913/964049‑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: