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.

GGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACC  >  NC_000913/120057‑120250
                                                                                                 |                                                                                                 
   ttCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACg                                                                                             <  2:3782052/100‑1 (MQ=255)
   ttCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACg                                                                                             <  2:1347723/100‑1 (MQ=255)
   ttCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACg                                                                                             <  2:184103/100‑1 (MQ=255)
   ttCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGCTGGTGTAAAGATTATTGCCCTCACCCTGTAc                                                                                               >  1:2858847/1‑99 (MQ=255)
      aaTGGCGAAGACAACCTCGTTCGTGCTCGCTGTGGTTGTTGCGCTGCTAACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACgtgt                                                                                          <  1:3018565/100‑1 (MQ=255)
       aTGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTg                                                                                         <  2:2983599/100‑1 (MQ=255)
        tGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGGTTATTGCCCTCACCCTGTACGGGTGa                                                                                        >  2:2317591/1‑100 (MQ=255)
        tGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGa                                                                                        >  1:3761069/1‑100 (MQ=255)
        tGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGGGTAAAGATTATTGCCCTCACCCTGTACGGGTg                                                                                         >  1:1742971/1‑99 (MQ=255)
          gCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAgg                                                                                      <  2:506478/100‑1 (MQ=255)
          gCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAgg                                                                                      <  2:1491931/100‑1 (MQ=255)
           cGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCGTACAGGGTGAggg                                                                                     <  1:118960/100‑1 (MQ=255)
           cGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCGTACAGGGTGAggg                                                                                     <  1:2995452/100‑1 (MQ=255)
           cGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCGTACAGGGTGAggg                                                                                     <  1:2774984/100‑1 (MQ=255)
           cGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCGTACAGGGTGAggg                                                                                     <  1:2411077/100‑1 (MQ=255)
           cGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCGTACAGGGTGAggg                                                                                     <  2:1873545/100‑1 (MQ=255)
           cGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCGTACAGGGTGAggg                                                                                     <  2:3441777/100‑1 (MQ=255)
           cGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCGTACAGGGTGAggg                                                                                     <  2:848252/100‑1 (MQ=255)
            gaggaaAACCTCGTGCGTGGTCGAGGTGGCGATTGCACTACCCACGATTTACGGGGCGTGGGTGTAACGATTATTGCCCTCACCCGTACAGGGTGAGGGc                                                                                    <  1:1261846/97‑1 (MQ=255)
            gaagaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCGTACAGGGTGAGGGc                                                                                    <  2:3640862/100‑1 (MQ=255)
            gaagaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCGTACAGGGTGAGGGc                                                                                    <  1:3224614/100‑1 (MQ=255)
            gaagaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCGTACAGGGTGAGGGc                                                                                    <  1:3041583/100‑1 (MQ=255)
            gaagaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGGGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGc                                                                                    >  2:2436192/1‑100 (MQ=255)
             aagaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCg                                                                                   >  2:1603456/1‑100 (MQ=255)
              agaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGt                                                                                  <  2:83341/100‑1 (MQ=255)
                aaaaCCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGGAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTag                                                                                >  2:1119305/1‑100 (MQ=255)
                 aaaCCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTaga                                                                               <  1:2546678/100‑1 (MQ=255)
                   aCCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTagag                                                                              <  1:1966682/99‑1 (MQ=255)
                    ccTCGTTCGTGGTCGTGGTGGTGATTGCACTACTCACGATTTACAGGGCGTTGGTGTGAAGACTATTGACCTCACCCTGTACGGGTGAGGGCGTagagag                                                                            >  1:692990/1‑100 (MQ=255)
                     cTCGTTCGTGGGCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAgagaga                                                                           <  2:903944/100‑1 (MQ=255)
                       ctttcgtGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGTGCGTTGGTGTAAAGATTATTGCCCTCGCCCTGTACGGGTGAGGGCGTAGAGAGAtt                                                                         <  1:160694/98‑1 (MQ=255)
                       cgttcgtGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGAtt                                                                         >  1:1150862/1‑100 (MQ=255)
                       cgttcgtGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGAtt                                                                         <  2:589881/100‑1 (MQ=255)
                           tgtggtcgtggtGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATg                                                                     <  1:822649/99‑1 (MQ=255)
                           cgtggtcgtggtGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTagaga                                                                             <  1:3364938/92‑1 (MQ=255)
                              ggtcgtggtGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATgcgc                                                                  <  2:585733/100‑1 (MQ=255)
                                    ggtggtgATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTAc                                                            >  1:588398/1‑100 (MQ=255)
                                     gtggtgATTGCACTACTTACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACg                                                           <  2:3760390/100‑1 (MQ=255)
                                     gtggtgATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACg                                                           <  2:3760596/100‑1 (MQ=255)
                                      tggtgATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACgg                                                          <  1:2096951/100‑1 (MQ=255)
                                       ggtgATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTCGAGAGATTAATGCGCTTTTACGGc                                                         >  2:2600272/1‑100 (MQ=255)
                                         tgATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGCAAGATTAATGCGATTTTACGGCtt                                                       >  1:3501139/1‑100 (MQ=255)
                                            ttGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTgg                                                    <  1:1662015/100‑1 (MQ=255)
                                             ttcACTCCTGACGATTTACGGGGCTTTGGTGTAAAGATTATTGCCCGCACCCGGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTg                                                     <  1:3332601/96‑1 (MQ=255)
                                               cACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCgg                                                 <  1:1500440/100‑1 (MQ=255)
                                               cACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCgg                                                 <  2:1192273/100‑1 (MQ=255)
                                                 ctactGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCgggt                                               >  2:2202430/1‑98 (MQ=255)
                                                 ctactGACGATTTACGGGGCGTTGGGGTAAAGATTATTGCCCTCACTCTGTACGGGTGTGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGtt                                               <  1:1534191/100‑1 (MQ=255)
                                                  tactGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGttt                                              <  1:1267839/100‑1 (MQ=255)
                                                   actGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGtttt                                             >  2:1407451/1‑100 (MQ=255)
                                                    ctGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTtc                                            >  2:1348528/1‑100 (MQ=255)
                                                     tGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTtct                                           <  2:1398950/100‑1 (MQ=255)
                                                      gACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTtctc                                          <  2:2324947/100‑1 (MQ=255)
                                                      gACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTtctc                                          >  2:2360661/1‑100 (MQ=255)
                                                      gACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTtctc                                          <  1:369245/100‑1 (MQ=255)
                                                       aCGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTctct                                         <  1:2739600/100‑1 (MQ=255)
                                                       aCGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTctct                                         <  2:2852886/100‑1 (MQ=255)
                                                         gATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCttt                                       >  2:3326758/1‑100 (MQ=255)
                                                          aTTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTa                                      <  2:796210/100‑1 (MQ=255)
                                                           tttACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCttt                                       <  1:3337135/98‑1 (MQ=255)
                                                             tACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAAc                                   <  1:703502/100‑1 (MQ=255)
                                                              aCGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACa                                  <  2:1225031/100‑1 (MQ=255)
                                                               cGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAg                                 <  2:1603722/100‑1 (MQ=255)
                                                                    cgctGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGc                            >  2:2607222/4‑100 (MQ=255)
                                                                     gTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGcc                           >  2:1804378/1‑100 (MQ=255)
                                                                      ttGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCg                          <  2:3411899/100‑1 (MQ=255)
                                                                        ggTGTAAAGATTATTGTCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGAt                        <  1:2213401/100‑1 (MQ=255)
                                                                         gtgtAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATa                       <  1:299208/100‑1 (MQ=255)
                                                                           gtAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTATAGAGATTAATGCGCTTTTACGGCTTTGGCGGTGTTCTCTTTAAACAGATAGCCGATAcc                     >  1:1695679/1‑100 (MQ=255)
                                                                           gtAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATAcc                     <  1:3018740/100‑1 (MQ=255)
                                                                             aaaGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTa                   >  1:1312006/1‑100 (MQ=255)
                                                                             aaaGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTa                   >  1:3558172/1‑100 (MQ=255)
                                                                             aaaGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTa                   >  2:1990964/1‑100 (MQ=255)
                                                                             aaaGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTGGGCGGTTTTCTCTTTAAAGAGATAGCGGGTACCTa                   >  1:1803387/1‑100 (MQ=255)
                                                                                  ttattGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACg              >  1:968452/1‑100 (MQ=255)
                                                                                   tattGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGa             >  2:3597094/1‑100 (MQ=255)
                                                                                   tattGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGa             >  2:351686/1‑100 (MQ=255)
                                                                                    tacgCCCTCACCCGTACAGGGTGAGGGCGTAGAGAGATCAATGCGCTTTTACGGCTTTGGCGGGTTTGTCTTTAAACAGATAGCACATACCTAACACCAt            >  1:2334765/4‑100 (MQ=25)
                                                                                     ttGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATc           <  1:3173627/100‑1 (MQ=255)
                                                                                      tGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCa          >  1:1627454/1‑100 (MQ=255)
                                                                                      tGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCa          <  2:792331/100‑1 (MQ=255)
                                                                                      cgCCCTCACCCGTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCa          >  2:2871085/2‑100 (MQ=39)
                                                                                       gCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAg         <  1:951376/100‑1 (MQ=255)
                                                                                       gCCCTCACCCGTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAg         >  1:2381889/1‑100 (MQ=255)
                                                                                        cccTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGc        >  2:3330885/1‑100 (MQ=255)
                                                                                        cccTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGc        <  1:2237162/100‑1 (MQ=255)
                                                                                        cccTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGc        <  2:777116/100‑1 (MQ=255)
                                                                                        cccTCACCCGTACAGGGTGAGGGCGTAGAGAGCTTCATGCGCTTTTACGGCTTTGGCGGGTTTCTCTTTAAACAGATAGCCGATACCTCACACGCTCAGc        >  1:1457855/1‑100 (MQ=255)
                                                                                        cccTCACCCGTACAGGGTGAGGGCGTAGAGAGCTTAATGCGCTTTTACGGCCTTGGAGGATTTATCTTTAAACAGGTAGCCGATACCTAACCCGATCAGc        >  1:1353189/1‑100 (MQ=255)
                                                                                        cccTCACCCGTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGGCGATACCTAAAACGATCAGc        >  2:2284237/1‑100 (MQ=255)
                                                                                        cccTCACCCGTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGc        >  2:347204/1‑100 (MQ=255)
                                                                                        cccTCACCCGTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGc        >  1:1607268/1‑100 (MQ=255)
                                                                                        cccTCACCCGTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGc        >  2:2284731/1‑100 (MQ=255)
                                                                                        cccTCACCCGTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGc        >  2:2145455/1‑100 (MQ=255)
                                                                                        cccTCACCCGTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGc        >  2:1718482/1‑100 (MQ=255)
                                                                                        cccTCACCCGTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGCTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGc        >  2:3245331/1‑100 (MQ=255)
                                                                                         ccTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGcc       >  2:2295888/1‑100 (MQ=255)
                                                                                          cTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGccn      <  1:2565813/100‑2 (MQ=255)
                                                                                          cTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCa      <  1:354258/100‑1 (MQ=255)
                                                                                              cccTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATAcc  <  2:3081886/100‑1 (MQ=255)
                                                                                                 |                                                                                                 
GGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACC  >  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: