Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,650,088 C→A 13.2% intergenic (+27/+107) PP_2320 → / ← oprI hypothetical protein/major outer membrane lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,650,0880CA13.2% 98.0 / 7.4 38intergenic (+27/+107)PP_2320/oprIhypothetical protein/major outer membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (15/18);  new base A (3/2);  total (18/20)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.93e-01

ACTGGGACATGGTGCGTGACGTGGTGGCTGCAGAAGATGGTATGCCGGTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCG  >  NC_002947/2649940‑2650233
                                                                                                                                                    |                                                                                                                                                 
aCTGGGACATGGTGCGTGACGTGGTGGCTGCAGAAGATGGTATGCCGGTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGcc                                                                                                                                                   <  1:43181/149‑1 (MQ=255)
    ggACATGGTGCGTGACGTGGTGGCTGCAGAAGATGGTATGCCGGTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCt                                                                                                                                                 <  1:211006/147‑1 (MQ=255)
    ggACATGGTGCGTGACGTGGTGGCTGCAGAAGATGGTATGCCGGTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTa                                                                                                                                               <  4:257546/149‑1 (MQ=255)
      aCATGGTGCGTGACGTGGTGGCTGCAGAAGATGGTATGCCGGTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGt                                                                                                                                             >  4:362706/1‑149 (MQ=255)
      aCATGGTGCGTGACGTGGTGGCTGCAGAAGATGCTATGCCGGTAGAAATTGCCGTGCCGGTGGAAAACGAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGt                                                                                                                                             <  7:74051/149‑1 (MQ=255)
         tGGTGCGTGACGTGGTGGCTGCAGAAGATGGTATGCCGGTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCa                                                                                                                                          <  8:226398/149‑1 (MQ=255)
               gTGACGTGGTGGCTGCAGAAGATGGTATGCCGGTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCACTAAGGGCAGGCCtg                                                                                                                                    <  2:212072/149‑2 (MQ=255)
                tGACGTGGTGGCTGCAGAAGATGGTATGCCGGTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCttt                                                                                                                                   <  7:232914/149‑1 (MQ=255)
                       gtggCTGCAGAAGATGGTATGCCGGTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTAc                                                                                                                             <  3:362706/148‑1 (MQ=255)
                                          tGCCGGTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCg                                                                                                          <  6:154954/148‑1 (MQ=255)
                                               gTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTc                                                                                                     >  6:133692/1‑148 (MQ=255)
                                                   aaaTTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAgg                                                                                                <  2:159819/149‑1 (MQ=255)
                                                    aaTTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGc                                                                                               <  4:353080/149‑1 (MQ=255)
                                                            tgccgGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATaaaa                                                                                        >  3:215802/1‑148 (MQ=255)
                                                             gccgGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGACCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATaaaa                                                                                        >  6:234316/1‑147 (MQ=255)
                                                              ccgGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATaaaaaa                                                                                      <  1:221916/148‑1 (MQ=255)
                                                                    gAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGAcc                                                                               >  6:379111/1‑149 (MQ=255)
                                                                      gaaCCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGc                                                                                               >  4:238807/2‑131 (MQ=255)
                                                                      caaCCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGc                                                                                               <  3:238807/130‑1 (MQ=255)
                                                                            gggggCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCcac                                                                            >  4:134310/1‑144 (MQ=255)
                                                                            gggggCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCcac                                                                            <  3:134310/144‑1 (MQ=255)
                                                                                       cGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGcttct                                                            <  5:133692/149‑1 (MQ=255)
                                                                                       cGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGcttct                                                            <  5:379111/149‑1 (MQ=255)
                                                                                         aTGGTGACGAGCATTCCGCCCGAACTGCAGTAATGTAATTTCGGCGTACCCAGGACTGCCCTTAGTGCCGGCCTTTTCGTACCTGCTTACCGCGTTTGCCGGTTTCACAGGCAATAAAAAAGCCGACCAACAAGTGTGTCGGcttctt                                                           >  2:8610/1‑148 (MQ=255)
                                                                                                      ttCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGa                                               <  7:128055/147‑1 (MQ=255)
                                                                                                      ttCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGa                                               >  8:128055/1‑147 (MQ=255)
                                                                                                          gCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGatta                                         >  7:20455/1‑149 (MQ=255)
                                                                                                                          ggTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCggctggct                         <  8:20455/149‑1 (MQ=255)
                                                                                                                           gTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCtt                        >  4:188515/1‑149 (MQ=255)
                                                                                                                              aTTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTg                      <  4:215802/148‑1 (MQ=255)
                                                                                                                                ttCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTc                    >  7:331272/1‑148 (MQ=255)
                                                                                                                                         aaaaGGCCTGCACTAAGGGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGACGGTTTCAAAGTCAAT‑AAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGTATTATTACTTTCGGCTGGCTTTGTACAGCATACg           >  8:314698/4‑147 (MQ=255)
                                                                                                                                          ccAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACg           >  6:30684/1‑147 (MQ=255)
                                                                                                                                           aaGGCCTGCACTAAGGGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAg        >  8:295352/2‑149 (MQ=255)
                                                                                                                                           aaGGCCTGCACTAAGGGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAg        >  8:233481/2‑149 (MQ=255)
                                                                                                                                           aaGGCCTGCACTAAGGGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAAT‑AAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAgg       <  7:314698/148‑1 (MQ=255)
                                                                                                                                            aGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAgg       <  3:188515/149‑1 (MQ=255)
                                                                                                                                               ccTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAAAAATCCGTGAGGATTATTAATTGCGGATGGCTTTGTACAGCATACGAATGgcg    >  5:216900/1‑149 (MQ=255)
                                                                                                                                                  gCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGgcgcg  >  4:223001/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
ACTGGGACATGGTGCGTGACGTGGTGGCTGCAGAAGATGGTATGCCGGTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCG  >  NC_002947/2649940‑2650233

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