Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 120,154 Δ1 bp 16.9% intergenic (+19/+24) 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,1540T.16.9% 571.8 / 117.0 201intergenic (+19/+24)ampE/aroPampicillin resistance inner membrane protein; putative signaling protein in beta‑lactamase regulation/aromatic amino acid transporter
Reads supporting (aligned to +/- strand):  ref base T (74/91);  new base . (14/20);  total (89/112)
Fisher's exact test for biased strand distribution p-value = 8.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.99e-01

Only 100 of 205 total aligned reads displayed.

GCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGG  >  NC_000913/120055‑120253
                                                                                                   |                                                                                                    
gCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCt                                                                                                      <  1:1221246/100‑1 (MQ=255)
  ggTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGt                                                                                                    >  1:2855607/1‑100 (MQ=255)
   gTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGt                                                                                                    <  1:783516/99‑1 (MQ=255)
    tttCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAc                                                                                                  <  2:1072817/100‑1 (MQ=255)
     ttCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑g                                                                                                <  2:1756389/100‑1 (MQ=255)
     ttCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑g                                                                                                <  2:3774532/100‑1 (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)
           ggCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCC‑GTACAGGGTGAg                                                                                          <  2:1534278/100‑1 (MQ=255)
            gCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCC‑GTACAGGGTGacg                                                                                         <  1:4367886/100‑3 (MQ=255)
            gCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCC‑GTACAGGGTGAgg                                                                                         <  1:138312/100‑1 (MQ=255)
             cGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTTGTGTAAAGATTATTGCCCTCACCCAGTAC‑GGGTGAggg                                                                                        <  1:2708905/100‑1 (MQ=255)
             cGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAggg                                                                                        >  2:4211682/1‑100 (MQ=255)
             cGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAggg                                                                                        <  2:1796924/100‑1 (MQ=255)
             cGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAggg                                                                                        <  2:949258/100‑1 (MQ=255)
             cGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCC‑GTACAGGGTGAggg                                                                                        <  2:3320089/100‑1 (MQ=255)
              gaagaaAGGCTCGTGCGTGGTCGTGGTGGGGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAggg                                                                                        <  1:3992175/99‑1 (MQ=255)
              gaagaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCC‑GTACAGGGTGAggg                                                                                        <  1:2859792/99‑1 (MQ=255)
              gaagaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCC‑GTACAGGGTGAGGGc                                                                                       <  2:3037029/100‑1 (MQ=255)
              gaagaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCC‑GTACAGGGTGAGGGc                                                                                       <  2:3583370/100‑1 (MQ=255)
              gaagaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCC‑GTACAGGGTGAGGGc                                                                                       <  2:3660599/100‑1 (MQ=255)
              gaagaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCC‑GTACAGGGTGAGGGc                                                                                       <  1:2451696/100‑1 (MQ=255)
               aagaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCg                                                                                      <  2:2815417/100‑1 (MQ=255)
                  aaaaCCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTag                                                                                   >  1:2363370/1‑100 (MQ=255)
                   aaaCCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTaga                                                                                  <  2:228650/100‑1 (MQ=255)
                     aCCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTagaga                                                                                >  2:2001634/1‑100 (MQ=255)
                     aCCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTagaga                                                                                <  2:167686/100‑1 (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                                                                            >  1:823108/1‑100 (MQ=255)
                         cgttcgtGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGAtt                                                                            <  2:1577178/100‑1 (MQ=255)
                           ttcgtGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTaa                                                                          <  2:1972361/100‑1 (MQ=255)
                               tggtcgtggtGGTGATTGCACTACTTACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATgcg                                                                      <  1:2069018/100‑1 (MQ=255)
                               tggtcgtggtGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATgcg                                                                      <  2:2080604/100‑1 (MQ=255)
                                     tggtggtgATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTa                                                                <  2:458397/100‑1 (MQ=255)
                                     tggtggtgATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCCGTAC‑GGGTGAGGGCGTAGAGAGACTAATGCGCTTTTa                                                                >  1:695894/1‑100 (MQ=255)
                                       gtggtgATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACg                                                              <  1:2514587/100‑1 (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                                                      >  1:1837684/1‑100 (MQ=255)
                                                gCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCACTCACCCTGTAC‑GGGTGAGGGCGTAGATAGATTAATGCGCTTTTAAGGCTTTGGct                                                     >  1:2665031/1‑99 (MQ=255)
                                                   ctactGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGtt                                                  >  2:1975462/1‑100 (MQ=255)
                                                   ctactGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGtt                                                  <  2:3604567/100‑1 (MQ=255)
                                                    tactGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGttt                                                 <  1:1184713/100‑1 (MQ=255)
                                                     actGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGtttt                                                <  2:4133742/100‑1 (MQ=255)
                                                       tGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTtct                                              >  2:1729138/1‑100 (MQ=255)
                                                        gACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTtctc                                             >  2:2998270/1‑100 (MQ=255)
                                                         aCGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTctct                                            >  1:2601328/1‑100 (MQ=255)
                                                         aCGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTctct                                            <  1:3199883/100‑1 (MQ=255)
                                                           gATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCttt                                          >  2:1171823/1‑100 (MQ=255)
                                                            aTTTACGGGTCGTTGGTGTAAAGATTATTGCCCTCGCCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTa                                         <  1:4193610/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)
                                                                    ggCGTTGGTGTAAAGATTATTGCCCTTACCTTGTAT‑GGGTGAGGGTGTAGAGATATTAATGCGTTTTTACGGTTTTGGCGGTTTTCTCTTTAAACAGATa                                 <  2:1121113/100‑1 (MQ=255)
                                                                    ggCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATa                                 >  1:1064016/1‑100 (MQ=255)
                                                                    ggCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATa                                 <  2:2761787/100‑1 (MQ=255)
                                                                    ggCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCACTTTAAAAAGAtc                                 <  1:3265474/100‑2 (MQ=255)
                                                                       gTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGcc                              <  1:426853/100‑1 (MQ=255)
                                                                        ttGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGACGGCGTAGAGAGATTAATGCCCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCg                             <  1:2398880/100‑1 (MQ=255)
                                                                          ggTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGAt                           >  1:3875800/1‑100 (MQ=255)
                                                                            tgtAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATAc                         >  2:720388/1‑100 (MQ=255)
                                                                             gtAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATAcc                        >  2:4102676/1‑100 (MQ=255)
                                                                               aaaGATTATTGCCCTCACCCTGTAC‑GGGTGAGGNCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCNGTTTTNTCTNNAAACNGNTNNCCGATACCTa                      >  2:2151058/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)
                                                                                  gATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAaca                   <  2:1874562/100‑1 (MQ=255)
                                                                                   ccgatctcccTCACCC‑GTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAac                    <  1:3282151/91‑1 (MQ=255)
                                                                                    ttattGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACg                 <  2:57206/100‑1 (MQ=255)
                                                                                        tGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCa             <  2:1011246/100‑1 (MQ=255)
                                                                                        tGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCa             <  2:778741/100‑1 (MQ=255)
                                                                                        cgCCCTCACCC‑GTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCa             >  2:184035/2‑100 (MQ=39)
                                                                                        cgCCCTCACCC‑GTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGGTTTCTCTTTAAACCGATAGCCGATACTTAAAACGATc              >  1:1341688/2‑99 (MQ=37)
                                                                                         gCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCa             >  2:2563541/1‑99 (MQ=255)
                                                                                         gCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCa             <  2:3726596/99‑1 (MQ=255)
                                                                                         gCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAg            >  1:1132778/1‑100 (MQ=255)
                                                                                          cccTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGc           >  1:28694/1‑100 (MQ=255)
                                                                                          cccTCACCC‑GTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAacagatcggaa           >  2:3282167/1‑92 (MQ=255)
                                                                                          cccTCACCC‑GTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGc           >  1:3111919/1‑100 (MQ=255)
                                                                                          cccTCACCC‑GTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTGTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGc           >  1:922146/1‑100 (MQ=255)
                                                                                           ccTCACCCTGTGC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGATTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGcc          >  2:2953398/1‑100 (MQ=255)
                                                                                           ccTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGcc          <  2:1038566/100‑1 (MQ=255)
                                                                                           ccTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGcc          <  2:2345236/100‑1 (MQ=255)
                                                                                           ccTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGcc          <  2:1706101/100‑1 (MQ=255)
                                                                                            cTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCa         >  2:4315943/1‑100 (MQ=255)
                                                                                            cTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCa         >  1:634739/1‑100 (MQ=255)
                                                                                             tCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGcc          <  2:2203246/98‑1 (MQ=255)
                                                                                              cACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACCGATAGCCGATACCTAACACGATCAGCCATa       >  1:3639578/1‑100 (MQ=255)
                                                                                                cccTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATAcc     >  1:25683/1‑100 (MQ=255)
                                                                                                cccTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATAcc     >  1:1397548/1‑100 (MQ=255)
                                                                                                   tGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCggg  <  2:2293631/100‑1 (MQ=255)
                                                                                                   tGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCggg  <  2:1430324/100‑1 (MQ=255)
                                                                                                   |                                                                                                    
GCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGG  >  NC_000913/120055‑120253

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