Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,054,917 A→G 10.8% Q80R (CAG→CGG)  exbD → TonB‑gated outer membrane transporter gating inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,054,9170AG10.8% 92.4 / 3.0 37Q80R (CAG→CGG) exbDTonB‑gated outer membrane transporter gating inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (14/19);  new base G (2/2);  total (16/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.57e-01

CTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCTTGGCCCGATGCTCGACGCCAAGACCAAGGGTGACAAGGAAACCACCATCTTTTTCCAGGCTGACAAAGGCGTGGATTACGGCGACCTGATGGAAGTGATGAACAACATGCGCGCGGCCGGCTACC  >  NC_002947/6054769‑6055063
                                                                                                                                                    |                                                                                                                                                  
ctgATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCTACCAGAAGCTGTATGTCGGCGACGATCAg                                                                                                                                                   <  1:250679/150‑1 (MQ=255)
ctgATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGCCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAg                                                                                                                                                   <  2:466664/150‑1 (MQ=255)
                     gcTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCtt                                                                                                                              <  1:57407/150‑1 (MQ=255)
                            tGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCTTGGCCCGa                                                                                                                       <  2:400289/150‑1 (MQ=255)
                               cTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCTTGGCCCGATGc                                                                                                                    <  2:371062/150‑1 (MQ=255)
                                        aCATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCGGGTTGCTGCACCCGACCAGCTTGGCCCGATGCTCGACGCCa                                                                                                           >  1:510795/1‑150 (MQ=255)
                                        aCATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCGGGTTGCTGCACCCGACCAGCTTGGCCCGATGCTCGACGCCa                                                                                                           >  1:510743/1‑150 (MQ=255)
                                                            gCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGc                                                                                                                                <  2:74976/109‑1 (MQ=255)
                                                            gCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGc                                                                                                                                >  1:74976/1‑109 (MQ=255)
                                                             ccTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGCTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCTTGGCCCGATGCTCGACGCCAAGACCAAGGGTGACAAGGaaa                                                                                      <  2:386551/150‑1 (MQ=255)
                                                               tCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCTTGGCCCGATGCTCGACGCCAAGACCAAGGGTGACAAGGAAacc                                                                                    >  1:174871/1‑150 (MQ=255)
                                                                 aaCCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTAGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATAAGGTTGCTGCACCCGACCAGCTTGGCCCGATGCTCGACGCCAAGACCAAGGGTGACAAGGAAaccac                                                                                  <  2:102296/150‑1 (MQ=255)
                                                                    cGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCTTGGCCCGATGCTCGACGCCAAGACCAAGGGTGACAAGGAAACCACCAt                                                                               >  1:476988/1‑150 (MQ=255)
                                                                      ccAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCTTGGCCCGATGCTCGACGCCAAGACCAAGGGTGACAAGGAAACCACCATCt                                                                             >  1:379836/1‑150 (MQ=255)
                                                                               cgcCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCTTGGCCCGATGCTCGACGCCAAGACCAAGGGTGACAAGGAAACCACCATCTTTTTCCAgg                                                                    <  2:194431/150‑1 (MQ=255)
                                                                                  cAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCTTGGCCCGATGCTCGACGCCAAGACCAAGGGTGACAAGGAAACCACCATCTTTTTCCAGGCTg                                                                 >  1:121817/1‑150 (MQ=255)
                                                                                        cGGAGAAACCGGTGCTCGGCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCTTGGCCCGATGCTCGACGCCAAGACCAAGGGTGACAAGGAAACCACCATCTTTTTCCAGGCTGACAAAg                                                           <  2:174871/150‑1 (MQ=255)
                                                                                         ggAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCTTGGCCCGATGCTCGACGCCAAGACCAAGGGTGACAAGGAAACCACCATCTTTTTCCAGGCTGACAAAgg                                                          >  1:508160/1‑150 (MQ=255)
                                                                                         ggAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCTTGGCCCGATGCTCGACGCCAAGACCAAGGGTGACAAGGAAACCACCATCTTTTTCCAGGCTGACAAAgg                                                          >  2:145427/1‑150 (MQ=255)
                                                                                               aCCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCTTGGCCCGATGCTCGACGCCAAGACCAAGGGTGACAAGGAAACCACCATCTTTTTCCAGGCTGACAAAGGCGTGGa                                                    <  1:321583/150‑1 (MQ=255)
                                                                                                      ttCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCTTGGCCCGATGCTCGACGCCAAGACCAAGGGTGACAAGGAAac                                                                                     >  1:11191/1‑110 (MQ=255)
                                                                                                      ttCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCTTGGCCCGATGCTCGACGCCAAGACCAAGGGTGACAAGGAAac                                                                                     <  2:11191/110‑1 (MQ=255)
                                                                                                         gTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCTTGGCCCGATGCTCGACGCCAAGACCAAGGGTGACAAGGAAACCACCATCTTTTTCCAGGCTGACAAAGGCGTGGATTACGGCGAc                                          >  2:227121/1‑150 (MQ=255)
                                                                                                              cGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCTTGGCCCGATGCTCGACGCCAAGACCAAGGGTGACAAGGAAACCACCATCTTTTTCCAGGCTGACAAAGGCGTGGATTACGGCGACCTGAt                                     <  2:476988/150‑1 (MQ=255)
                                                                                                               gTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCTTGGCCCGATGCTCGACGCCAAGACCAAGGGTGACAAGGAAACCACCATCTTTTTCCAGGCTGACAAAGGCGTGGATTACGGCGACCTGATg                                    <  1:52798/150‑1 (MQ=255)
                                                                                                                    ggCGGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCTTGGCCCGATGCTCGACGCCAAGACCAAGGGTGACAAGGAAACCACCATCTTTTTCCAGGCTGACAAAGGCGTGGATTACGGCGACCTGATGGAAGt                               <  2:128950/150‑1 (MQ=255)
                                                                                                                       cGACCAGAAGCTGTATGTCGGCGACGATCGGGTTGCTGCACCCGACCAGCTTGGCCCGATGCTCGACGCCAAGACCAAGGGTGACAAGGAAACCACCATCTTTTTCCAGGCTGACAAAGGCGTGGATTACGGCGACCTGATGGAAGtgat                            <  2:510743/150‑1 (MQ=255)
                                                                                                                       cGACCAGAAGCTGTATGTCGGCGACGATCGGGTTGCTGCACCCGACCAGCTTGGCCCGATGCTCGACGCCAAGACCAAGGGTGACAAGGAAACCACCATCTTTTTCCAGGCTGACAAAGGCGTGGATTACGGCGACCTGATGGAAGtgat                            <  2:510795/150‑1 (MQ=255)
                                                                                                                        gACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCTTGGCCCGATGCTCGACGCCAAGACCAAGGGTGACAAGGAAACCACCATCTTTTTCCAGGCTGACAAAGGCGTGGATTACGGCGACCTGATGGAAGtgatg                           >  2:235540/1‑150 (MQ=255)
                                                                                                                           cAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCTTGGCCCGATGCTCGACGCCAAGACCAAg                                                                                                   >  2:138641/1‑75 (MQ=255)
                                                                                                                           cAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCTTGGCCCGATGCTCGACGCCAAGACCAAg                                                                                                   <  1:138641/75‑1 (MQ=255)
                                                                                                                           cAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCTTGGCCCGATGCTCGACGCCAAGACCAAGGGTGACAAGGAAACCACCATCTTTTGCCAGGCTGACAAAGGCGTGGATTACGGCGACCTGATGGAAGTGATGaac                        >  1:177796/1‑150 (MQ=255)
                                                                                                                                 cTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCTTGGCCCGATGCTCGACGCCAAGACCAAGGGTGACAAGGAAACCACCATCTTTTTCCAGGCTGACAAAGGCGTGGATTACGGCGACCTTATGGAAGTGATGAACCACAta                  >  2:247266/1‑149 (MQ=255)
                                                                                                                                  tGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCa                                                                                                                                  >  2:154808/1‑37 (MQ=255)
                                                                                                                                  tGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCa                                                                                                                                  <  1:154808/37‑1 (MQ=255)
                                                                                                                                                gaTCAGGTTGCTGCACCCGACCAGCTTGGCCCGATGCTCGACGCCAAGACCAAGGGTGACAAGGAAACCACCATCTTTTTCCAGGCTGACAAAGGCGTGGATTACGGCGACCTGATGGAAGTGATGAACAACATGCGCGCGGCCGGCTAc   <  1:235540/150‑1 (MQ=255)
                                                                                                                                                 aTCAGGTTGCTGCACCCGACCAGCTTGGCCCGATGCGCGACGCCAAGACCAAGGGTGACAAGGAAACCACCATCTTTTTCCAGGCTGACAAAGGCGTGGATTACGGCGACCTGATGGAAGTGATGAACAACATGCGCGCGGCCGGCTAcc  <  1:247266/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CTGATCATCTTCATGGTCGCCGCTCCCTTGGCTACGGTCGACATCAAGGTCGACCTGCCGGCCTCAACCGCCAAACCGGCGCCAAGGCCGGAGAAACCGGTGTTCGTCAGCGTCAAGGCCGACCAGAAGCTGTATGTCGGCGACGATCAGGTTGCTGCACCCGACCAGCTTGGCCCGATGCTCGACGCCAAGACCAAGGGTGACAAGGAAACCACCATCTTTTTCCAGGCTGACAAAGGCGTGGATTACGGCGACCTGATGGAAGTGATGAACAACATGCGCGCGGCCGGCTACC  >  NC_002947/6054769‑6055063

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