Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,662,472 C→T 13.8% A71A (GCC→GCT)  PP_4969 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,662,4720CT13.8% 64.9 / 3.9 29A71A (GCC→GCT) PP_4969hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (14/11);  new base T (2/2);  total (16/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 = 4.09e-01

CTGCCCAGCCGGATGCCGGCCTGACCGGTTGCGCCGCCAAGCGCAGCGCCATCGAAAACCAGCTGAAAATTGCCCGTGACCACGGCAACAGCGACCAGGTTGGAGGGCTGGAGGAAGCCCTGCGCGGTGTGGACAACTGCACCGACGCCAGCCTGCGCAAGGAACGTGAGCAGAAGGTGCTCGATGCGCGTCACGAAGTGGCGCAGCGGGAAAAGGACTTGAAGAAGGCAGAGAAGAAAGGCGACGCCGAGAAGATCAACAAGCGCAAGGACAAGCTGGCCGAGTCGCGTAAAGAG  >  NC_002947/5662324‑5662619
                                                                                                                                                    |                                                                                                                                                   
cTGCCCAGCCGGATGCCGGCCTGACCGGTTGCGCCGCCAAGCGCAGCGCCATCGAAAACCAGCTGAAAATTGCCCGTGACCACGGCAACAGCGACCAGGTTGGAGGGCTGGAGGAAGCCCTGCGCGGTGTGGACAACTGCACCGACGCCa                                                                                                                                                    >  2:491268/1‑150 (MQ=255)
                                       agcgcagcgcCATCGAAAACCAGCTGAAAATTGCCCGTGACCACGGCAACAGCGACCAGGTTGGAGGGCTGGAGGAAGCCCTGCGCGGTGTGGACAACTGCACCGACGCCAGCCTGCGCAAGGAACGTGAGCAGAAGGTGCTCGATgcgc                                                                                                             >  2:70108/1‑150 (MQ=255)
                                           cagcgcCATCGAAAACCAGCTGAAAATTGCCCGTGACCACGGCAACAGCGACCAGGTTGGAGGGCTGGAGGAAGCCCTGCGCGGTGTGGACAACTGCACCGACGCCAGCCTGCGCAAGGAACGTGAGCAGAAGGTGCTCGATGCGCGTCa                                                                                                         <  1:491268/150‑1 (MQ=255)
                                                ccATCGAAAACCAGCTGAAAATTGCCCGTGACCACGGCAACAGCGACCAGGTTGGAGGGCTGGAGGAAGCCCTGCGCGGTGTGGACAACTGCACCGACGCCAGCCTGCGCAAGGAACGTGa                                                                                                                                 >  2:431428/1‑121 (MQ=255)
                                                ccATCGAAAACCAGCTGAAAATTGCCCGTGACCACGGCAACAGCGACCAGGTTGGAGGGCTGGAGGAAGCCCTGCGCGGTGTGGACAACTGCACCGACGCCAGCCTGCGCAAGGAACGTGa                                                                                                                                 <  1:431428/121‑1 (MQ=255)
                                                  aTCGAAAACCAGCTGAAAATTGCCCGTGACCACGGCAACAGCGACCAGGTTGGAGGGCTGGAGGAAGCCCTGCGCGGTGTGGACAACTGCACCGACGCCAGCCTGCGCAAGGAACGTGAGCAGAAGGTGCTCGATGCGCGTCACGAAGTg                                                                                                  <  1:512259/150‑1 (MQ=255)
                                                     gAAAACCAGCTGAAAATTGCCCGTGACCACGGCAACAGCGACCAGGTTGGAGGGCTGGAGGAAGCCCTGCGCGGTGTGGACAACTGCACCGACGCCAGCCTGCGCAAGGAACGTGAGCAGAAGGTGCTCGATGCGCGTCACGAAGTGgcg                                                                                               >  1:182841/1‑150 (MQ=255)
                                                               tGAAAATTGCCCGTGACCACGGCAACAGCGACCAGGTTGGAGGGCTGGAGGAAGCCCTGCGCGGTGTGGACAACTGCACCGACGCCAGCCTGCGCAAGGAACGTGAGCAGAAGGTGCTCGATGCGCGTCACGAAGTGGCGCAGCGGGaaa                                                                                     >  1:113177/1‑150 (MQ=255)
                                                                     ttGCCCGTGACCACGGCAACAGCGACCAGGTTGGAGGGCTGGAGGAAGCCCTGCGCGGTGTGGACAACTGCACCGACGCCAGCCTGCGCAAGGAACGTGAGCAGAAGGTGCTCGATGCGCGTCACGAAGTGGCGCAGCGGGAAAAGGACt                                                                               <  1:533727/150‑1 (MQ=255)
                                                                                   ggCAACAGCGACCAGGTTGGAGGGCTGGAGGAAGCCCTGCGCGGTGTGGACAACTGCACCGACGCCAGCCTGCGCAAGGAACGTGAGCAGAAGGTGCTCGATGCGCGTCACGAAGTGGCGCAGCGGGAAAAGGACTTGAAGAAGGCagag                                                                 >  2:527464/1‑150 (MQ=255)
                                                                                                  gTTGGAGGGCTGGAGGAAGCCCTGCGCGGTGTGGACAACTGCACCGACGCCAGCCTGCGCAAGGAACGTGAGCAGAAGGTGCTCGATGCGCGTCACGAAGTGGCGCAGCGGGACAAGGACTTGAAGGAGGCAGAGAAGAAAGGGGACGcc                                                  >  1:197010/1‑150 (MQ=255)
                                                                                                      gAGGGCAGGAGGAAGTCCTGCGCGGGGGGGACAACTGCACCGACGCCAGCCTGCGCAAGGAACGTGAGCAGAAGGTGCTCGATGCGCGTCACGAAGTGGCGCAGCGGGAAAAGGACTTGAAGAAGGCAGACAAGAAAGGCGACGCCgaga                                              >  2:37542/1‑150 (MQ=255)
                                                                                                                 gaAGCCCTGCGCGGTGTGGACAACTGCACCGACGCCAGCCTGCGCAAGGAACGTGAGCAGAAGGTGCTCGATGCGCGTCACGAAGTGGCGCAGCGGGAAAAGGACTTGAAGAAGGCAGAGAAGAAAg                                                          <  1:147569/127‑1 (MQ=255)
                                                                                                                 gaAGCCCTGCGCGGTGTGGACAACTGCACCGACGCCAGCCTGCGCAAGGAACGTGAGCAGAAGGTGCTCGATGCGCGTCACGAAGTGGCGCAGCGGGAAAAGGACTTGAAGAAGGCAGAGAAGAAAg                                                          >  2:147569/1‑127 (MQ=255)
                                                                                                                 gaAGCCCTGCGCGGTGTGGACAACTGCACCGACGCCAGCCTGCGCAAGGAACGTGAGCAGAAGGTGCTCGATGCGCGTCACGAAGTGGCGCAGCGGGAAAAGGACTTGAAGAAGGCAGAGAAGAAAGGCGACGCCGAGAAGATCAACAAg                                   >  1:522193/1‑150 (MQ=255)
                                                                                                                  aaGCCCTGCGCGGTGTGGACAACTGCACCGACGCCAGCCTGCGCAAGGAACGTGAGCAGAAGGTGCTCGATGCGCGTCACGAAGTGGCGCAGCGGGAAAAGGACTTGAAGAAGGCAGAGAAGAAAGGCGACGCCGAGAAGATCAACAAgc                                  <  1:70108/150‑1 (MQ=255)
                                                                                                                       cTGCGCGGTGTGGACAACTGCACCGACGTCAGCCTGCGCAAGGAACGTGAGCAGAAGGTGCTCGATGCGCGTCACGAAGTGGCGCAGCGGGAAAAGGAATTGAAGAAGGCAGAGAAGAAAGGCGACGCCGAGAAGTTCAACAAGCGCAAg                             <  2:362301/150‑1 (MQ=255)
                                                                                                                          cgcgGTGTGGACAACTGCACCGACGCCAGCCTGCGCAAGGAACGTGAGCAGAAGGTGCTCGATGCGCGTCACGAAGTGGCGCAGCGGGAAAAGGACTTGAAGAAGGCAGAGAAGAAAGGCGACGCCGAGAAGATCAACAAGCGCAAGGAc                          <  2:522193/150‑1 (MQ=255)
                                                                                                                                  ggACAACTGCACCGACGCCAGCCTGCGCAAGGAACGTGAGCAGAAGGTGCTCGATGCGCGTCACGAAGTGGCGCAGCGGGAAAAGGACTTGAAGAAGGCAGAGAAGAAAGGCGACGCCGAGAAGATCAACAAGCGCAAGGACAAGCTGGc                  <  1:527464/150‑1 (MQ=255)
                                                                                                                                      aaCTGCACCGACGCCAGCCTGCGCAAGGAACGTGAGCAGAAGGTGCTCGATGCGCGTCACGAAGTGGCGCAGCGGGAAAAGGACTTGAAGAAGGCAGAGAAGAAAGGCGACGCCGAGAAGATCAACAAGCGCAAGGACAAGCTGGCCGAg              >  2:519650/1‑150 (MQ=255)
                                                                                                                                         tGCACCGACGCCAGCCTGCGCAAGGAACGTGAGCAGAAGGTGCTCGATGCGCGTCACGAAGTGGCGCAGCg                                                                                          <  2:307751/71‑1 (MQ=255)
                                                                                                                                         tGCACCGACGCCAGCCTGCGCAAGGAACGTGAGCAGAAGGTGCTCGATGCGCGTCACGAAGTGGCGCAGCg                                                                                          >  1:307751/1‑71 (MQ=255)
                                                                                                                                            aCCGACGCTAGCCTGCGCAAGGAACGTGAGCAGAAGGTGCTCGATGCGCGTCACGAAGTGGCGCAGCGGGAAAAGGACTTGAAGAAGGCAGAGAAGAAAGGCGACGCCGAGAAGATcaacaa                                    <  1:390633/122‑1 (MQ=255)
                                                                                                                                            aCCGACGCTAGCCTGCGCAAGGAACGTGAGCAGAAGGTGCTCGATGCGCGTCACGAAGTGGCGCAGCGGGAAAAGGACTTGAAGAAGGCAGAGAAGAAAGGCGACGCCGAGAAGATcaacaa                                    >  2:390511/1‑122 (MQ=255)
                                                                                                                                            aCCGACGCTAGCCTGCGCAAGGAACGTGAGCAGAAGGTGCTCGATGCGCGTCACGAAGTGGCGCAGCGGGAAAAGGACTTGAAGAAGGCAGAGAAGAAAGGCGACGCCGAGAAGATcaacaa                                    >  2:390633/1‑122 (MQ=255)
                                                                                                                                            aCCGACGCTAGCCTGCGCAAGGAACGTGAGCAGAAGGTGCTCGATGCGCGTCACGAAGTGGCGCAGCGGGAAAAGGACTTGAAGAAGGCAGAGAAGAAAGGCGACGCCGAGAAGATcaacaa                                    <  1:390511/122‑1 (MQ=255)
                                                                                                                                              cGACGCCAGCCTGCGCAAGGAACGTGAGCAGAAGGTGCTCGATGCGCGTCACGAAGTGGCGCAGCGGGAAAAGGACTTGAAGAAGGCAGAGAAGAAAGGCGACGCCGAGAAGATCAACAAGCGCAAGGACAAGCTGGCCGAGTCGCGTaa      >  1:305909/1‑150 (MQ=255)
                                                                                                                                               gACGCCAGCCTGCGCAAGGAACGTGAGCAGAAGGTGCTCGATGCGCGTCACGAAGTGGCGCAGCGGGAAAAGGACTTGAAGAAGGCAGAGAAGAAAGGCGACGCCGAGAAGATCAACAAGCGCAAGGACAAGCTGGCCGAGTCGCGTaaa     <  2:305909/150‑1 (MQ=255)
                                                                                                                                                  gCCAGCCTGCGCAAGGAACGTGAGCAGAAGGTGCTCGATGCGCGTCACGAAGTGGCGCAGCGGGAAAAGGACTTGAAGAAGGCAGAGAAGAAAGGCGACGCCGAGAAGATCAACAAGCGCAAGGACAAGCTGGCCGAGTCGCGTAAagag  >  1:528640/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CTGCCCAGCCGGATGCCGGCCTGACCGGTTGCGCCGCCAAGCGCAGCGCCATCGAAAACCAGCTGAAAATTGCCCGTGACCACGGCAACAGCGACCAGGTTGGAGGGCTGGAGGAAGCCCTGCGCGGTGTGGACAACTGCACCGACGCCAGCCTGCGCAAGGAACGTGAGCAGAAGGTGCTCGATGCGCGTCACGAAGTGGCGCAGCGGGAAAAGGACTTGAAGAAGGCAGAGAAGAAAGGCGACGCCGAGAAGATCAACAAGCGCAAGGACAAGCTGGCCGAGTCGCGTAAAGAG  >  NC_002947/5662324‑5662619

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