Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,425,830 T→C 36.4% C33R (TGC→CGC)  PP_3915 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,425,8300TC36.4% 3.3 / 7.1 11C33R (TGC→CGC) PP_3915hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (5/2);  new base C (2/2);  total (7/4)
Fisher's exact test for biased strand distribution p-value = 5.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCTCGCCGGCGGCATGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCTGCTGGGATGTTCGTTTCGAGAGCCATTCGAATGGGATGGATGGCGACTACACGATCGGCATTGAGATCGTCGGCCACTACATGGGCAAACCATGTGCGCGGGTGCTCGGCGAGAACTACAACGAGAACCTGCG  >  NC_002947/4425696‑4425963
                                                                                                                                      |                                                                                                                                     
gCCTCGCCGGCGGCATGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGCCCGTGAGCGACTGAACGCCCTGGAGTCAtgctg                                                                                                                                      >  1:204840/1‑136 (MQ=255)
                                                 gaaAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCTGCTGGGATGTTCGTTTCGAGAGCCATTCGAATGGGATGGATGGCGACTac                                                                                     >  1:227741/1‑136 (MQ=255)
                                                              aCCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCTGCTGGGATGTTCGTTTCGAGAGc                                                                                                                >  1:53941/1‑96 (MQ=255)
                                                              aCCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCTGCTGGGATGTTCGTTTCGAGAGc                                                                                                                <  2:53941/96‑1 (MQ=255)
                                                              aCCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCCGCTGGGATGTTCGTTTCGAGAGc                                                                                                                >  1:284346/1‑96 (MQ=255)
                                                              aCCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCCGCTGGGATGTTCGTTTCGAGAGc                                                                                                                >  1:285736/1‑96 (MQ=255)
                                                              aCCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCCGCTGGGATGTTCGTTTCGAGAGc                                                                                                                <  2:284346/96‑1 (MQ=255)
                                                              aCCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCCGCTGGGATGTTCGTTTCGAGAGc                                                                                                                <  2:285736/96‑1 (MQ=255)
                                                                                    tGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCTGCTGGGATGTTCGTTTCGAGAGCCATTCGAATGGGATGGATGGCGACTACACGATCGGCATTGAGATCGTCGGCCACTACATggg                                                  >  2:295213/1‑136 (MQ=255)
                                                                                              gcgcCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCTGCTGGGATGTTCGTTTCGAGAGCCATTCGAATGGGATGGATGGCGACTACGCGATCGGCATTGAGATCGTCGGCCACTACATGGGCAAACCAtgt                                        <  1:59706/136‑1 (MQ=255)
                                                                                              gcgcCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCTGCTGGGATGTTCGTTTCGAGAGCCATTCGAATGGGATGGATGGCGACTACACGATCGGCATTGAGATCGTCGGCCACTACATGGGCAATCCAtgt                                        >  1:83787/1‑136 (MQ=255)
                                                                                                                            ggAGTCATGCTGCTGGGATGTTCGTTTCGAGAGCCATTCGAATGGGATGGATGGCGACTACACGATCGGCATTGAGATCGTCGGCCACTACATGGGCAAACCATGTGCGCGGGTGCTCGGCGAGAACTACAACgag          >  2:50087/1‑136 (MQ=255)
                                                                                                                                      tgctgGGATGTTCGTTTCGAGAGCCATTCGAATGGGATGGATGGCGACTACACGATCGGCATTGAGATCGTCGGCCACTACATGGGCAAACCATGTGCGCGGGTGCTCGGCGAGAACTACAACGAGAACCTgcg  <  1:108249/134‑1 (MQ=255)
                                                                                                                                      tgctgGGATGTTCGTTTCGAGAGCCATTCGAATGGGATGGATGGCGACTACACGATCGGCATTGAGATCGTCGGCCACTACATGGGCAAACCATGTGCGCGGGTGCTCGGCGAGAACTACAACGAGAACCTgcg  >  2:108249/1‑134 (MQ=255)
                                                                                                                                      |                                                                                                                                     
GCCTCGCCGGCGGCATGGATTGATGAGAGGTATCAGCTGTGAGTCGGAAGAAAAGCCACTTCACCATTGTTTCCAGCACAGAGCTGGAGGAATTGCGCCGCGACCGTGAGCGACTGAACGCCCTGGAGTCATGCTGCTGGGATGTTCGTTTCGAGAGCCATTCGAATGGGATGGATGGCGACTACACGATCGGCATTGAGATCGTCGGCCACTACATGGGCAAACCATGTGCGCGGGTGCTCGGCGAGAACTACAACGAGAACCTGCG  >  NC_002947/4425696‑4425963

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