Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,650,134 T→A 18.7% intergenic (+30/+104) 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,1340TA18.7% 73.4 / 12.0 32intergenic (+30/+104)PP_2320/oprIhypothetical protein/major outer membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (13/13);  new base A (2/4);  total (15/17)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

TGGTGCGTGACGTGGTGGCTGCAGAAGATGGTATGCCGGTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTG  >  NC_002947/2649992‑2650283
                                                                                                                                              |                                                                                                                                                      
tGGTGCGTGACGTGGTGGCTGCAGAAGATGGTATGCCGGTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAg                                                                                                                                                 <  2:318632/150‑1 (MQ=255)
  gTGCGTGACGTGGTGGCTGCAGAAGATGGTATGCCGGTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGCGAACGCCGTCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGc                                                                                                                                               <  2:543357/150‑1 (MQ=255)
  gTGCGTGACGTGGTGGCTGCAGAAGATGGTATGCCGGTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGc                                                                                                                                               <  1:332704/150‑1 (MQ=255)
    gCGTGACGTGGTGGCTGCAGAAGATGGTATGCCGGTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCACTAAGGGCAGGCCt                                                                                                                                             <  2:643871/150‑1 (MQ=255)
      gTGACGTGGTGGCTGCAGAAGATGGTATGCCGGTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCACTAAGGGCAGGCCtgg                                                                                                                                           <  1:458759/150‑3 (MQ=255)
       tatcttggtggCTGCAGAAGATGGTATGCCGGTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCtttt                                                                                                                                          <  1:264669/145‑1 (MQ=255)
        gACGTGGTGGCTGCAGAAGATGGTATGCCGGTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTc                                                                                                                                         <  1:717881/150‑1 (MQ=255)
             ggtggCTGCAGAAGATGGTATGCCGGTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGGAATTTGGGCGGACCCAGGACTGCCCTTAGTGCAGGCCTTTTcagatg                                                                                                                                    >  2:264669/1‑145 (MQ=255)
                     cAGAAGATGGTATGCCGGTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGAcc                                                                                                                            >  2:123857/1‑150 (MQ=255)
                             ggTATGCCGGTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGc                                                                                                                    <  2:598052/150‑1 (MQ=255)
                             ggTATGCCGGTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGc                                                                                                                    >  1:598052/1‑150 (MQ=255)
                              gTATGCCGGTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGcc                                                                                                                   <  1:123857/150‑1 (MQ=255)
                                                        gTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGa                                                                                        >  1:148150/1‑150 (MQ=255)
                                                               aCCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACaa                                                                                 <  2:509738/150‑1 (MQ=255)
                                                                  aGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCcacttgtg                                                                              >  2:151081/1‑145 (MQ=255)
                                                                     gggCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACTTGTGGGt                                                                           <  1:151081/150‑1 (MQ=255)
                                                                            aCCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCg                                                                                                                  <  1:262149/105‑1 (MQ=255)
                                                                            aCCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCg                                                                                                                  >  2:262149/1‑105 (MQ=255)
                                                                             ccGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACctgct                                                                                                                                <  1:12612/90‑1 (MQ=255)
                                                                             ccGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACctgct                                                                                                                                >  2:12612/1‑90 (MQ=255)
                                                                              cGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTGCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGcttctt                                                                  >  1:108457/1‑150 (MQ=255)
                                                                                        gAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGt                                                        <  2:148150/150‑1 (MQ=255)
                                                                                                   ccGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGAttatta                                             >  1:79690/1‑150 (MQ=255)
                                                                                                             gTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGCTCACGAACGACTAAGAAAACCAATCCACTAGTGGGTACGCTTCTCAATAGCGTGGGATGATTATTACTTGCggctg                                   >  1:233015/1‑150 (MQ=255)
                                                                                                                   tAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGt                                                        >  1:617092/1‑123 (MQ=255)
                                                                                                                   tAATTTCGGCGTACCCACGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGt                                                        <  2:617092/123‑1 (MQ=255)
                                                                                                                                 ccAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAg               >  2:585316/1‑150 (MQ=255)
                                                                                                                                 aaaGGCCTGCACTAAGGGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCttt                              <  2:49828/133‑1 (MQ=255)
                                                                                                                                 aaaGGCCTGCACTAAGGGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCttt                              >  1:49828/3‑135 (MQ=255)
                                                                                                                                  aaGGCCTGCACTAAGGGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGc                                                                                           <  2:179779/72‑1 (MQ=37)
                                                                                                                                  aaGGCCTGCACTAAGGGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGc                                                                                           >  1:179779/2‑73 (MQ=37)
                                                                                                                                            cTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGt    >  2:359606/1‑150 (MQ=255)
                                                                                                                                              tAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTg  <  1:359606/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                      
TGGTGCGTGACGTGGTGGCTGCAGAAGATGGTATGCCGGTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTG  >  NC_002947/2649992‑2650283

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