Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 455,991 G→T 20.0% R223R (CGG→AGG)  PP_0375 ← prolyl oligopeptidase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947455,9910GT20.0% 48.2 / 2.6 20R223R (CGG→AGG) PP_0375prolyl oligopeptidase family protein
Reads supporting (aligned to +/- strand):  ref base G (10/6);  new base T (2/2);  total (12/8)
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.99e-01

CAGGTACAGGTGCCCAATTGCCAAGGCGCAGCAGCATGGTCCGCCGGTGCAGCAGGCAACGGCTGCCATTGGCCATCGACCTCACCCCAGGGTTGCCAGAAGCCATTGCGGTCGGAAAGGCAGTAGAGCCGACCTTCGGCATCGAAGCGTGGCTGTTGCAGCGATTCTTCGGCGGCGGCTACACATTGCGCCGGTCCCCAGTGGCCGCTGGCATCGCGCACGCGGCACACCAACCGGGTGACGGTCCAGGGCTGCGCCGGGCGG  >  NC_002947/455861‑456124
                                                                                                                                  |                                                                                                                                     
cAGGTACAGGTGCCCAATTGCCAAGGCGCAGCAGCATGGTCCGCCGGTGCAGCAGGCAACGGCTGCCATTGGCCAACGACCTCACCCCAGGGTTGACAGAAGCCATTGCGGTCGGAAAGGAAGTAGATCCGACCta                                                                                                                                  >  2:497056/1‑135 (MQ=255)
            cccAATTGCCAAGGCGCAGCAGCATGGTCCGCCGGTGCAGCAGGCAACGGCTGCCATTGGCCATCGACCTCACCCCAGGGTTGCCAGAAGCCATTGCGGTCGGAAAGGCAGTAGAGCCGACCTTCGGCATCGAAGc                                                                                                                      >  1:525910/1‑136 (MQ=255)
             ccAATTGCCAAGGCGCAGCAGCATGGTCCGCCGGTGCAGCAGGCAACGGCTGCCATTGGCCATCGACCTCACCCCAGGGTTGCCAGAAGCCATTGCGGTCGGAAAGGCAGTAGAGCCGACCTTCGGCATCGAAGCg                                                                                                                     <  2:115979/136‑1 (MQ=255)
               aaTTGCCAAGGCGCAGCAGCATGGTCCGCCGGTGCAGCAGGCAACGGCTGCCATTGGCCATCGACCTCACCCCAGGGTTGCCAGAAGCCATTGCGGTCGGAAAGGCAGTAGAGCCTACCTTCGGCATCGAAGCGTg                                                                                                                   >  1:484729/1‑136 (MQ=255)
               aaTTGCCAAGGCGCAGCAGCATGGTCCGCCGGTGCAGCAGGCAACGGCTGCCATTGGCCATCGACCTCACCCCAGGGTTGCCAGAAGCCATTGCGGTCGGAAAGGCAGTAGAGCCGACCTTCGGCATCGAAGCGTg                                                                                                                   >  1:435356/1‑136 (MQ=255)
                                    tGGTCCGCCGGTGCAGCAGGCAACGGCTGCCATTGGCCATCGACCTCACCCCAGGGTTGCCAGAAGCCATTGCGGTCGGAAAGGCAGTAGAGCCGACCTTCGGCATCGAAGCGTGGCTGTTGCAGCGATTCTTcgg                                                                                              <  1:473458/136‑1 (MQ=255)
                                       tCCGCCGGTGCAGCAGGCAACGGCTGCCATTGGCCATCGACCTCACCCCAGGGTTGCCAGAAGCCATTGCGGTCGGAAAGGCAGTAGAGCCGACCTTCGGCATCGAAGCGTGGCTGTTGCAGCGATTCTTcggcgg                                                                                           >  2:459334/1‑136 (MQ=255)
                                          gccgGTGCAGCAGGCAACGGCTGCCATTGGCCATCGACCTCACCCCAGGGTTGCCAGAAGCCATTGCGGTCGGAAAGGCAGTAGAGCCGACCTTCGGCATCGAAGCGTGGCTGTTGCAGCGATTCTTcggcggcgg                                                                                        >  1:425387/1‑136 (MQ=255)
                                                                              aCCTCACCCCAGGGTTGCCAGAAGGCATTGCGGTCGGAAAGGCAGTAGAGCCGACCTTCGGCATCGAAGCGTGGCTGTTGCAGCGATTCTTCGGCGGCGGCTACACATTGCGCCGGTCCCCAGTGGCCGCTGGCAt                                                    >  2:367162/1‑136 (MQ=255)
                                                                               ccTCACCCCAGGGTTGCCAGAAGCCATTGCGGTCGGAAAGGCAGTAGAGCCTACCTTCGGCATCGAAGCGTGGCTGTTGCAGCGATTCTTCGGCGGCGGCTACACATTGCGCCGGTCCCCAGTGGCCGCTGGCATc                                                   >  2:124510/1‑136 (MQ=255)
                                                                                     cccAGGGTTGCCAGAAGCCATTGCGGTCGGAAAGGCAGTAGAGCCGACCTTCGGCATCGAAGCGTGGCTGTTGCAGCGATTCTTCGGCGGCGGCTACACATTGCGCCGGTCCCCAGTGGCCGCTGGCATCGCGCAc                                             >  1:271545/1‑136 (MQ=255)
                                                                                      ccAGGGTTGCCAGAAGCCATTGCGGTCGGAAAGGCAGTAGAGCCGACCTTCGGCATCGAAGCGTGGCTGTTGCAGCGATTCTTCGGCGGCGGCTACACATTGCGCCGGTCCCCAGTGGCCGCTGGCATCGCGCAcg                                            <  2:525910/136‑1 (MQ=255)
                                                                                                cAGAAGCCATTGCGGTCGGAAAGGCAGTAGAGCCGACCTTCGGCATCGAAGCGTGGCTGTTGCAGCGATTCTTCGGCGGCGGCTACACATTGCGCCGGTCCCCAGTGGCCGCTGGCATCGCGCACGCGGCACACCa                                  >  2:256513/1‑136 (MQ=255)
                                                                                                  gAAGCCATTGCGGTCGGAAAGGCAGTAGAGCCGACCTTCGGCATCGAAGCGTGGCTGTTGCAGCGATTCTTCGGCGGCGGCTACACATTGCGCCGGTCCCCAGTGGCCGCTGGCATCGCGCACGCGGCACACCAAc                                >  1:25112/1‑136 (MQ=255)
                                                                                                      ccATTGCGGTCGGAAAGGCAGTAGAGCCTACCTTCGGCATCGAAGCGTGGCTGTTGCAGCGATTCTTCGGCGGCGGCTACACATTGCGCCGGTCCCCAGTGGCCGCTGGCATCGCGCACGCGGCACACCAACCggg                            <  1:124510/136‑1 (MQ=255)
                                                                                                      ccATTGCGGTCGGAAAGGCAGTAGAGCCGACCTTCGGCATCGAAGCGTGGCTGTTGCAGCGATTCTTCGGCGGCGGCTACACATTGCGCCGGTCCCCAGTGGCCGCTGGCATCGCGCACGCGGCACACCAACCggg                            <  2:25112/136‑1 (MQ=255)
                                                                                                           gCGGTCGGAAAGGCAGTAGAGCCTACCTTCGGCATCGAAGCGTGGCTGTTGCAGCGATTCTTCGGCGGCGGCTACACATTGCGCCGGTCCCCAGTGGCCGCTGGCATCGCGCACGCGGCACACCAACCGGGTGACg                       <  2:484729/136‑1 (MQ=255)
                                                                                                           gCGGTCGGAAAGGCAGTAGAGCCGACCTTCGGCATCGAAGCGTGGCTGTTGCAGCGATTCTTCGGCGGCGGCTACACATTGCGCCGGTCCCCAGTGGCCGCTGGCATCGCGCACGCGGCACACCAACCGGGTGACg                       <  2:435356/136‑1 (MQ=255)
                                                                                                                          gTAGAGCCGACCTTCGGCATCGAAGCGTGGCTGTTGCAGCGATTCTTCGGCGGCGTCTACACATTGCGCCGGTCCCCAGTGGCCGCTGGCATCGCGCACGCGGCACACCAACCGGGTGACGGTCCAGGGCTGCGcc        <  1:497056/136‑1 (MQ=255)
                                                                                                                                ccGACCTTCGGCATCGAAGCGTGGCTGTTGCAGCGATTCTTCGGCGGCGGCTACACATTGCGCCGGTCCCCAGTGGCCGCTGGCATCGCGCACGCGGCACACCAACCGGGTGACGGTCCAGGGCTGCGCCGGGCgg  >  1:216801/1‑136 (MQ=255)
                                                                                                                                  |                                                                                                                                     
CAGGTACAGGTGCCCAATTGCCAAGGCGCAGCAGCATGGTCCGCCGGTGCAGCAGGCAACGGCTGCCATTGGCCATCGACCTCACCCCAGGGTTGCCAGAAGCCATTGCGGTCGGAAAGGCAGTAGAGCCGACCTTCGGCATCGAAGCGTGGCTGTTGCAGCGATTCTTCGGCGGCGGCTACACATTGCGCCGGTCCCCAGTGGCCGCTGGCATCGCGCACGCGGCACACCAACCGGGTGACGGTCCAGGGCTGCGCCGGGCGG  >  NC_002947/455861‑456124

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: