Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,650,204 A→T 7.8% intergenic (+100/+34) 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,2040AT 7.8% 136.9 / 2.1 50intergenic (+100/+34)PP_2320/oprIhypothetical protein/major outer membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (26/19);  new base T (2/2);  total (28/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCGGTCTGCTGAGCCTTCTGAGCGGCTGCCATTGCGTCGTCAGCCTTACGGTAGGCTTC  >  NC_002947/2650055‑2650351
                                                                                                                                                     |                                                                                                                                                   
aCCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACaa                                                                                                                                                     <  1:1332249/150‑1 (MQ=255)
  cAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGACCTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGGCCCACAAGt                                                                                                                                                   >  1:561739/1‑150 (MQ=255)
  cAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGt                                                                                                                                                   <  1:17555/150‑1 (MQ=255)
           gCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGct                                                                                                                                          <  1:802818/150‑1 (MQ=255)
            cACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGctt                                                                                                                                         >  2:797187/1‑150 (MQ=255)
                  tGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAAc                                                                                                                                   <  1:781452/150‑1 (MQ=255)
                       aCGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATcc                                                                                                                              >  1:133323/1‑150 (MQ=255)
                         gAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGTCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGt                                                                                                                            >  2:230829/1‑150 (MQ=255)
                         gAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGt                                                                                                                            >  1:781823/1‑150 (MQ=255)
                         gAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGt                                                                                                                            <  1:735095/150‑1 (MQ=255)
                           gCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGa                                                                                                                          <  1:40970/150‑1 (MQ=255)
                                  gCCCGAACTGCAGGAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGCCCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGattat                                                                                                                   >  2:767569/1‑150 (MQ=255)
                                   cccGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGattatt                                                                                                                  <  2:133323/150‑1 (MQ=255)
                                    ccGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTaa                                                                                                                                    <  1:1044469/131‑1 (MQ=255)
                                    ccGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTaa                                                                                                                                    >  2:1044469/1‑131 (MQ=255)
                                    ccGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGAttatta                                                                                                                 <  2:781823/150‑1 (MQ=255)
                                          tGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCg                                                                                                           >  2:1280690/1‑150 (MQ=255)
                                              gTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGa                                                                                                                          <  1:1137339/131‑1 (MQ=255)
                                              gTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGa                                                                                                                          >  2:1137339/1‑131 (MQ=255)
                                               tAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCggctgg                                                                                                      <  1:1280690/150‑1 (MQ=255)
                                                   gTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACTTGTGGGTCGGctttttta                                                                                                                                    <  1:1319814/116‑3 (MQ=39)
                                                   gTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACTTGTGGGTCGGctttttta                                                                                                                                    >  2:1319814/1‑111 (MQ=255)
                                                     aaTTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACCTGCGGCTGGCTTTGt                                                                                                <  2:561739/150‑1 (MQ=255)
                                                                   aaGGCCTGCACTAAGGGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTc                                                                                               >  2:152518/2‑137 (MQ=255)
                                                                   aaGGCCTGCACTAAGGGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTc                                                                                               <  1:152518/136‑1 (MQ=255)
                                                                           cccTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGt                                                                                                                            >  2:253212/1‑100 (MQ=255)
                                                                           cccTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGt                                                                                                                            <  1:253212/100‑1 (MQ=255)
                                                                             cTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGt                                                                        >  1:1248732/1‑150 (MQ=255)
                                                                              ttAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGtt                                                                       >  1:614041/1‑150 (MQ=255)
                                                                                    cAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCTTTGGCTTc                                                                 <  2:945372/150‑1 (MQ=255)
                                                                                     aGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACTTGTGGGTc                                                                                                                                              >  1:1410398/1‑72 (MQ=39)
                                                                                     aGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACTTGTGGGTc                                                                                                                                              <  2:1410398/72‑1 (MQ=39)
                                                                                                    tgctgACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCGGTCTGCTGa                                                 <  1:230829/150‑1 (MQ=255)
                                                                                                    tgctgACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTAGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCGGTCTGCTGa                                                 >  2:1151870/1‑150 (MQ=255)
                                                                                                         aCCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGtt                                                                       >  1:486286/1‑123 (MQ=255)
                                                                                                         aCCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAACCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGtt                                                                       <  2:486286/123‑1 (MQ=255)
                                                                                                               tttGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCGGTCTGCTGAGCCTTCTGAGc                                      >  2:869508/1‑150 (MQ=255)
                                                                                                                    cGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCGGTCTGCTGAGCCTTCTGAGCGGCTg                                 >  2:235913/1‑150 (MQ=255)
                                                                                                                          cAAAAGCAACAAAAAAGCCGACCCAAACGCGGGTCGGCTTCTTAACACTCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGCCTTCGTCAGCGGTCTGCTGAGCCTTCTGAGCGGCTGCCATTg                           <  2:227656/150‑1 (MQ=255)
                                                                                                                           aaaGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCGGTCTGCTGAGCCTTCTGAGCGGCTGCCATTGc                          >  1:227656/1‑150 (MQ=255)
                                                                                                                            aaGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGGCCAGCATACGCAGGGCGCGCTCGTTGGATTGGTCAACAGTCTGCTGGGCCTTCTGAGCGGCTGCCATTGcg                         >  1:1027068/1‑150 (MQ=255)
                                                                                                                                   tAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCGGTCTGCTTAGCCTTCTGAGCGGCTGCCATTGCGTCGTCAg                  >  2:145956/1‑150 (MQ=255)
                                                                                                                                      aaaaGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCGGTCTGCTGAGCCTTCTGAGCGGCTGCCATTGCGTCGTCAGCCt               >  1:421942/1‑150 (MQ=255)
                                                                                                                                             gACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCg                                                                         <  1:787042/85‑1 (MQ=255)
                                                                                                                                             gACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCg                                                                         >  2:787042/1‑85 (MQ=255)
                                                                                                                                             gACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCGGTCTGCTGAGCCTTCTGAGCGGCTGCCATTGCGTCGTCAGCCTTACGGTa        >  2:1061841/1‑150 (MQ=255)
                                                                                                                                                ccACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCGGTCTGCTGAGCCTTCTGAGCGGCTGCCATTGCGTCGTCAGCCTTACGGTAGGc     >  1:196589/1‑150 (MQ=255)
                                                                                                                                                ccACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTAACTGAGGCTGCCAATGTGCAGCATACGTAGCGCTTGCGCGTTGTCGTCGTCAGCGGTCTGCTGAGCCTTCTGAGCGGCTGCCATTGCGTCGTCAGCCTTACGGTAGGc     >  1:1242781/1‑150 (MQ=255)
                                                                                                                                                 cacaAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCGGTCTGCTGAGCCTTCTGAGCGGCTGCCATTGCGTCGTCAGCCTTACGGTAGGCt    >  2:1298465/1‑150 (MQ=255)
                                                                                                                                                   caAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCGGTCTGCTGAGCCTTCTGAGCGGCTGCCATTGCGTCGTCAGCCTTACGGTAGGCTTc  <  1:797187/150‑1 (MQ=255)
                                                                                                                                                   caAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCGGTCTGCTGAGCCTTCTGAACGGCTGCCATTGCGTCGTCAGCCTTACGGTAGGCTTc  <  1:1151870/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
ACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCGGTCTGCTGAGCCTTCTGAGCGGCTGCCATTGCGTCGTCAGCCTTACGGTAGGCTTC  >  NC_002947/2650055‑2650351

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: ‑