Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,018,313 G→T 23.5% P234H (CCC→CAC)  PP_1793 ← group 2 family glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,018,3130GT23.5% 38.4 / 3.0 17P234H (CCC→CAC) PP_1793group 2 family glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base G (6/7);  new base T (2/2);  total (8/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.22e-01

CGGGCTGCCTACCAGGTGGGTGTTGTCGGCGAAGCGGATATGCACCTGCGCCTTGGGGTTTGGCACACTTACACGCAAGCCGCCATGAGTGGTGCAAAGGCCGGCAGCGCTGAGTAAGGGGTGGGGTATGCTGGCGACGGCGTCAACGTTCACGCTGTCAGCCTCGATCTTCAAGGACGCAGGAGCCTGCAGGTTGCCAAGGTCAATGGCCCATCCCTGGATCTCTCTGTGTGTCTCCAGGTACCGCACCACACCT  >  NC_002947/2018189‑2018444
                                                                                                                            |                                                                                                                                   
cGGGCTGCCTACCAGGTGGGTGTTGTCGGCGAAGCGGATATGCACCTGCGCCTTGGGGTTTGGCACACTTACACGCAAGCCGCCATGAGTGGTGCAAAGGCCGGCAGCGCTGAGTAAGGGGTGGTGTATGCTGGCg                                                                                                                          >  1:269406/1‑136 (MQ=255)
cGGGCTGCCTACCAGGTGGGTGTTGTCGGCGAAGCGGATATGCACCTGCGCCTTGGGGTTTGGCACACTTACACGCAAGCCGCCATGAGTGGTGCAAAGGCCGGCAGCGCTGAGTAAGGGGTGGTGTATGCTGGCg                                                                                                                          >  1:49496/1‑136 (MQ=255)
                       tgtCGGCGAAGCGGATATGCACCTGCGCCTTGGGGTTTGGCACACTTACACGCAAGCCGCCATGAGTGGTGCAAAGGCCGGCAGCGCTGAGTAAGGGGTGGGGTATGCTGGCGACGGCGTCAACGTTCACGCTGTc                                                                                                   >  1:136275/1‑136 (MQ=255)
                                aGCGGATATGCACCTGCGCCTTGGGGTTTGGCACACTTACACGCAAGCCGCCATGAGTGGTGCAAAGGCCGGCAGCGCTGAGTAAGGGGTGGGGTATGCTGGCGACGGCGTCAACGTTCACGCTGTCAGCCTCGAt                                                                                          >  1:242064/1‑136 (MQ=255)
                                aGCGGATATGCACCTGCGCCTTGGGGTTTGGCACACTTACACGCAAGCCGCCATGAGTGGTGCAAAGGCCGGCAGCGCTGAGTAAGGGGTGGGGTATGCTGGCGACGGCGTCAACGTTCACGCTGTCAGCCTCGAt                                                                                          >  1:25052/1‑136 (MQ=255)
                                                        ggTTTGGCACACTTACACGCAAGCCGCCATGAGTGGTGCAAAGGCCGGCAGCGCTGAGTAAGGGGTGGGGTATGCTGGCGACGGCGTCAACGTTCACGCTGTCAGCCTCGATCtt                                                                                       >  1:362841/1‑115 (MQ=255)
                                                        ggTTTGGCACACTTACACGCAAGCCGCCATGAGTGGTGCAAAGGCCGGCAGCGCTGAGTAAGGGGTGGGGTATGCTGGCGACGGCGTCAACGTTCACGCTGTCAGCCTCGATCtt                                                                                       <  2:362841/115‑1 (MQ=255)
                                                                  acTTACACGCAAGCCGCCATGAGTGGTGCAAAGGCCGGCAGCGCTGAGTAAGGGGTGGGGTATGCTGGCGACGGCGTCAACGTTCACGCTGTCAGCCTCGATCTTCAAGGACGCAGGAGCCTGCAGGTTGCCAAgg                                                        >  2:418070/1‑136 (MQ=255)
                                                                             aGCCGCCATGAGTGGTGCAAAGGCCGGCAGCGCTGAGTAAGGGGTGGGGTATGCTGGCGACGGCGTCAACGTTCACGCTGTCAGCCTCGATCTTCAAGGACGCAGGAGCCTGCAGGTTGCCAAGGTCAATGGCCCa                                             <  2:25052/136‑1 (MQ=255)
                                                                               ccgccATGAGTGGTGCAAAGGCCGGCAGCGCTGAGTAAGGGGTGGGGTATGCTGGCGACGGCGTCAACGTTCACGCTGTCAGCCTCGATCTTCAAGGACGCAGGAGCCTGCAGGTTGCCAAGGTCAATGGcc                                               >  2:373648/1‑132 (MQ=255)
                                                                               ccgccATGAGTGGTGCAAAGGCCGGCAGCGCTGAGTAAGGGGTGGGGTATGCTGGCGACGGCGTCAACGTTCACGCTGTCAGCCTCGATCTTCAAGGACGCAGGAGCCTGCAGGTTGCCAAGGTCAATGGcc                                               <  1:373648/132‑1 (MQ=255)
                                                                                   cATGAGTGGTGCAAAGGCCGGCAGCGCTGAGTAAGGGGTGGGGTATGCTGGCGACGGCGTCAACGTTCACGCTGTCAGCCTCGATCTTCAAGGACGCAGGAGCCTGCAGGTTGCCAAGGTCAATGGCCCATCCCTg                                       <  2:136275/136‑1 (MQ=255)
                                                                                   cATGAGTGGTGCAAAGGCCGGCAGCGCTGAGTAAGGGGTGGGGTATGCTGGCGACGGCGTCAACGTTCACGCTGTCAGCCACGATCTTCAAGGACGCAGGAGCCTGCAGGTTGCCAAGGTCAATGGCCCATCCCTg                                       <  2:380272/136‑1 (MQ=255)
                                                                                       aGTGGTGCAAAGGCCGGCAGCGCTGAGTAAGGGGTGGTGTATGCTGGCGACGGCGTCAACGTTCACGCTGTCAGCCTCGATCTTCAAGGACGCAGGAGCCTGCAGGTTGCCAAGGTCAATGGCCCATCCCTGGAtc                                   <  2:269406/136‑1 (MQ=255)
                                                                                       aGTGGTGCAAAGGCCGGCAGCGCTGAGTAAGGGGTGGTGTATGCTGGCGACGGCGTCAACGTTCACGCTGTCAGCCTCGATCTTCAAGGACGCAGGAGCCTGCAGGTTGCCAAGGTCAATGGCCCATCCCTGGAtc                                   <  2:49496/136‑1 (MQ=255)
                                                                                              cAAAGGCCGGCAGCGCTGAGTAAGGGGTGGGGTATGCTGGCGACGGCGTCAACGTTCACGCTGTCAGCCTCGATCTTCAAGGACGCAGGAGCCTGCAGGTTGCCAAGGTCAATGGCCCATCCCTGGATCTCTCtgt                            <  2:213372/136‑1 (MQ=255)
                                                                                                                   aaGGGGTGGGGTATGCTGGCGACGGCGTCAACGTTCACGCTGTCAGCCTCGATCTTCAAGGACGCAGGAGCCTGCAGGTTGCCAAGGTCAATGGCCCATCCCTGGATCTCTCTGTGTGTCTCCAGGTACCGcacca       <  2:339277/136‑1 (MQ=255)
                                                                                                                        gtggggtATGCTGGCGACGGCGTCAACGTTCACGCTGTCAGCCTCGATCTTCAAGGACGCAGGAGCCTGCAGGTTGCCAAGGTCAATGGCCCATCCCTGGATCTCTCTGTGTGTCTCCAGGTACCGCACCACACCt  >  2:47851/1‑136 (MQ=255)
                                                                                                                            |                                                                                                                                   
CGGGCTGCCTACCAGGTGGGTGTTGTCGGCGAAGCGGATATGCACCTGCGCCTTGGGGTTTGGCACACTTACACGCAAGCCGCCATGAGTGGTGCAAAGGCCGGCAGCGCTGAGTAAGGGGTGGGGTATGCTGGCGACGGCGTCAACGTTCACGCTGTCAGCCTCGATCTTCAAGGACGCAGGAGCCTGCAGGTTGCCAAGGTCAATGGCCCATCCCTGGATCTCTCTGTGTGTCTCCAGGTACCGCACCACACCT  >  NC_002947/2018189‑2018444

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