Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,139,718 T→A 7.4% intergenic (+233/+41) yegH → / ← asmA inner membrane protein/suppressor of OmpF assembly mutants; putative outer membrane protein assembly factor; inner membrane‑anchored periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,139,7180TA 7.4% 192.6 / 5.7 54intergenic (+233/+41)yegH/asmAinner membrane protein/suppressor of OmpF assembly mutants; putative outer membrane protein assembly factor; inner membrane‑anchored periplasmic protein
Reads supporting (aligned to +/- strand):  ref base T (25/25);  new base A (2/2);  total (27/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.44e-01

TGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTG  >  NC_000913/2139581‑2139864
                                                                                                                                         |                                                                                                                                                  
tggAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTa                                                                                                                                          >  3:1001813/1‑148 (MQ=255)
   agaTTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCa                                                                                                                                       <  4:1001813/148‑1 (MQ=255)
      ttGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTAc                                                                                                                                                  >  2:242166/1‑134 (MQ=255)
      ttGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTAc                                                                                                                                                  <  1:242166/134‑1 (MQ=255)
                 tGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGc                                                                                                                                    >  1:799071/1‑137 (MQ=255)
                 tGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGc                                                                                                                                    <  2:799071/137‑1 (MQ=255)
                        ttCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCa                                                                                                                     >  3:276488/1‑145 (MQ=255)
                        ttCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCa                                                                                                                     <  4:276488/145‑1 (MQ=255)
                             cGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTaa                                                                                                             >  1:99006/1‑148 (MQ=255)
                                   gAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACa                                                                                                       <  2:99006/148‑1 (MQ=255)
                                         tCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTc                                                                                                                                             >  2:759415/1‑104 (MQ=255)
                                         tCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTc                                                                                                                                             <  1:759415/104‑1 (MQ=255)
                                           cATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATc                                                                                                                                        <  3:898280/107‑1 (MQ=255)
                                           cATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATc                                                                                                                                        >  4:898280/1‑107 (MQ=255)
                                                            tCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATAAGACGCGTAAAGCGTCGCATCAGGc                                                                                                                                    >  1:156822/1‑94 (MQ=25)
                                                            tCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATAAGACGCGTAAAGCGTCGCATCAGGc                                                                                                                                    <  2:156822/94‑1 (MQ=25)
                                                             cAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTAcg                                                                                                                                                 >  2:447918/1‑80 (MQ=255)
                                                             cAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTAcg                                                                                                                                                 <  1:447918/80‑1 (MQ=255)
                                                             cAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGTCCTGTAGGCCTgataag                                                                                                                    >  1:408594/1‑104 (MQ=38)
                                                             cAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGTCCTGTAGGCCTgataag                                                                                                                    <  2:408594/109‑6 (MQ=38)
                                                             cAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCt                                                                             <  1:45569/148‑1 (MQ=255)
                                                             cAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCt                                                                             >  2:45569/1‑148 (MQ=255)
                                                                  gAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCa                                                                                                                     >  1:706060/1‑103 (MQ=255)
                                                                  gAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCa                                                                                                                     <  2:706060/103‑1 (MQ=255)
                                                                       gAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTa                                                                                                              >  1:896970/1‑105 (MQ=255)
                                                                       gAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTa                                                                                                              <  2:896970/105‑1 (MQ=255)
                                                                        aaCAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTAcg                                                                                                                                                 >  2:87372/1‑69 (MQ=255)
                                                                        aaCAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTAcg                                                                                                                                                 <  1:87372/69‑1 (MQ=255)
                                                                          cAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCaa                                                                                          <  3:445966/122‑1 (MQ=255)
                                                                          cAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCaa                                                                                          >  4:445966/1‑122 (MQ=255)
                                                                              cTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTa                                                                                                                                          <  1:282641/70‑1 (MQ=255)
                                                                              cTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTa                                                                                                                                          >  2:282641/1‑70 (MQ=255)
                                                                              cTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACAtcttct                                                                                                 <  2:309547/111‑1 (MQ=255)
                                                                              cTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACAtcttct                                                                                                 >  1:309547/1‑111 (MQ=255)
                                                                                tCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCttt                                                                                                                                                    >  1:570553/1‑58 (MQ=255)
                                                                                tCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCttt                                                                                                                                                    <  2:570553/58‑1 (MQ=255)
                                                                                tCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTAc                                                                                                                                >  2:684479/1‑78 (MQ=255)
                                                                                tCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTAc                                                                                                                                <  1:684479/78‑1 (MQ=255)
                                                                                tCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACa                                                                                                       <  4:104624/103‑1 (MQ=255)
                                                                                tCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACa                                                                                                       >  3:104624/1‑103 (MQ=255)
                                                                                              tCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTa                                                                                                                                                   <  1:61839/45‑1 (MQ=255)
                                                                                              tCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTa                                                                                                                                                   >  2:61839/1‑45 (MQ=255)
                                                                                              tCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTaa                                                                                                             >  4:1074436/1‑83 (MQ=255)
                                                                                              tCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTaa                                                                                                             <  3:1074436/83‑1 (MQ=255)
                                                                                              tCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACa                                                                                                       >  1:134124/1‑89 (MQ=255)
                                                                                              tCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACa                                                                                                       <  2:134124/89‑1 (MQ=255)
                                                                                                      taATGATGCCGCCACGAACGCCTGATGCGACGCTTTAcg                                                                                                                                                 >  3:19988/1‑39 (MQ=255)
                                                                                                      taATGATGCCGCCACGAACGCCTGATGCGACGCTTTAcg                                                                                                                                                 <  4:19988/39‑1 (MQ=255)
                                                                                                                          ccTGATAAGACGCGTAAAGCGTCGCATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGc                     >  2:133450/1‑143 (MQ=255)
                                                                                                                          ccTGATAAGACGCGTAAAGCGTCGCATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGc                     <  1:485251/143‑1 (MQ=255)
                                                                                                                            tGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAActtct                                                                                     <  2:104768/77‑1 (MQ=255)
                                                                                                                            tGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAActtct                                                                                     >  1:104768/1‑77 (MQ=255)
                                                                                                                                  gACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCt                                                                                                                  <  1:652141/42‑1 (MQ=255)
                                                                                                                                  gACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCt                                                                                                                  >  2:652141/1‑42 (MQ=255)
                                                                                                                                   aCGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCttt                                                         >  3:752198/1‑98 (MQ=255)
                                                                                                                                   aCGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCttt                                                         <  4:752198/98‑1 (MQ=255)
                                                                                                                                        ttACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTg  >  3:1047576/1‑148 (MQ=255)
                                                                                                                                         |                                                                                                                                                  
TGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTG  >  NC_000913/2139581‑2139864

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: