Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,633,616 G→T 36.0% W12L (TGG→TTG)  PP_2304 → peptidyl‑prolyl cis‑trans isomerase D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,633,6160GT36.0% 6.6 / 7.2 14W12L (TGG→TTG) PP_2304peptidyl‑prolyl cis‑trans isomerase D
Reads supporting (aligned to +/- strand):  ref base G (4/5);  new base T (2/3);  total (6/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.28e-01

ACACGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTGGATTGCCAAGACCATCATCGGCCTTATCGTCGTGTTGATGGCATTGACCGGTTTCGAAGCCATTTTCCAGGCCGCCACTCATAGCCAGGACGCCGCCAAGGTGAACGGC  >  NC_002947/2633473‑2633725
                                                                                                                                               |                                                                                                             
acacGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTGGATTGc                                                                                                         <  1:338113/150‑1 (MQ=255)
     tttGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTGGATTGCCAAGa                                                                                                    <  1:294495/150‑1 (MQ=255)
             ggaaaGACGCATCGCCGGCTGCGCCTTTCTTTTATAAGCATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTGGGATTGCCAAGCCcctcagc                                                                                            >  2:64442/3‑148 (MQ=255)
               aaaGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTGGATTGCCAAGACCATCATCgg                                                                                          >  2:300165/1‑150 (MQ=255)
                           tCGGATGCGCCTTTCTTTTATCAGACTTGTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGCCGCATGCTGCAAAATCTCAGTTTCTATACGCCAGTTTTGATTGCCAAGACCATCATCGGCCTTATcgtcgt                                                                              <  1:300165/150‑1 (MQ=255)
                               aTGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTGGATTGCCAAGACCATCATCGGCCTTATCGTCGTGTTg                                                                          >  2:51504/1‑150 (MQ=255)
                                                           ccACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTGGATTGCCAAGACCATCATCGGCCTTATCGTCGTGTTGATGGCATTGACCGGTTTCAAAGCCAttt                                              >  2:7995/1‑150 (MQ=255)
                                                                                 cattGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTGGATTGCCAAGACCATCATCTGCCTTATCGTCGTGTTGATGGCATTGACCGGTTTCGAAGCCATTTTCCAGGCCGCCACTCATAGCCa                        <  2:120633/147‑1 (MQ=255)
                                                                                                ttCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTGGATTGCCAAGACCATCATCGGCCTTATCGTCGTGTTGATGGCATTGACCGGTTTCGAAGCCATTTTCCAGGCCGCCACTCATAGCCAGGACGCCGCCAAGGt         <  2:200941/150‑1 (MQ=255)
                                                                                                       ggACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTGGATTGCCAAGACCATCATCGGCCTTATCGTCGTGTTGATGGCATTGACCGGTTTCGAAGCCATTTTCCAGGCCGCCACTCATAGCCAGGACGCCGCCAAGGTGAACGGc  <  1:51504/150‑1 (MQ=255)
                                                                                                                                      cacaAGGTTTGATTGCCAAGACCATCATCGGCCTTATCGTCGTGTTGATGGCATTGACCGGTTTCGAAGCCATTTTCCAGGCCGCCACTCATAGCCAGGa                     >  1:335019/1‑100 (MQ=255)
                                                                                                                                      cacaAGGTTTGATTGCCAAGACCATCATCGGCCTTATCGTCGTGTTGATGGCATTGACCGGTTTCGAAGCCATTTTCCAGGCCGCCACTCATAGCCAGGa                     >  1:335173/1‑100 (MQ=255)
                                                                                                                                      cacaAGGTTTGATTGCCAAGACCATCATCGGCCTTATCGTCGTGTTGATGGCATTGACCGGTTTCGAAGCCATTTTCCAGGCCGCCACTCATAGCCAGGa                     <  2:335019/100‑1 (MQ=255)
                                                                                                                                      cacaAGGTTTGATTGCCAAGACCATCATCGGCCTTATCGTCGTGTTGATGGCATTGACCGGTTTCGAAGCCATTTTCCAGGCCGCCACTCATAGCCAGGa                     <  2:335173/100‑1 (MQ=255)
                                                                                                                                               |                                                                                                             
ACACGTTTGCTAAGAAAAGGCGCATCCTCGGATGCGCCTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTGGATTGCCAAGACCATCATCGGCCTTATCGTCGTGTTGATGGCATTGACCGGTTTCGAAGCCATTTTCCAGGCCGCCACTCATAGCCAGGACGCCGCCAAGGTGAACGGC  >  NC_002947/2633473‑2633725

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