Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,139,654 A→G intergenic (+169/+105) 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,6540AG97.6% 138.0 / ‑3.5 41intergenic (+169/+105)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 A (1/0);  new base G (17/23);  total (18/23)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.32e-01

AATAGCCGAAAGGAATAAGATGAGGGAACCGCGCCATCGGCGATTTACCGCCGGGAGCCCGGGTTGCCAGGGTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACG  >  NC_000913/2139506‑2139800
                                                                                                                                                    |                                                                                                                                                  
aaTAGCCGAAAGGAATAAGATGAGGGAACCGCGCCATCGGCGATTTACCGCCGGGAGCCCGGGTTGCCAGGGTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAGCa                                                                                                                                                  <  1:120939/151‑1 (MQ=255)
       gAAAGGAATAAGATGAGGGAACCGCGCCATCGGCGATTTACCGCCGGGAGCCCGGGTTGCCAGGGTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAGCAGCCTTCg                                                                                                                                           <  1:606027/151‑1 (MQ=255)
        aaaGGAATAAGATGAGGGAACCGCGCCATCGGCGATTTACCGCCGGGAGCCCGGGTTGCCAGGGTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAGCAGCCTTCGa                                                                                                                                          <  1:407450/151‑1 (MQ=255)
                    tGAGGGAACCGCGCCATCGGCGATTTACCGCCGGGAGCCCGGGTTGCCAGGGTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAGCAGCCTTCGATCTTACATATTc                                                                                                                              >  1:584212/1‑151 (MQ=255)
                      aGGGAACCGCGCCATCGGCGATTTACCGCCGGGAGCCCGGGTTGCCAGGGTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAGCAGCCTTCGATCTTACATATTc                                                                                                                              >  2:41379/1‑149 (MQ=255)
                       gggAACCGCGCCATCGGCGATTTACCGCCGGGAGCCCGGGTTGCCAGGGTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAGCAGCCTTCGATCTTACATATTc                                                                                                                              >  2:517723/1‑148 (MQ=255)
                            ccGCGCCATCGGCGATTTACCGCCGGGAGCCCGGGTTGCCAGGGTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAGCAGCCTTCGATCTTACATATTCAGAatat                                                                                                                       >  1:560720/1‑150 (MQ=255)
                            ccGCGCCATCGGCGATTTACCGCCGGGAGCCCGGGTTGCCAGGGTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAGCAGCCTTCGATCTTACATATTCAGAata                                                                                                                        >  2:205842/1‑149 (MQ=255)
                             cgcgCCATCGGCGATTTACCGCCGGGAGCCCGGGTTGCCAGGGTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAGCAGCCTTCGATCTTACATATTCAGAATATaa                                                                                                                     <  1:208311/151‑1 (MQ=255)
                             cgcgCCATCGGCGATTTACCGCCGGGAGCCCGGGTTGCCAGGGTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAGCAGCCTTCGATCTTACATATTCAGAATATaa                                                                                                                     <  2:570678/151‑1 (MQ=255)
                              gcgcCATCGGCGATTTACCGCCGGGAGCCCGGGTTGCCAGGGTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAGCAGCCTTCGATCTTACATATTCAGAATATAat                                                                                                                    <  1:429505/151‑1 (MQ=255)
                                 ccATCGGCGATTTACCGCCGGGAGCCCGGGTTGCCAGGGTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAGCAGCCTTCGATCTTACATATTCAGAATATAatgat                                                                                                                 <  2:585760/151‑1 (MQ=255)
                                 ccATCGGCGATTTACCGCCGGGAGCCCGGGTTGCCAGGGTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAGCAGCCTTCGATCTTACATATTCAGAATATAatga                                                                                                                  <  2:361023/150‑1 (MQ=255)
                                        cGATTTACCGCCGGGAGCCCGGGTTGCCAGGGTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAGCAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCa                                                                                                          <  1:339528/151‑1 (MQ=255)
                                              aCCGCCGGGAGCCCGGGTTGCCAGGGTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAGCAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACg                                                                                                    <  2:629213/151‑1 (MQ=255)
                                                         cccGGGTTGCCAGGGTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAGCAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGa                                                                                          >  1:121193/1‑150 (MQ=255)
                                                                tGCCAGGGTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAGCAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTa                                                                                   <  2:272183/150‑1 (MQ=255)
                                                                    aGGGTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAGCAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTAcgcg                                                                               >  1:125741/1‑150 (MQ=255)
                                                                         ggtggAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAGCAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTAt                                                                         <  1:138039/151‑1 (MQ=255)
                                                                            ggAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAGCAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATc                                                                        >  2:479069/1‑149 (MQ=255)
                                                                               gaTTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAGCAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGc                                                                    <  2:393594/150‑1 (MQ=255)
                                                                                   gAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTAc                                                                >  1:295950/1‑150 (MQ=255)
                                                                                       cACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAGCAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGAc                                                           <  2:329964/151‑1 (MQ=255)
                                                                                        aCCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAGCAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACt                                                          >  2:491320/1‑151 (MQ=255)
                                                                                         cccTGGCACGTTCACCGGCTTGAGCGTTCCCTAGTAGCGAGGAACATCAAGTGAACGGAGCAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTg                                                         >  1:555407/1‑151 (MQ=255)
                                                                                           cTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAGCAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGt                                                       >  2:432997/1‑151 (MQ=255)
                                                                                                  gTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAGCAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTgg                                                >  2:305773/1‑151 (MQ=255)
                                                                                                   ttCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAGCAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTgg                                                >  2:423465/1‑150 (MQ=255)
                                                                                                        cGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAGCAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCt                                          <  1:491330/151‑1 (MQ=255)
                                                                                                        cGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAGCAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCt                                          >  1:352447/1‑151 (MQ=255)
                                                                                                          gCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAGCAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTAc                                        >  2:119834/1‑151 (MQ=255)
                                                                                                              gAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAGCAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACAtct                                    <  2:152973/151‑1 (MQ=255)
                                                                                                               aGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAGCAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACAtctt                                   <  1:428470/151‑1 (MQ=255)
                                                                                                                gCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAGCAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACAtcttc                                  <  2:183278/151‑1 (MQ=255)
                                                                                                                 cGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAGCAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACAtcttct                                 >  1:284031/1‑151 (MQ=255)
                                                                                                                          gtagCGAGGAACATCAAGTGAACGGAGCAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAAct                        <  2:352440/151‑1 (MQ=255)
                                                                                                                              cGAGGAACATCAAGTGAACGGAGCAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAActtctt                    <  2:67106/151‑1 (MQ=255)
                                                                                                                                aGGAACATCAAGTGAACGGAGCAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTca                  <  1:279973/151‑1 (MQ=255)
                                                                                                                                    aCATCAAGTGAACGGAGCAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATc              <  2:578547/151‑1 (MQ=255)
                                                                                                                                       tCAAGTGAACGGAGCAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCtt            <  2:57033/150‑1 (MQ=255)
                                                                                                                                                cGGAGCAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACg  >  1:89148/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AATAGCCGAAAGGAATAAGATGAGGGAACCGCGCCATCGGCGATTTACCGCCGGGAGCCCGGGTTGCCAGGGTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACG  >  NC_000913/2139506‑2139800

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