Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 120,154 Δ1 bp 16.2% 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.2% 114.0 / 17.5 37intergenic (+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 (17/14);  new base . (4/2);  total (21/16)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.56e-01

CTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCAC  >  NC_000913/120006‑120300
                                                                                                                                                    |                                                                                                                                                   
ctctctGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCt                                                                                                                                                     >  3:117974/1‑149 (MQ=255)
   tctGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTa                                                                                                                                                  <  2:479833/149‑1 (MQ=255)
         gcgTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTg                                                                                                                                           >  2:259604/1‑149 (MQ=255)
          cgTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGa                                                                                                                                          <  4:117974/149‑1 (MQ=255)
           gTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAg                                                                                                                                         >  2:379208/1‑149 (MQ=255)
           gTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAg                                                                                                                                         >  1:623652/1‑149 (MQ=255)
            tGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAgg                                                                                                                                        <  1:379208/149‑1 (MQ=255)
            tGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCC‑GTACAGGGTGAgg                                                                                                                                        <  1:140630/149‑1 (MQ=255)
                                      tgtttatggCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTAc                                                                                                              >  4:584404/8‑149 (MQ=255)
                                              gCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTgg                                                                                                      >  3:297284/1‑149 (MQ=255)
                                                 gCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCgg                                                                                                   >  4:583412/1‑149 (MQ=255)
                                                 gCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCgg                                                                                                   >  2:557200/1‑149 (MQ=255)
                                                        cAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTtctc                                                                                            <  2:623652/149‑1 (MQ=255)
                                                           tGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCttt                                                                                         <  4:438960/149‑1 (MQ=255)
                                                               gaagaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTa                                                                                                                                   >  4:494757/1‑103 (MQ=255)
                                                               gaagaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTa                                                                                                                                   <  3:494757/103‑1 (MQ=255)
                                                                  gaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGa                                                                                  >  2:304266/1‑149 (MQ=255)
                                                                  gaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGa                                                                                  >  2:39237/1‑149 (MQ=255)
                                                                       ccTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGcc                                                                             <  3:321645/149‑1 (MQ=255)
                                                                        cTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCg                                                                            <  1:259604/149‑1 (MQ=255)
                                                                        cTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCg                                                                            >  4:447597/1‑149 (MQ=255)
                                                                                       ggtggtgATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATc                                                             <  3:447597/149‑1 (MQ=255)
                                                                                               ttGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATAc                                                     <  1:39237/149‑1 (MQ=255)
                                                                                                  cACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAc                                                                                                                                                 >  2:138164/1‑55 (MQ=255)
                                                                                                  cACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAc                                                                                                                                                 <  1:138164/55‑1 (MQ=255)
                                                                                                         gACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAgg                                           >  3:106981/1‑149 (MQ=255)
                                                                                                            gATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGtat                                        >  1:80948/1‑149 (MQ=255)
                                                                                                                               tAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTg                     >  1:221767/1‑149 (MQ=255)
                                                                                                                                      ctacgCCCTCACCC‑GTACAGGGTGAGGGCATAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGtcat              >  1:145569/5‑149 (MQ=255)
                                                                                                                                        ttGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGtcatca            <  4:297284/149‑1 (MQ=255)
                                                                                                                                          gCCCTCACCC‑GTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGc          >  1:163185/1‑149 (MQ=255)
                                                                                                                                            ccTCACCC‑GTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGa                                                               <  3:402192/94‑1 (MQ=255)
                                                                                                                                            ccTCACCC‑GTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGa                                                               >  4:402192/1‑94 (MQ=255)
                                                                                                                                            ccTCACCC‑GTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCAt        >  1:489039/1‑149 (MQ=255)
                                                                                                                                             cTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATa       >  2:143858/1‑149 (MQ=255)
                                                                                                                                                  ccTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGa                                                                           >  1:467368/1‑76 (MQ=255)
                                                                                                                                                  ccTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGa                                                                           <  2:467368/76‑1 (MQ=255)
                                                                                                                                                  ccTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCAc  <  1:304266/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CTCTCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAC‑GGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCAC  >  NC_000913/120006‑120300

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑