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

GCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCGGTCTGCTGAGCCTTCTGAGCGGCTGCCATTGCGTCGTCAGCCTTACGGTAGGCTTCG  >  NC_002947/2650066‑2650352
                                                                                                                                          |                                                                                                                                                    
gCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACTTGTGGGTCGGct                                                                                                                                           <  2:206142/150‑1 (MQ=255)
  aCCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGcttc                                                                                                                                         <  1:592461/150‑1 (MQ=255)
   ccGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAGAAGCCGACCCACAAGTGGGTCGGcttct                                                                                                                                        <  2:754842/150‑1 (MQ=255)
        ggTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACTTGCTGACTGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACTTGTGGGTCGGcttttttatt                                                                                                                                   <  2:601106/150‑8 (MQ=255)
              gAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACTTGTGGGTCGGct                                                                                                                                           >  1:283974/1‑136 (MQ=255)
              gAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACTTGTGGGTCGGct                                                                                                                                           <  2:283974/136‑1 (MQ=255)
                gCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTaa                                                                                                                                     >  1:30747/1‑140 (MQ=255)
                gCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTaa                                                                                                                                     <  2:30747/140‑1 (MQ=255)
                gCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGa                                                                                                                           <  1:519426/150‑1 (MQ=255)
                                                       aaaGGCCTGCACTAAGGGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAg                                                                                    >  1:192321/3‑150 (MQ=255)
                                                                  cTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGt                                                                         <  1:572908/150‑1 (MQ=255)
                                                                              cTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGattatt                                                                                                                   <  1:128610/96‑1 (MQ=255)
                                                                              cTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGattatt                                                                                                                   >  2:128610/1‑96 (MQ=255)
                                                                                 ttCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCCTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCgg                                                          >  1:109692/1‑150 (MQ=255)
                                                                                                          ggTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTAATGGGTAGGGGCGGGGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCGGTCTGCTGAGCCTTCTGAGCGGCTGc                                 >  2:229121/1‑150 (MQ=255)
                                                                                                                               gCCGACCCACTTGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCGGTCTGCTGAGCCTTCTGAGCGGCTGCCATTGCGTCGTCAGCCTTACg            >  1:511744/1‑150 (MQ=255)
                                                                                                                                 cGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTc                                                                                                <  1:559687/64‑1 (MQ=255)
                                                                                                                                 cGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTc                                                                                                >  2:559687/1‑64 (MQ=255)
                                                                                                                                  gACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCGGTCTGCTGAGCCTTCTGAGCGGCTGCCATTGCGTCGTCAGCCTTACGGTa         <  2:192321/150‑1 (MQ=255)
                                                                                                                                         aaGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCGGTCTGCTGAGCCTTCTGAGCGGCTGCCATTGCGTCGTCAGCCTTACGGTAGGCTTCg  >  1:765958/1‑150 (MQ=255)
                                                                                                                                          |                                                                                                                                                    
GCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCGGTCTGCTGAGCCTTCTGAGCGGCTGCCATTGCGTCGTCAGCCTTACGGTAGGCTTCG  >  NC_002947/2650066‑2650352

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: