Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 120,158:1 +A 17.3% intergenic (+23/+20) ampE → / ← aroP ampicillin resistance inner membrane protein; putative signaling protein in beta‑lactamase regulation/aromatic amino acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913120,1581.A17.3% 578.4 / 146.8 196intergenic (+23/+20)ampE/aroPampicillin resistance inner membrane protein; putative signaling protein in beta‑lactamase regulation/aromatic amino acid transporter
Reads supporting (aligned to +/- strand):  ref base . (73/89);  new base A (14/20);  total (87/109)
Fisher's exact test for biased strand distribution p-value = 7.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.17e-01

Only 100 of 206 total aligned reads displayed.

TTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCA  >  NC_000913/120059‑120257
                                                                                                    |                                                                                                   
tttCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAc                                                                                                      <  2:1072613/100‑1 (MQ=255)
 ttCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑g                                                                                                    <  2:1756389/100‑1 (MQ=255)
 ttCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑g                                                                                                    <  2:3774532/100‑1 (MQ=255)
  tCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGACTGCACTACTGACGCTTTACGGGGCGTTGGTGTAAAGATTATTGCCCTAACCCTGTAC‑gg                                                                                                   >  1:4245479/1‑100 (MQ=255)
   cAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑ggg                                                                                                  <  2:3971802/100‑1 (MQ=255)
    aaTGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGt                                                                                                 >  2:2848328/1‑100 (MQ=255)
     aTGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTg                                                                                                <  1:729520/100‑1 (MQ=255)
      tGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGa                                                                                               <  1:2130052/100‑1 (MQ=255)
       ggCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAg                                                                                              >  1:2761772/1‑100 (MQ=255)
       ggCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAg                                                                                              >  2:103839/1‑100 (MQ=255)
       ggCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCC‑GTACAGGGTGAg                                                                                              <  2:1534278/100‑1 (MQ=255)
        gCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCC‑GTACAGGGTGAgg                                                                                             <  1:138312/100‑1 (MQ=255)
         cGAAGAAATCCTCGTTCGTGGTCGGGGTGCTGATTGCACTACTGACGATTTACGGGGCGTTGGGGTCAAGATTATTGCCCTCACCC‑GTACAGGGTGAggg                                                                                            <  1:2638095/100‑1 (MQ=255)
         cGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTTGTGTAAAGATTATTGCCCTCACCCAGTAC‑GGGTGAggg                                                                                            <  1:2708905/100‑1 (MQ=255)
         cGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAggg                                                                                            <  2:1796924/100‑1 (MQ=255)
          gaagaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCC‑GTACAGGGTGAGGGc                                                                                           <  1:2451696/100‑1 (MQ=255)
          gaagaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCC‑GTACAGGGTGAGGGc                                                                                           <  2:3660599/100‑1 (MQ=255)
          gaagaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCC‑GTACAGGGTGAGGGc                                                                                           <  2:3583370/100‑1 (MQ=255)
          gaagaaAACCTCGTTCGTGGTCGTAGTGGTGATTGCACTACTGGCGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCC‑GTACAGGGTGAGGGc                                                                                           <  2:775415/100‑1 (MQ=255)
           aagaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCg                                                                                          <  2:2815417/100‑1 (MQ=255)
              aaaaCCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTag                                                                                       >  1:2363370/1‑100 (MQ=255)
              aaaaCCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTag                                                                                       <  2:197986/100‑1 (MQ=255)
               aaaCCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTaga                                                                                      <  2:228650/100‑1 (MQ=255)
               aaaCCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATCGCCCTCACCC‑GTACAGGGTGAGGGcaataa                                                                                      <  1:937173/100‑6 (MQ=37)
                 aCCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTagaga                                                                                    <  2:167686/100‑1 (MQ=255)
                 aCCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTagaga                                                                                    >  1:389127/1‑100 (MQ=255)
                  ccTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTagagag                                                                                   >  2:831371/1‑100 (MQ=255)
                  ccTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATCATTGCCCTCACCC‑GTACAGGGTGAGGGcaataatct                                                                                   <  2:3111198/100‑9 (MQ=37)
                    tcgtttgtGGTCGTGGTGGTGATTGCTTTAGTGACGATTTACGGGGCGTGGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTagaga                                                                                    <  1:1456883/97‑1 (MQ=255)
                    tcgttcgtGGTCGTGGTGGTGATTGGGCTCCTGACGATTTACGTGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTagag                                                                                     <  1:3571544/96‑1 (MQ=255)
                    tcgttcgtGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGAt                                                                                 <  2:3805443/100‑1 (MQ=255)
                    tcgttcgtGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGAt                                                                                 <  2:2885702/100‑1 (MQ=255)
                     cgttcgtGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGAtt                                                                                <  2:1577178/100‑1 (MQ=255)
                         cgtggtcgtggtGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATg                                                                            <  2:1855016/100‑1 (MQ=255)
                          gtggtcgtggtGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATgc                                                                           >  1:650328/1‑100 (MQ=255)
                          ggggtagtggtGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATg                                                                            <  1:2593696/97‑1 (MQ=255)
                           tggtcgtggtGGTGATTGCACTACTTACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATgcg                                                                          <  1:2069018/100‑1 (MQ=255)
                            ggtcgtggtGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATgcgc                                                                         <  2:2041808/100‑1 (MQ=255)
                                 tggtggtgATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCttttt                                                                    >  2:1974014/1‑99 (MQ=255)
                                 tggtggtgATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCCGTAC‑GGGTGAGGGCGTAGAGAGACTAATGCGCTTTTa                                                                    >  1:695894/1‑100 (MQ=255)
                                  ggtggtgATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTAc                                                                   >  1:2864845/1‑100 (MQ=255)
                                       tgATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGCGAGGGCGTAGAGAGATTAATGCCCTTTTGCGGCtt                                                              >  1:3453282/1‑100 (MQ=255)
                                         aTTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTg                                                            <  1:2191977/100‑1 (MQ=255)
                                          ttGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTgg                                                           <  2:1604005/100‑1 (MQ=255)
                                           tGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGc                                                          <  2:1833678/100‑1 (MQ=255)
                                           tGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGc                                                          >  1:1837684/1‑100 (MQ=255)
                                            gCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCACTCACCCTGTAC‑GGGTGAGGGCGTAGATAGATTAATGCGCTTTTAAGGCTTTGGct                                                         >  1:2665031/1‑99 (MQ=255)
                                             cACTACTGACGATTTACGGGGCGTTGGTTTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCgg                                                        <  2:4142984/100‑1 (MQ=255)
                                               ctactGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGtt                                                      <  2:3604567/100‑1 (MQ=255)
                                                 actGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGgttt                                                    >  1:539843/1‑100 (MQ=255)
                                                  ctGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTtc                                                   <  2:3040842/100‑1 (MQ=255)
                                                   tGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTtct                                                  >  2:1729138/1‑100 (MQ=255)
                                                    gACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTtctc                                                 >  2:3036456/1‑100 (MQ=255)
                                                     aCGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTctct                                                <  2:637638/100‑1 (MQ=255)
                                                     aCGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTctct                                                <  1:3199883/100‑1 (MQ=255)
                                                     aCGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTctct                                                <  2:4107268/100‑1 (MQ=255)
                                                        aTTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTa                                             <  1:285573/100‑1 (MQ=255)
                                                         tttACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGTCGGCTTTCCCTTGaa                                            >  1:1955392/1‑100 (MQ=255)
                                                           tACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAAc                                          <  1:3587953/100‑1 (MQ=255)
                                                            aCGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACa                                         >  2:1322055/1‑100 (MQ=255)
                                                            aCGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACa                                         <  1:3127186/100‑1 (MQ=255)
                                                               gggCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGCTTAATGCGCTTTTACGGCTTTGGCGGGTTTGTCTTTAAACAGAt                                      >  1:3726572/1‑100 (MQ=255)
                                                                ggCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATa                                     >  1:1064016/1‑100 (MQ=255)
                                                                ggCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATa                                     <  2:1120949/100‑1 (MQ=255)
                                                                ggCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATa                                     <  2:2761787/100‑1 (MQ=255)
                                                                ggCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCACTTTAAAAAGAtc                                     <  1:3265474/100‑2 (MQ=255)
                                                                  cGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGc                                   >  2:2102614/1‑100 (MQ=255)
                                                                   gTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGcc                                  <  2:3316475/100‑1 (MQ=255)
                                                                    ttGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGACGGCGTAGAGAGATTAATGCCCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCg                                 <  1:2398880/100‑1 (MQ=255)
                                                                    ttGGTGTAAAGATTAGTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGATAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCg                                 <  1:2517086/100‑1 (MQ=255)
                                                                       gtgtAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATa                              >  2:4257020/1‑100 (MQ=255)
                                                                       gtgtAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATa                              <  2:194864/100‑1 (MQ=255)
                                                                        tgtAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATAc                             >  2:1157687/1‑100 (MQ=255)
                                                                         gtAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATAcc                            >  1:1205127/1‑100 (MQ=255)
                                                                           aaaGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCt                           >  2:3137820/1‑99 (MQ=255)
                                                                           aaaGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTa                          <  1:2309482/100‑1 (MQ=255)
                                                                              gATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAaca                       >  2:2639699/1‑100 (MQ=255)
                                                                               ccgatctcccTCACCC‑GTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAac                        <  1:3282151/91‑1 (MQ=255)
                                                                                  attGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGCGATTAATGCGCTTTTACGGCTTTGGCGGGTTTCTCTTTACACATATAGCCGATACCtca                         >  1:4299717/1‑92 (MQ=255)
                                                                                    cgCCCTCACCC‑GTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGGTTTCTCTTTAAACCGATAGCCGATACTTAAAACGATc                  >  1:1341688/2‑99 (MQ=37)
                                                                                     gCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCa                 <  2:3726596/99‑1 (MQ=255)
                                                                                     gCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCa                 >  2:2563541/1‑99 (MQ=255)
                                                                                     gCCCTCACCC‑GTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAAGAGATAGCCGATACCTAACACGACCAg                >  1:1146781/1‑100 (MQ=255)
                                                                                     gCCCTCACCC‑GTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAg                >  2:4322065/1‑100 (MQ=255)
                                                                                      cccTCACCC‑GTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGCTACCTAACACGATCAGc               >  1:4067208/1‑100 (MQ=255)
                                                                                      cccTCACCC‑GTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGc               >  2:474520/1‑100 (MQ=255)
                                                                                       ccTCACCCTGTGC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGATTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGcc              >  2:2953398/1‑100 (MQ=255)
                                                                                       ccTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGcc              >  1:808043/1‑100 (MQ=255)
                                                                                       ccTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGcc              <  2:2345236/100‑1 (MQ=255)
                                                                                       ccTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGcc              <  1:864772/100‑1 (MQ=255)
                                                                                       ccTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGcc              <  2:1706101/100‑1 (MQ=255)
                                                                                         tCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGcc              <  2:2203246/98‑1 (MQ=255)
                                                                                            cccTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATAcc         >  1:1397548/1‑100 (MQ=255)
                                                                                               tGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCggg      <  2:1430324/100‑1 (MQ=255)
                                                                                               tGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCggg      <  2:2293631/100‑1 (MQ=255)
                                                                                                gTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGa     <  1:1097597/100‑1 (MQ=255)
                                                                                                 tAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGAt    >  2:4102527/1‑100 (MQ=255)
                                                                                                 tAC‑GGGTGAGCGCGTAGAGAGATTAATGCGCTTTTCCGGGTTTCGCGGTTTTCTCTTTAAACAGATACCCGATACCTAACACGCTCAGCCATACCGGGAt    <  1:2996274/100‑1 (MQ=255)
                                                                                                  aC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGAGAGCCGATACCTAACACGATCAGCCATACCGGGATc   <  1:3014188/100‑1 (MQ=255)
                                                                                                   c‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCa  <  2:3265521/100‑1 (MQ=255)
                                                                                                    |                                                                                                   
TTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCA  >  NC_000913/120059‑120257

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: ‑