Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,293,270 C→A 19.2% A289A (GCC→GCA PP_1128 → OmpA family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,293,2700CA19.2% 54.6 / 5.0 26A289A (GCC→GCAPP_1128OmpA family protein
Reads supporting (aligned to +/- strand):  ref base C (10/11);  new base A (3/2);  total (13/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

CCTGAACCGCCGCGTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGTTGACAAGAAAGTAGCTGC  >  NC_002947/1293129‑1293413
                                                                                                                                             |                                                                                                                                               
ccTGAACCGCCGCGTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCtggtgg                                                                                                                                         >  1:633172/1‑150 (MQ=255)
               ggAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCaa                                                                                                                          <  2:233088/150‑1 (MQ=255)
                    tGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCGGCCTTGGCAGCtgg                                                                                                                                            >  2:81713/1‑127 (MQ=255)
                                          cATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCt                                                                                                       >  2:383014/1‑142 (MQ=255)
                                          cATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCt                                                                                                       <  1:383014/142‑1 (MQ=255)
                                           aTGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAg                                                                                              >  2:283218/1‑150 (MQ=255)
                                            tGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCCCCCCCGCCTCGATGGTGGGGGTGCAGGATCCTTACGCACCAGCGGCCAAGGCCGCTGGTGGCAAACGTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAgg                                                                                             >  2:524241/1‑150 (MQ=255)
                                             ggTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGc                                                                                            <  2:60594/150‑1 (MQ=255)
                                             ggTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGc                                                                                            <  1:592540/150‑1 (MQ=255)
                                                   gctgctGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAAc                                                                                      >  1:70553/1‑150 (MQ=255)
                                                            cAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGCCACCAGCGGCCTTGGCAGCtggtgcctt                                                                                                                                      >  2:210967/1‑89 (MQ=25)
                                                            cAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGCCACCAGCGGCCTTGGCAGCtggtgcctt                                                                                                                                      <  1:210967/93‑5 (MQ=25)
                                                              aaTACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGa                                                                           <  1:315096/150‑1 (MQ=255)
                                                                        ggCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCaa                                                                 <  2:70553/150‑1 (MQ=255)
                                                                         gccgccCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCa                                                                      >  1:485316/1‑144 (MQ=255)
                                                                         gccgccCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCa                                                                      <  2:485316/144‑1 (MQ=255)
                                                                                            cGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCa                                             >  1:156625/1‑150 (MQ=255)
                                                                                                         ggCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCaa                                <  2:156625/150‑1 (MQ=255)
                                                                                                                        cAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGa                                                                           <  1:636211/92‑1 (MQ=255)
                                                                                                                        cAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGa                                                                           >  2:636211/1‑92 (MQ=255)
                                                                                                                        cAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGTTGa                 <  1:283218/150‑1 (MQ=255)
                                                                                                                           gccACCAGCGGCCTTGGCAGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGTTGACaa              >  1:343586/3‑150 (MQ=255)
                                                                                                                                 aGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAgg                                                                                             <  1:411462/65‑1 (MQ=255)
                                                                                                                                 aGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAgg                                                                                             >  2:411462/1‑65 (MQ=255)
                                                                                                                                   cggccttggCAGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGTTGACAAGAAAGTAg      <  2:343586/143‑1 (MQ=255)
                                                                                                                                       cAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGTTGACAAGAAAGTAGCTGc  >  2:275290/1‑150 (MQ=255)
                                                                                                                                             |                                                                                                                                               
CCTGAACCGCCGCGTGGAAATGTTGCTGACCCCGCAAAACACCATGGTCGCGCTGCTGGCCAAATACCAGCAGGCCGCCCCGGCACCCGCCCCGGCCTCGATGGTGGCGGTGCAGGATGCCAAGGCACCAGCTGCCAAGGCCGCTGGTGGCAAACCTGCCGTCAAGGCTGCGGCGAAGAAACCTGTAGCCAAGGCCAAAACGCCTGCGAAGACTGCAGCCAAGAAAAAAGCTGCGCCGGCCAAGCCTGCCGCCAAGAAAGCTGCCGTTGACAAGAAAGTAGCTGC  >  NC_002947/1293129‑1293413

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