Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,449,251 A→C 46.4% S18R (AGC→CGC)  PP_5556 → bacteriophage lysis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,449,2510AC46.4% ‑5.3 / 14.6 39S18R (AGC→CGC) PP_5556bacteriophage lysis protein
Reads supporting (aligned to +/- strand):  ref base A (13/6);  new base C (15/2);  total (28/11)
Fisher's exact test for biased strand distribution p-value = 2.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

TCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCG‑CCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGG  >  NC_002947/3449121‑3449377
                                                                                                                                   |                                                                                                                              
acACCAAACGTATCTTTGGTGGTACCAACGGTCTGGAGGATCGCAAGGCCCTATACAAGCGAGCGCGCGTGGTGCTGCTGTGAGCGCTGGGCCCACCCGCCTGGGCGCGGTGGCGCTGCTGCTGGCCG‑CCAgcgcg                                                                                                                           <  2:360888/135‑2 (MQ=255)
    cAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCG‑CCCGCGCTggcg                                                                                                                       >  2:525635/1‑136 (MQ=255)
               aTGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCG‑CCACCGCTGGCGGCGCCCGGGcc                                                                                                            >  1:454836/1‑135 (MQ=255)
                        gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCG‑CCAGCGCTGGCGGCGCCCGGGCAGCCTGGgtg                                                                                                   >  1:495856/1‑136 (MQ=255)
                              gTCTGGAGGACCCCCAGGCCCTCTACCCGCGCGCGCTCGCGGTGCTGCAGGGCGCGCTTGGCCCCCCCGCGGGGGCGCGGGGGCGCTGCGGCCGGCCG‑CCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAg                                                                                             <  2:495856/136‑1 (MQ=255)
                                    aggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCG‑CCACCGCACGCGGCGCCCCGGCACCCTGGGTGTGGCAGGCCacc                                                                                       >  1:65577/1‑134 (MQ=255)
                                           cAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCG‑CCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCttaa                                                                                >  2:195787/1‑133 (MQ=255)
                                            aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCG‑CCCGCGCTGGCGGCGCCCGGGCAGCCGGGGTGTGGGAGGCCAACGCCTATaa                                                                               >  1:201997/1‑136 (MQ=255)
                                            aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCG‑CCAGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCCACGCCGATaa                                                                               >  1:363632/1‑136 (MQ=255)
                                                    taacacGCGTGCGCTCCTGGTGCTGCTGTGCGCGCTGGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCG‑CCAGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGc                                                                       <  2:221670/134‑1 (MQ=255)
                                                             agcgCTCGAGGTGCTGCAGTGAGCGCTCGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCG‑CCCGCGCTGGCGGCGCCCGGGCCGCCTGGGTGTGGGAGGCCCACGCCGATAAGGCGCAGCTCGCGGATc                                                              >  2:39130/1‑136 (MQ=255)
                                                                gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCG‑CCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCAACGCCGATAAGGCGCCGCTCGCTGATCAgg                                                           >  1:128546/1‑136 (MQ=255)
                                                                gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCG‑CCCGCGCTGGCGGCGCCCGGGCAGCCGGGGTGTGGCAGGCCCACGCCGATAAGGCGCCGCTCGCTGATCAgg                                                           >  1:131247/1‑136 (MQ=255)
                                                                gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCG‑CCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAgg                                                           >  2:143007/1‑136 (MQ=255)
                                                                gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCG‑CCACCGCTGGCGGCGCCCGGGCCGCCGGGGTGGGGCAGGCCAACGCCAAAAAGGCGCAGCTCGCTGATCAgg                                                           >  1:221670/1‑136 (MQ=255)
                                                                  tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCG‑CCCGCGCTGGCGGCGCCCGGGCCGCCTGGGTGTGGCCGGCCAACGCCTATAAGGCGCAGCTCGCTGCTCAGGcc                                                         >  1:525239/1‑136 (MQ=255)
                                                                  tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCG‑CCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCCACGCCAATAAGGCGCAGCTCGCTGATCAGGcc                                                         >  1:104695/1‑136 (MQ=255)
                                                                  tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCG‑CCCGCGCGGGCGGCGCCCGGGCAGCCGGGGTGTGGCGGGCCAACGCCGATAAGCCGCCGCGCGCTGATCAGGcc                                                         >  2:472699/1‑136 (MQ=255)
                                                                  tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCG‑CCCGCGCCGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCCCCGCCTCAAAGGAGCACCACGCGGATCAGGcc                                                         >  1:452183/1‑136 (MQ=255)
                                                                  tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCG‑CCAGCGCGGGCGGCGCCCCGGCCGCCGGGGGGTGGGAGGCCAACGCCGAAAAGGCGCCGCGCGCAGATCACGcc                                                         >  1:476209/1‑136 (MQ=255)
                                                                               gTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCG‑CCCGCGCGGGCGGCGCCCCGGCCGCCGGGGTGTGGCCGGCCCACGCCGATAAGGCGCCGCGCGCGGATCAGGCCGCGGACTGCgggc                                            >  2:291293/1‑134 (MQ=255)
                                                                                   gcgcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCG‑CCAGCGCTGGCGGCGCCCGGGCCGCCAGGGTGTGGCAGGCCAACGCCGATAAGGCGCAGCTCGCTGATCAGGCCCATGACTTCGGCAAGCa                                        >  1:294690/1‑136 (MQ=255)
                                                                                     gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCG‑CCAGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCCACGCCGATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGt                                      >  2:110232/1‑136 (MQ=255)
                                                                                     gcTGGGGCCCCCCGCCTGGTCGCGGGGGCGCGGCTGGGGGCCG‑CCCGGGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGt                                      <  2:26000/136‑1 (MQ=255)
                                                                                       ttGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCG‑CCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCAACGCCGATAAGGCGCAGCTCGCTGCTCAGGCCGCTGACTTCGGCAAGCAGTTg                                    >  1:40268/1‑136 (MQ=255)
                                                                                       ttGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCG‑CCAGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGGAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTg                                    >  1:505/1‑136 (MQ=255)
                                                                                       ttGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCG‑CCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGGAGGCCAACGCCGATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCCAGCAGTTg                                    >  2:8342/1‑136 (MQ=255)
                                                                                         gggCCCCCCGCCTGGTCGCGGTGGAGCTGCGGGTGGCCG‑CCGGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGc                                  <  2:551432/136‑1 (MQ=255)
                                                                                          ggCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCG‑CCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGGAGGCCAACGCCTATAAGGCGCAGCTCGCGGATCCGGCCGATGACTTCGGCAAGCAGTTgggg                                 >  2:125848/1‑134 (MQ=255)
                                                                                          ggCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCG‑CCAGCGCTGGCGGCGCCCGGGCAGCCGGGGTGTGGGCGGCCAACGCCGATAAGGCGCCGCTCGCGGATCAGGCCGCTGTCTTCGGCCAGCAGTTgggg                                 >  1:326323/1‑134 (MQ=255)
                                                                                                  gCCTGGTCGCCGTGGCGCTGCTGCTGGCCG‑CCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGGAGGCCCACCCCTATAAGGCGCAGCTCGCGGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAggt                         >  1:232868/1‑135 (MQ=255)
                                                                                                  gCCTGGTCGCCGTGGCGCCGCTGCTGGCCG‑CCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCCGGCCAACGCCGATAAGGCGCCGCTCGCTGATCAGGCCGCTGCCTTCGGCCAGCAGTTGGCGGAGAAggc                         >  1:496907/1‑135 (MQ=255)
                                                                                                       gtggcggTGGCGCTGCTGCTGGCCG‑CCAGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcg                    <  1:487722/130‑1 (MQ=255)
                                                                                                       gTCGCCGTGGCGCTGCTGCTGGCCC‑CCAGCCCTGGCGGCGCCCCGGCCGCCGGGGTGTGGGAGGCCCACCCCTATAAGCGGCCCGACGCGGAAAAGCCCCATGAATAAGGGAACAAGTTGGCGGAGAAGGAAcggg                    >  2:267198/1‑134 (MQ=255)
                                                                                                        gcggggTGGCGCTGGTGC‑GCCCGCCCTGCGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                   <  2:104695/131‑1 (MQ=255)
                                                                                                        gcGCCGTGGGGCTGCTGCTGGCCG‑CCAGCGGTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                   <  1:229809/135‑1 (MQ=255)
                                                                                                                       ggcgggCCG‑CCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCg    <  2:306435/132‑1 (MQ=255)
                                                                                                                        gcgggCCG‑CGCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCgg   <  1:291293/133‑1 (MQ=255)
                                                                                                                         ctgGCCG‑CCTGCGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCggg  <  2:306250/136‑1 (MQ=255)
                                                                                                                                   |                                                                                                                              
TCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCG‑CCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGG  >  NC_002947/3449121‑3449377

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑