Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,016,205 T→A 14.1% K63M (AAG→ATG)  PP_1795 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,016,2050TA14.1% 101.3 / 6.8 41K63M (AAG→ATG) PP_1795hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (21/14);  new base A (2/3);  total (24/17)
Fisher's exact test for biased strand distribution p-value = 6.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

CGGGGGATGTGTCGATGCGGTCGATGCGCAGGTCGACGGGTTCTATCACTTCGTCGAACATATGGTATTCGGTGATCAGCACGGTGGTTACGCCACCTTTCACTCGCGCCACAATGGTGTCGCCCAGCTTGGTCAGCATGTTCTTGCTGACGCCCTCGGCGGGGAAGAAATCGATGGAGCCTTTGCCCAGTGGTGCTTGCTGCAAGGCAATGCATACGGGCTTGTCTGCCGGGGGCTGCGAGGCGTAGCGGAAGATGTACATGCCCGGTTCGAGCAGCGTAACCCGAGAGC  >  NC_002947/2016062‑2016352
                                                                                                                                               |                                                                                                                                                   
cGGGGGATGTGTCGATGCGGTCGATGCGCAGGTCGACGGGTTCTATCACTTCGTCGAACATATGGTATTCGGTGATCAGCACGGTGGTTACGCCACCTTTCACTCGCGCCACAATGGTGTCGCCCAGCTTGGTCAGCATGTTCTTGCTg                                                                                                                                                <  6:323823/149‑1 (MQ=255)
      aTGTGTCGATGCGGTCGATGCGCAGGTCGACGGGTTCTATCACTTCGTCGAACATATGGTATTCGGTGATCAGCACGGTGGTTACGCCACCTTTCACTCGCGCCACAATGGTGTCGCCCAGCTTGGTCAGCAAGAACATGCTGAccaag                                                                                                                                          <  2:173770/149‑5 (MQ=255)
         tgtCGATGCGGTCGATGCGCAGGTCGACGGGTTCTATCACTTCGTCGAACATATGGTATTCGGTGATCAGCACGGTGGTTACGCCACCTTTCGCTCGCGCCACAATGGTGTCGCCCAGCTTGGTCAGCATGTTCTTGCTGACGCCCTcg                                                                                                                                       >  8:23777/1‑149 (MQ=255)
                  ggtcgatgcgCAGGTCGACGGGTTCTATCACTTCGTCGAACATATGGTATTCGGTGATCAGCACGGTGGTTACGCCACCTTTCACTCGCGCCACAATGGTGTCGCCCAGCTTGGTCAGCATGTTCTTGCTGACGCCCTCGGCgggg                                                                                                                                 >  1:36649/1‑146 (MQ=255)
                       atgcgCAGGTCGACGGGTTCTATCACTTCGTCGAACATATGGTATTCGGTGATCAGCACGGTGGTTACGCCACCTTTCACTCGCGCCACAATGGTGTCGCCCAGCTTGGTCAGCATGTTCtt                                                                                                                                                    <  3:260769/122‑1 (MQ=255)
                       atgcgCAGGTCGACGGGTTCTATCACTTCGTCGAACATATGGTATTCGGTGATCAGCACGGTGGTTACGCCACCTTTCACTCGCGCCACAATGGTGTCGCCCAGCTTGGTCAGCATGTTCtt                                                                                                                                                    >  4:260769/1‑122 (MQ=255)
                           gcAGGTCGACGGGTTCTATCACTTCGTCGAACATATGGTATTCGGTGATCAGCACGGTGGTTACGCCACCTTTCACTCGCGCCACAATGGTGTCGCCCAGCTTGGTCAGCATGTTCTTGCTGACGCCCTCGGCGGGGAAGAAATCGATg                                                                                                                     >  8:69809/1‑149 (MQ=255)
                                   aCGGTTTCTATCACTTAGTCGAACATATGGTATTCGGTGATCAGCACGGTGGTTACGCCACCTTTCACTCTCGCCACAATGGTGTCGCCCAGCTTGGTCAGCAAGAACATGCTGAc                                                                                                                                              <  3:210234/116‑1 (MQ=255)
                                   aCGGGTTCTATCACTTCGTCGAACATATGGTATTCGGTGATCAGCACGGTGGTTACGCCACCTTTCACTCGCGCCACAATGGTGTCGCCCAGCTTGGTCAGCAAGAACATGCTGAc                                                                                                                                              >  4:210234/1‑116 (MQ=255)
                                      ggTTCTATCACTTCGTCGAACATATGGTATTCGGTGATCAGCACGGTGGTTACGCCACCTTTCACTCGCGCCACAATGGTGTCGCCCAGCTTGGTCAGCATGTTCTTGCTGACGCCCTCGGCGGGGAAGAAATCGATGGAGCCTTTGcc                                                                                                          >  6:52868/1‑149 (MQ=255)
                                          cTATCACTTCGTCGAACATATGGTATTCGGTGATCAGCACGGTGGTTACTCCACATTTCACTCGCGCCACAATGGTGTCGCCCAGCTTGGTCAGCATGTTCCTGATGACGCCCTCGCCGGGGAAGAAATCGATGGAGCCTTTGCCCAgt                                                                                                      >  5:223742/1‑149 (MQ=255)
                                            aTCACTTCGTCGAACATATGGTATTCGGTGATCAGCACGGTGGTTACGCCACCTTTCACTCGCGCCACAATGGTGTCGCCCAGCTTGGTCAGCATGTTCTTGCTGACGCCCTCGGCGGGGAAGAAATCGATGGAGCCTTTGCCCAgtgg                                                                                                    >  2:178827/1‑149 (MQ=255)
                                            aTCACTTCGTCGAACATATGGTATTCAGTGATCAGCACGGTGGTTACGCCACCTTTCACTCGCGCCACAATGGTGTCGCCCAGCTTGGTCAGCATGTTCTTGCTGACGCCCTCGGCGGGGAAGAAATCGATGGAGCCTTTGCCCAgtgg                                                                                                    >  3:278799/1‑149 (MQ=255)
                                             tCACTTCGTCGAACATATGGTATTCGGTGATCAGCACGGTGGTTACGCCACCTTTCGCTCGCGCCACAATGGTGTCGCCCAGCTTGGTCAGCATGTTCTTGCTGACGCCCTCGGCGGGGAAGAAATCGATGGAGCCTTTGCCCAgt                                                                                                      <  7:23777/146‑1 (MQ=255)
                                                cTTCGTCGAACATATGGTATTCGGTGATCAGCACGGTGGTTACGCCACCTTTCACTCGCGCCACAATGGTGTCGCCCAGCTTGGTCAGCATGTTCTTGCTGACGCCCTCGGCGGGGAAGAAATCGATGGAGCCTTTGCCCAGTGGtgct                                                                                                >  6:266894/1‑149 (MQ=255)
                                                   cgtcgAACATATGGTATTCGGTGATCAGCACGGTGGTTACGCCACCTTTCACTCGCGCCACAATGGTGTCGCCCAGCTTGGTCAGCATGTTCTTGCTGACGCCCTCGGCGGGGAAGAAATCGATGGAGCCTTTGCCCAGTGGtgcttg                                                                                              <  1:178827/148‑1 (MQ=255)
                                                                   ttCGGTGATCAGCACGGTGGTTACGCCACCTTTCACTCGCGCCACAATGGTGTCGCCCAGCTTGGTCAGCATGTTCTTGCTGACGCCCTCGGCGGGGAAGAAATCGATGGAGCCTTTGCCCAGTGGTGCTTGCTGCAAGGCAATGCATa                                                                             >  4:26658/1‑149 (MQ=255)
                                                                        tGATCAGCACGGTGGTTACGCCACCTTTCACTCGCGCCACAATGGTGTCGCCCAGCTTGGTCAGCATGTTCTTGCTGACGCCCTCGGCGGGGAAGAAATCGATGGAGCCTTTGCCCAGTGGTGCTTGCTGCAAGGCAATGCATACGGGc                                                                        <  4:122596/149‑1 (MQ=255)
                                                                        tGATCAGCACGGTGGTTACGCCACCTTTCACTCGCGCCACAATGGTGTCGCCCAGCTTGGTCAGCATGTTCTTGCTGACGCCCTCGGCGGGGAAGAAATCGATGGAGCCTTTGCCCAGTGGTGCTTGCTGCAAGGCAATGCATACGGGc                                                                        <  2:36649/149‑1 (MQ=255)
                                                                                      gtTACGCCACCTTTCACTCGCGCCACAATGGTGTCGCCCAGCTTGGTCAGCATGTTCTTGCTGACGCCCTCGGCGGGGAAGAAATCGa                                                                                                                       <  8:215331/88‑1 (MQ=255)
                                                                                      gtTACGCCACCTTTCACTCGCGCCACAATGGTGTCGCCCAGCTTGGTCAGCATGTTCTTGCTGACGCCCTCGGCGGGGAAGAAATCGa                                                                                                                       >  7:215331/1‑88 (MQ=255)
                                                                                       ttACGCCACCTTTCACTCGCGCCACAATGGTGTCGCCCAGCTTGGTCAGCATGTTCTTGCTGACGCCCTCGGCGGGGAAGAAATCGATGGAGCCTTTGCCCAGTGGTGCTTGCTGCAAGGCAATGCATACGGGCTTGTCTGCCGGGGGc                                                         >  2:40147/1‑149 (MQ=255)
                                                                                                      cTCGCGCCACAATGGTGTCGCCCAGCTTGGTCAGCATGTTCTTGCTGACGCCCTCGGCGGGGAAGAAATCGATGGAGCCTTTGCCCAGTGGTGCTTGCTGCAAGGCAATGCATACGGGCTTGTCTGCCGGGGGCTGCGAGGCGTAGCg                                           >  6:150425/1‑148 (MQ=255)
                                                                                                         gcgcCAGAATGGTGTAGCCCAGCTTGGTCAGGATGTTCTTGCTGACGCCCTCGGCGGTGAAGAAATGGATGGAGCCTTTGCCCATTGGTGCTTGCTGCAAGGCAATGCATACGTGCTTGTCTGCCGGGGGCTGCGAGTCGTAGCGGAAg                                       <  2:108487/149‑1 (MQ=255)
                                                                                                         gcgcCACAATGGTGTCGCCCAGCTTGGTCAGCATGTTCTTGCTGACGCCCTCGGCGGGGAAGAAATCGATGGAGCCTTTGCCCAGTGGTGCTTGCTGCAAGGCAATGCATACGGGCTTGTCTGCCGGGGGCTGCGAGGCGTAGCGGa                                         <  1:40147/147‑1 (MQ=255)
                                                                                                             cacaATGGTGTCGCCCAGCTTGGTCAGCATGTTCTTGCTGACGCCCTCGGCGGGGAAGAAATCGATGGAGCCTTTGCCCAGTGGTGCTTGCTGCAAGGCAATGCATACGGGCTTGTCTGCCGGGGGCTGCGAGGCGTAGCGGaa                                        <  8:11006/144‑1 (MQ=255)
                                                                                                             cacaATGGTGTCGCCCAGCTTGGTCAGCATGTTCTTGCTGACGCCCTCGGCGGGGAAGAAATCGATGGAGCCTTTGCCCAGTGGTGCTTGCTGCAAGGCAATGCATACGGGCTTGTCTGCCGGGGGCTGCGAGGCGTAGCGGaa                                        >  7:11006/1‑144 (MQ=255)
                                                                                                             cacaATGGTGTCGCCCAGCTTGGTCAGCATGTTCTTGCTGACGCCCTCGGCGGGGAAGAAATCGATGGAGCCTTTGCCCAGTGGTGCTTGCTGCAAGGCAATGCATACGGGCTTGTCTGCCGGGGGCTGCGAGGCGTAGCGGAAGATg                                    >  8:21612/1‑148 (MQ=255)
                                                                                                               caATGGTGTCGCCCAGCTTGGTCAGCATGTTCTTGCTGACGCCCTCGGCGGGGAAGAAATCGATGGAGCCTTTGCCCAGTGGTGCTTGCTGCAAGGCAATGCATACGGGCTTGTCTGCCGGGGGCTGCGAGGCGTAGCGGAAGATGTAc                                 <  5:266894/149‑1 (MQ=255)
                                                                                                                aaTGGTGTCGCCCAGCTTGGTCAGCATGTTCTTGCTGACGCCCTCGGCGGGGAAGAAATCGATGGAGCCTTTGCCCAGTGGTGCTTGCTGCAAGGCAATGCATACGGGCTTGTCTGCCGGGGGCTGCGAGGCGTAGCGGAAGATGTACa                                >  2:189801/1‑149 (MQ=255)
                                                                                                                      gtCGCCCAGCTTGGTCAGCATGTTCTTGCTGACGCCCTCGGCGGGGAAGAAATCGATGGAGCCTTTGCCCAGTGGTGCTTGCTGCAAGGCAATGCATACGGGCTTGTCTGCCGGGGGCTGCGAGGCGTAGCGTAAGATGTACATGCCCg                          >  7:260448/1‑149 (MQ=255)
                                                                                                                         gCCCAGCTTGGTCAGCATGTTCTTGCTGACGCCCTCGGCGGGGAAGAAATCGATGGAGCCTTTGCCCAGTGGTGCTTGCTGCAAGGCAATGCATACGGGCTTGTCTGCCGGGGGCTGCGAGGCGTAGCGGAAGATGTACATGCCCGGtt                       <  1:115366/149‑1 (MQ=255)
                                                                                                                         gCCCAGCTTGGTCAGCATGTTCTTGCTGACGCCCTCGGCGGGGAAGAAATCGATGGAGCCTTTGCCCAGTGGTGCTTGCTGCAAGGCAATGCATACGGGCTTGTCTGCCGGGGGCTGCGAGGCGTAGCGGAAGATGTACATGCCCGGtt                       >  2:115366/1‑149 (MQ=255)
                                                                                                                         gCCCAGCTTGGTCAGCATGTTCTTGCTGACG‑CCTCGGCGGGGAAGAAATCGATGGAGCCTTTGCCCAGTGGTGCTTGCTGCAAGGCAATGCATACGGGCTTGTCTGCCGGGGGCTGCGAGGCGTAGCGGAAGATGTACATGCCCGGTTc                      >  4:304371/1‑149 (MQ=255)
                                                                                                                            cAGCTTGGTCAGCATGTTCTTGCTGACGCCCTCGGCGGGGAAGAAATCGATGGAGCCTTTGCCCAGTGGTGCTTGCTGCAAGGCAATGCATACGGGCTTGTCTGCCGGGGGCTGCGAGGCGTAGCGGAAGATGTACATGCCCGGTTCGa                    <  7:21612/149‑1 (MQ=255)
                                                                                                                                ttGGTCAGCATGTTCTTGCTGACGCCCTCGGCGGGGAAGAAATCGATGGAGCCTTTGCCCAGTGGTGCTTGCTGCAAGGCAATGCATACGGGCTTGt                                                                    >  7:159068/1‑97 (MQ=255)
                                                                                                                                ttGGTCAGCATGTTCTTGCTGACGCCCTCGGCGGGGAAGAAATCGATGGAGCCTTTGCCCAGTGGTGCTTGCTGCAAGGCAATGCATACGGGCTTGt                                                                    <  8:159068/97‑1 (MQ=255)
                                                                                                                                   gTCAGCAAGAACATGCTGACGCCCTCGGCGGGGAAGAAATCGATGGAGCCTTTGCCCAGTGGTGCTTGCTGCAAGGCAATGCATACGGGCTTGTCTGCCGGGGGCTGCGAGGCg                                                <  4:93836/114‑1 (MQ=255)
                                                                                                                                   gTCAGCAAGAACATGCTGACGCCCTCGGCGGGGAAGAAATCGATGGAGCCTTTGCCCAGTGGTGCTTGCTGCAAGGCAATGCATACGGGCTTGTCTGCCGGGGGCTGCGAGGCg                                                >  3:93836/1‑114 (MQ=255)
                                                                                                                                    tCAGCATGTTCTTGCTGACGCCCTCGGCGGGGAAGAAATCGATGGAGCCTTTGCCCAGTGGTGCTTGCTGCAAGGCAATGCATACGGGCTTGTCTGCCGGGGGCTGCGAGGCGTAGCGGAAGATGTACATGCCCGGTTCGAGCAGCGt             >  4:308768/1‑148 (MQ=255)
                                                                                                                                          tGTTCTTGCTGACGCCCTCGGCGGGGAAGAAATCGATGGAGCCTTTGCCCAGTGGTGCTTGCTGCAAGGCAATGCATACGGGCTTGTCTGCCGGGGGCTGCGAGGCGTAGCGGAAGATGTACATGCCCGGTTCGAGCAGCGTAACCCga      <  3:308768/149‑1 (MQ=255)
                                                                                                                                           gTTCTTGCTGACGCCCTCGGCGGGGAAGAAATCGATGGAGCCTTTGCCCAGTGGTGCTTGCTGCAAGGCAATGCATACGGGCTTGTCTGCCGGGGGCTGCGAGGCGTAGCGGAAGATGTACATGCCCGGTTCGAGCAGCGTAAccc        >  6:31365/1‑146 (MQ=255)
                                                                                                                                              cTTGCTGACGCCCTCGGCGGGGAAGAAATCGATGGAGCCTTTGCCCAGTGGTGCTTGCTGCAAGGCAATGCATACGGGCTTGTCTGCCGGGGGCTGCGAGGCGTAGCGGAAGATGTACATGCCCGGTTCGAGCAGCGTAACCCGAGAGc  >  6:131498/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
CGGGGGATGTGTCGATGCGGTCGATGCGCAGGTCGACGGGTTCTATCACTTCGTCGAACATATGGTATTCGGTGATCAGCACGGTGGTTACGCCACCTTTCACTCGCGCCACAATGGTGTCGCCCAGCTTGGTCAGCATGTTCTTGCTGACGCCCTCGGCGGGGAAGAAATCGATGGAGCCTTTGCCCAGTGGTGCTTGCTGCAAGGCAATGCATACGGGCTTGTCTGCCGGGGGCTGCGAGGCGTAGCGGAAGATGTACATGCCCGGTTCGAGCAGCGTAACCCGAGAGC  >  NC_002947/2016062‑2016352

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