Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 120,154 Δ1 bp 17.7% 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.17.7% 441.9 / 95.4 158intergenic (+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 (58/70);  new base . (13/15);  total (71/87)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.58e-01

Only 100 of 158 total aligned reads displayed.

GGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACC  >  NC_000913/120057‑120250
                                                                                                 |                                                                                                 
   ttCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑g                                                                                             <  2:3782052/100‑1 (MQ=255)
   ttCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑g                                                                                             <  2:1347723/100‑1 (MQ=255)
   ttCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑g                                                                                             <  2:184103/100‑1 (MQ=255)
   ttCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGCTGGTGTAAAGATTATTGCCCTCACCCTGTAc                                                                                               >  1:2858847/1‑99 (MQ=255)
      aaTGGCGAAGACAACCTCGTTCGTGCTCGCTGTGGTTGTTGCGCTGCTAACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑gtgt                                                                                          <  1:3018565/100‑1 (MQ=255)
       aTGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTg                                                                                         <  2:2983599/100‑1 (MQ=255)
        tGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGGTTATTGCCCTCACCCTGTAC‑GGGTGa                                                                                        >  2:2317591/1‑100 (MQ=255)
        tGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGa                                                                                        >  1:3761069/1‑100 (MQ=255)
        tGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGGGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTg                                                                                         >  1:1742971/1‑99 (MQ=255)
          gCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAgg                                                                                      <  2:506478/100‑1 (MQ=255)
          gCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAgg                                                                                      <  2:1491931/100‑1 (MQ=255)
           cGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCC‑GTACAGGGTGAggg                                                                                     <  1:118960/100‑1 (MQ=255)
           cGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCC‑GTACAGGGTGAggg                                                                                     <  1:2995452/100‑1 (MQ=255)
           cGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCC‑GTACAGGGTGAggg                                                                                     <  1:2774984/100‑1 (MQ=255)
           cGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCC‑GTACAGGGTGAggg                                                                                     <  1:2411077/100‑1 (MQ=255)
           cGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCC‑GTACAGGGTGAggg                                                                                     <  2:1873545/100‑1 (MQ=255)
           cGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCC‑GTACAGGGTGAggg                                                                                     <  2:3441777/100‑1 (MQ=255)
           cGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCC‑GTACAGGGTGAggg                                                                                     <  2:848252/100‑1 (MQ=255)
            gaggaaAACCTCGTGCGTGGTCGAGGTGGCGATTGCACTACCCACGATTTACGGGGCGTGGGTGTAACGATTATTGCCCTCACCC‑GTACAGGGTGAGGGc                                                                                    <  1:1261846/97‑1 (MQ=255)
            gaagaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCC‑GTACAGGGTGAGGGc                                                                                    <  2:3640862/100‑1 (MQ=255)
            gaagaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCC‑GTACAGGGTGAGGGc                                                                                    <  1:3224614/100‑1 (MQ=255)
            gaagaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCC‑GTACAGGGTGAGGGc                                                                                    <  1:3041583/100‑1 (MQ=255)
            gaagaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGGGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGc                                                                                    >  2:2436192/1‑100 (MQ=255)
             aagaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCg                                                                                   >  2:1603456/1‑100 (MQ=255)
              agaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGt                                                                                  <  2:83341/100‑1 (MQ=255)
                aaaaCCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGGAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTag                                                                                >  2:1119305/1‑100 (MQ=255)
                 aaaCCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTaga                                                                               <  1:2546678/100‑1 (MQ=255)
                   aCCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTagag                                                                              <  1:1966682/99‑1 (MQ=255)
                    ccTCGTTCGTGGTCGTGGTGGTGATTGCACTACTCACGATTTACAGGGCGTTGGTGTGAAGACTATTGACCTCACCCTGTAC‑GGGTGAGGGCGTagagag                                                                            >  1:692990/1‑100 (MQ=255)
                     cTCGTTCGTGGGCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAgagaga                                                                           <  2:903944/100‑1 (MQ=255)
                       ctttcgtGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGTGCGTTGGTGTAAAGATTATTGCCCTCGCCCTGTAC‑GGGTGAGGGCGTAGAGAGAtt                                                                         <  1:160694/98‑1 (MQ=255)
                       cgttcgtGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGAtt                                                                         >  1:1150862/1‑100 (MQ=255)
                       cgttcgtGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGAtt                                                                         <  2:589881/100‑1 (MQ=255)
                           tgtggtcgtggtGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATg                                                                     <  1:822649/99‑1 (MQ=255)
                           cgtggtcgtggtGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTagaga                                                                             <  1:3364938/92‑1 (MQ=255)
                              ggtcgtggtGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATgcgc                                                                  <  2:585733/100‑1 (MQ=255)
                                    ggtggtgATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTAc                                                            >  1:588398/1‑100 (MQ=255)
                                     gtggtgATTGCACTACTTACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACg                                                           <  2:3760390/100‑1 (MQ=255)
                                     gtggtgATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACg                                                           <  2:3760596/100‑1 (MQ=255)
                                      tggtgATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACgg                                                          <  1:2096951/100‑1 (MQ=255)
                                       ggtgATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTCGAGAGATTAATGCGCTTTTACGGc                                                         >  2:2600272/1‑100 (MQ=255)
                                         tgATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGCAAGATTAATGCGATTTTACGGCtt                                                       >  1:3501139/1‑100 (MQ=255)
                                            ttGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTgg                                                    <  1:1662015/100‑1 (MQ=255)
                                             ttcACTCCTGACGATTTACGGGGCTTTGGTGTAAAGATTATTGCCCGCACCCGGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTg                                                     <  1:3332601/96‑1 (MQ=255)
                                               cACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCgg                                                 <  1:1500440/100‑1 (MQ=255)
                                               cACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCgg                                                 <  2:1192273/100‑1 (MQ=255)
                                                 ctactGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCgggt                                               >  2:2202430/1‑98 (MQ=255)
                                                 ctactGACGATTTACGGGGCGTTGGGGTAAAGATTATTGCCCTCACTCTGTAC‑GGGTGTGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGtt                                               <  1:1534191/100‑1 (MQ=255)
                                                  tactGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGttt                                              <  1:1267839/100‑1 (MQ=255)
                                                   actGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGtttt                                             >  2:1407451/1‑100 (MQ=255)
                                                    ctGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTtc                                            >  2:1348528/1‑100 (MQ=255)
                                                     tGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTtct                                           <  2:1398950/100‑1 (MQ=255)
                                                      gACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTtctc                                          <  2:2324947/100‑1 (MQ=255)
                                                      gACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTtctc                                          >  2:2360661/1‑100 (MQ=255)
                                                      gACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTtctc                                          <  1:369245/100‑1 (MQ=255)
                                                       aCGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTctct                                         <  1:2739600/100‑1 (MQ=255)
                                                       aCGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTctct                                         <  2:2852886/100‑1 (MQ=255)
                                                         gATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCttt                                       >  2:3326758/1‑100 (MQ=255)
                                                          aTTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTa                                      <  2:796210/100‑1 (MQ=255)
                                                           tttACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCttt                                       <  1:3337135/98‑1 (MQ=255)
                                                             tACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAAc                                   <  1:703502/100‑1 (MQ=255)
                                                              aCGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACa                                  <  2:1225031/100‑1 (MQ=255)
                                                               cGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAg                                 <  2:1603722/100‑1 (MQ=255)
                                                                    cgctGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGc                            >  2:2607222/4‑100 (MQ=255)
                                                                     gTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGcc                           >  2:1804378/1‑100 (MQ=255)
                                                                      ttGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCg                          <  2:3411899/100‑1 (MQ=255)
                                                                        ggTGTAAAGATTATTGTCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGAt                        <  1:2213401/100‑1 (MQ=255)
                                                                         gtgtAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATa                       <  1:299208/100‑1 (MQ=255)
                                                                           gtAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTATAGAGATTAATGCGCTTTTACGGCTTTGGCGGTGTTCTCTTTAAACAGATAGCCGATAcc                     >  1:1695679/1‑100 (MQ=255)
                                                                           gtAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATAcc                     <  1:3018740/100‑1 (MQ=255)
                                                                             aaaGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTa                   >  1:1312006/1‑100 (MQ=255)
                                                                             aaaGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTa                   >  1:3558172/1‑100 (MQ=255)
                                                                             aaaGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTa                   >  2:1990964/1‑100 (MQ=255)
                                                                             aaaGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTGGGCGGTTTTCTCTTTAAAGAGATAGCGGGTACCTa                   >  1:1803387/1‑100 (MQ=255)
                                                                                  ttattGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACg              >  1:968452/1‑100 (MQ=255)
                                                                                   tattGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGa             >  2:3597094/1‑100 (MQ=255)
                                                                                   tattGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGa             >  2:351686/1‑100 (MQ=255)
                                                                                    tacgCCCTCACCC‑GTACAGGGTGAGGGCGTAGAGAGATCAATGCGCTTTTACGGCTTTGGCGGGTTTGTCTTTAAACAGATAGCACATACCTAACACCAt            >  1:2334765/4‑100 (MQ=25)
                                                                                     ttGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATc           <  1:3173627/100‑1 (MQ=255)
                                                                                      tGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCa          >  1:1627454/1‑100 (MQ=255)
                                                                                      tGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCa          <  2:792331/100‑1 (MQ=255)
                                                                                      cgCCCTCACCC‑GTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCa          >  2:2871085/2‑100 (MQ=39)
                                                                                       gCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAg         <  1:951376/100‑1 (MQ=255)
                                                                                       gCCCTCACCC‑GTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAg         >  1:2381889/1‑100 (MQ=255)
                                                                                        cccTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGc        >  2:3330885/1‑100 (MQ=255)
                                                                                        cccTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGc        <  1:2237162/100‑1 (MQ=255)
                                                                                        cccTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGc        <  2:777116/100‑1 (MQ=255)
                                                                                        cccTCACCC‑GTACAGGGTGAGGGCGTAGAGAGCTTCATGCGCTTTTACGGCTTTGGCGGGTTTCTCTTTAAACAGATAGCCGATACCTCACACGCTCAGc        >  1:1457855/1‑100 (MQ=255)
                                                                                        cccTCACCC‑GTACAGGGTGAGGGCGTAGAGAGCTTAATGCGCTTTTACGGCCTTGGAGGATTTATCTTTAAACAGGTAGCCGATACCTAACCCGATCAGc        >  1:1353189/1‑100 (MQ=255)
                                                                                        cccTCACCC‑GTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGGCGATACCTAAAACGATCAGc        >  2:2284237/1‑100 (MQ=255)
                                                                                        cccTCACCC‑GTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGc        >  2:347204/1‑100 (MQ=255)
                                                                                        cccTCACCC‑GTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGc        >  1:1607268/1‑100 (MQ=255)
                                                                                        cccTCACCC‑GTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGc        >  2:2284731/1‑100 (MQ=255)
                                                                                        cccTCACCC‑GTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGc        >  2:2145455/1‑100 (MQ=255)
                                                                                        cccTCACCC‑GTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGc        >  2:1718482/1‑100 (MQ=255)
                                                                                        cccTCACCC‑GTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGCTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGc        >  2:3245331/1‑100 (MQ=255)
                                                                                         ccTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGcc       >  2:2295888/1‑100 (MQ=255)
                                                                                          cTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGccn      <  1:2565813/100‑2 (MQ=255)
                                                                                          cTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCa      <  1:354258/100‑1 (MQ=255)
                                                                                              cccTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATAcc  <  2:3081886/100‑1 (MQ=255)
                                                                                                 |                                                                                                 
GGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACC  >  NC_000913/120057‑120250

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