Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,257,199 C→T 20.0% Q217Q (CAG→CAA PP_4634 ← alpha/beta family hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,257,1990CT20.0% 35.6 / 4.7 20Q217Q (CAG→CAAPP_4634alpha/beta family hydrolase
Reads supporting (aligned to +/- strand):  ref base C (9/7);  new base T (2/2);  total (11/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 = 9.75e-01

CAATTTCCTCGGCCGAGATATCCACAGGTGCAATGCGCAGGTGCTTTCGCTCGGCACCCTGTAGATAACCCAACAAACGCTCGGTCGCCGCCACGGTCCCGAAAGGATCGTCAGTCAGACTGATGGCCAACGTCGTCGCACGCACCTGGGCGAAAGGTAGCGCCGCCAACGCCCGCCCGCTGGGGCGGTGCTCATAACGTGGGGTACGCGTGGTCCAGTCGTGCACCACGCCGGATGGCGTGTCTTCCAGCCACCCCAGGCGCTTGCCGGGAAAATAGCCGAATACCCGGGTC  >  NC_002947/5257054‑5257346
                                                                                                                                                 |                                                                                                                                                   
cAATTTCCTCGGCCGAGATATCCACAGGTGCAATGCGCAGGTGCTTTCGCTCGGCACCCTGTAGATAACCCAACAAACGCTCGGTCGCCGCCACGGTCCCGAAAGGATCGTCAGTCAGACTGATGGCCAACGTCGTCGCACGCACTTggg                                                                                                                                                 >  1:398130/1‑150 (MQ=255)
cAATTTCCTCGGCCGAGATATCCACAGGTGCAATGCGCAGGTGCTTTCGCTCGGCACCCTGTAGATAACCCAACAAACGCTCGGTCGCCGCCACGGTCCCGAAAGGATCGTCAGTCAGACTGATGGCCAACGTCGTCGCACGCACTTggg                                                                                                                                                 >  1:398132/1‑150 (MQ=255)
                gaTATCCACAGGTGCAATGCGCAGGTGCTTTCGCTCGGCACCCTGTAGATAACCCAACAAACGCTCGGTCGCCGCCACGGTCCCGAAAGGATCGTCAGTCAGACTGATGGCCAACGTCGTCGCACGCACCTGGGCGAAAGGTAGcgccgc                                                                                                                                 <  2:131126/150‑1 (MQ=255)
                 atatCCACAGGTGCAATGCGCAGGTGCTTTCGCTCGGCACCCTGTAGATAACCCAACAAACGCTCGGTCGCCGCCACGGTCCCGAAAGGATCGTCAGTCAGACTGATGGCCAACGTCGTCGCACGCACCTGGGCGAAAGGTAGCgccgcc                                                                                                                                >  2:280713/1‑150 (MQ=255)
                       acaGGTGCAATGCGCAGGTGCTTTCGCTCGGCACCCTGTAGATAACCCAACAAACGGTCGGTCGCCGCCACGGTCCCGAAAGGATCGTCAGTCAGACTGATGGCCAACGTCGTCGCACGCACCTGGGCGAAAGGTAGCGCCGCCAAcgcc                                                                                                                          >  2:329271/1‑150 (MQ=255)
                       acaGGTGCAATGCGCAGGTGCTTTCGCTCGGCACCCTGTAGATAACACAACAAACGCTCGGTCGCCGCCACGGTCCCGAAAGGATCGTCAGTAAGACTGATGGCCAACGTCGTCGCACGCACCTGGGCGAAAGGTAGCGCCGCCAAcgcc                                                                                                                          >  1:220492/1‑150 (MQ=255)
                        caGGTGCAATGCGCAGGTGCTTTCGCTCGGCACCCTGTAGATAACCCAACAAACGCTCGGTCGCCGCCACGGTCCCGAAAGGATCGTCAGTCAGACTGATGGCCAACGTCGTCGCACGCACTTGGGCGAAAGGTAGCGCCGCCAAcgccc                                                                                                                         <  2:398132/150‑1 (MQ=255)
                        caGGTGCAATGCGCAGGTGCTTTCGCTCGGCACCCTGTAGATAACCCAACAAACGCTCGGTCGCCGCCACGGTCCCGAAAGGATCGTCAGTCAGACTGATGGCCAACGTCGTCGCACGCACTTGGGCGAAAGGTAGCGCCGCCAAcgccc                                                                                                                         <  2:398130/150‑1 (MQ=255)
                                aTGCGCAGGTGCTTTCGCTCGGCACCCTGTAGATAACCCAACAAACGCTCGGTCGCCGCCACGGTCCCGAAAGGATCGTCAGTCAGACTGATGGCCAACGTCGTCGCACGCACCTGGGCGAAAGGTAGCGCCGCCAACGCCCGCCCGCTg                                                                                                                 <  2:229403/150‑1 (MQ=255)
                                            tttCGCTCGGCACCCTGTAGATAACCCAACAAACGCTCGGTCGCCGCCACGGTCCCGAAAGGATCGTCAGTCAGACTGATGGCCAACGTCGTCGCACGCACCTGGGCGAAAGGTAGCGCCGCCAACGCCCGCCCGCTGGGGCGGt                                                                                                          <  1:315998/145‑1 (MQ=255)
                                            tttCGCTCGGCACCCTGTAGATAACCCAACAAACGCTCGGTCGCCGCCACGGTCCCGAAAGGATCGTCAGTCAGACTGATGGCCAACGTCGTCGCACGCACCTGGGCGAAAGGTAGCGCCGCCAACGCCCGCCCGCTGGGGCGGt                                                                                                          >  2:315998/1‑145 (MQ=255)
                                              tCGCTCGGCACCCTGTAGATAACCCAACAAACGCTCGGTCGCCGCCACGGTCCCGAAAGGATCGTCAGTCAGACTGATGGCCAACGTCGTCGCACGCACCTGGGCGAAAGGTAGCGCCGCCAACGCCCGCCCGCTGGGGCGGTGCTCATa                                                                                                   <  2:37351/150‑1 (MQ=255)
                                                          cTGTAGATAACCCAACAAACGCTCGGTCGCCGCCACGGTCCCGAAAGGATCGTCAGTCAGACTGATGGCCAACGTCGTCGCACGCACCTGGGCGAAAGGTAGCGCCGCCAACGCCCGCCCGCTGGGGCGGTGCTCATAACGTGGGGTAcg                                                                                       >  1:119792/1‑150 (MQ=255)
                                                                                     cgccgcCACGGTCCCGAAAGGATCGTCAGTCAGACTGATGGCCAACGTCGTCGCACGCACCTGGGCGAAAGGTAGCGCCGCCAACGCCCGCCCGCTGGGGCGGTGCTCATAACGTGGGGTACGCGTGGTCCAGTCGTGCACCACGCCGGa                                                            <  1:280713/150‑1 (MQ=255)
                                                                                                tCCCGAAAGGATCGTCAGTCAGACTGATGGCCAACGTCGTCGCACGCACCTGGGCGAAAGGTAGCGCCGCCAACGCCCGCCCGCTGGGGCGGTGCTCATAACGTGGGGTACGCGTGGTCCAGTCGTGCACCACGCCGGATGGCGTGTCtt                                                 >  2:15129/1‑150 (MQ=255)
                                                                                                               cagtcagACTGATGGCCAACGTCGTCGCACGCACCTGGGCGAAAGGTAGCGCCGCCAACGCCCGCCCGCTGGGGCGGTGCTCATAACGTGGGGTACGCGTGGTCCAGTCGTGCACCACGCCGGATGGCGTGTCTTCCAGCCACCCCAGgc                                  <  1:329271/150‑1 (MQ=255)
                                                                                                                    agACTGATGGCCAACGTCGTCGCACGCACCTGGGCGAAAGGTAGCGCCGCCAACGCCCGCCCGCTGGGGCGGTGCTCATAACGTGGGGTACGCGTGGTCCAGTCGTGCACCACGCCGGATGGCGTGTCTTCCAGCCACCCCAGGCGCTTg                             >  2:201008/1‑150 (MQ=255)
                                                                                                                                     cgtcgCACGCACCTGGGCGAAAGGTAGCGCCGCCAACGCCCGCCCGCTGGGGCGGTGCTCATAACGTGGGGTACGCGTGGTCCAGTCGTGCACCACGCCGGATGGCGTGTCTTCCAGCCACCCCAGGCGCTTGCCGGGAAAATa                  >  2:258840/1‑144 (MQ=255)
                                                                                                                                     cgtcgCACGCACCTGGGCGAAAGGTAGCGCCGCCAACGCCCGCCCGCTGGGGCGGTGCTCATAACGTGGGGTACGCGTGGTCCAGTCGTGCACCACGCCGGATGGCGTGTCTTCCAGCCACCCCAGGCGCTTGCCGGGAAAATa                  <  1:258840/144‑1 (MQ=255)
                                                                                                                                               acCTGGGCGAAAGGTAGCGCCGCCAACGCCCGCCCGCTGGGGCGGTGCTCATAACGTGGGGTACGCGTGGTCCAGTCGTGCACCACGCCGGATGGCGTGTCTTCCAGCCACCCCAGGCGCTTGCCGGGAAAATAGCCGAATACCCGGGTc  >  2:285715/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CAATTTCCTCGGCCGAGATATCCACAGGTGCAATGCGCAGGTGCTTTCGCTCGGCACCCTGTAGATAACCCAACAAACGCTCGGTCGCCGCCACGGTCCCGAAAGGATCGTCAGTCAGACTGATGGCCAACGTCGTCGCACGCACCTGGGCGAAAGGTAGCGCCGCCAACGCCCGCCCGCTGGGGCGGTGCTCATAACGTGGGGTACGCGTGGTCCAGTCGTGCACCACGCCGGATGGCGTGTCTTCCAGCCACCCCAGGCGCTTGCCGGGAAAATAGCCGAATACCCGGGTC  >  NC_002947/5257054‑5257346

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