Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,533,716 T→G 16.6% I117M (ATT→ATG atoA → 3‑oxoacid CoA‑transferase subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,533,7160TG16.6% 47.7 / 2.5 24I117M (ATT→ATGatoA3‑oxoacid CoA‑transferase subunit A
Reads supporting (aligned to +/- strand):  ref base T (11/9);  new base G (2/2);  total (13/11)
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.83e-01

GAAGACCGGCAGATTCGCAAGATGATCGCCTCCTACGTGGGCGAAAACGCCGAGTTCGAACGCCAGTTGCTCAGTGGCGAGCTGGAAGTGGAACTGACCCCACAGGGCACCCTCGCCGAAAAAATGCGCGCAGGTGGTGCCGGCATTCCGGCGTTCTACACCGCCACCGGCTATGGCACCCCGGTTGCCGATGGCAAGGAAGTCCGCGAGTTCAAGGGCCGCAAGTACATCCTGGAAGAATCCATCACCGGTGACTTCGCCATCGTCAAGGGCTGGAAGGCC  >  NC_002947/3533570‑3533851
                                                                                                                                                  |                                                                                                                                       
gAAGCCCGGCAGATTCGCAAGATGATCGCCTCCTACGTGGGCGAAAACTCCGAGTTCGAACGCCAGTTGCTCAGTGGCGAGCTGGAAGTGGAACTGACCCCACAGGGCACCCTCGCCGAAAAAATGCGCGCAGGTGGTGCCGGCATTCCg                                                                                                                                      <  1:55414/150‑1 (MQ=255)
                    gatgatCGCCTCCTACGTGGGCGAAAACGCCGAGTTCGAACGCCAGTTGCTCAGTGGCGAGCTGGAAGTGGAACTGACCCCACAGGGCACCCTCGCCGAAAAAATGCGCGCAGGTGGTGCCGGCATTCCGGCGTTCTACACCGCCACCgg                                                                                                                  >  1:167298/1‑150 (MQ=255)
                              tcctACGTGGGCGAAAACGCCGAGTTCGAACGCCAGTTGCTCAGTGGCGAGCTGGAAGTGGAACTGACCCCACAGGGCACCCTCGCCGAAAAAATGCGCGCAGGTGGTGCCGGCATTCCGGCGTTCTACACCGCCACCGGCTATGGCAcc                                                                                                        >  2:302097/1‑150 (MQ=255)
                                                gCCGAGTTCGAACGCCAGTTGCTCAGTGGCGAGCTGGAAGTGGAACTGACCCCACAGGGCACCCTCGCCGAAAAAATGCGCGCAGGTGGTGCCGGCATTCCGGCGTTCTACACCGCCACCGGCTATGGCACCCCGGTTGCCGATGGCAAg                                                                                      >  1:174387/1‑150 (MQ=255)
                                                         gAACGCCAGTTGCTCAGTGGCGAGCTGGAAGTGGAACTGACCCCACAGGGCACCCTCGCCGAAAAAATGCGCGCAGGTGGTGCCGGCATTCCGGCGTTCTACACCGCCACCGGCTATGGCACCCCGGTTGCCGATGGCAAGGAAGTCcgc                                                                             >  2:356659/1‑150 (MQ=255)
                                                               cAGTTGCTCAGTGGCGAGCTGGAAGTGGAACTGACCCCACAGGGCACCCTCGCCGAAAAAATGCGCGCAGGTGGTGCCGGCATTCCGGCGTTCTACACCGCCACCGGCTATGGCACCCCGGTTGCCGATGGCAAGGAAGTCCGCGAGTTc                                                                       >  2:125595/1‑150 (MQ=255)
                                                                         gTGGCGAGCTGGAAGTGGAACTGACCCCACAGGGCACCCTCGCCGAAAAAATGCGCGCAGGTGGTGCCGGCATTCCGGCGTTCTACACCGCCACCGGCTATGGCACCCCGGTTGCCGATGGCAAGGAAGTCCGCGAGTTCAAGGGCCGCa                                                             >  1:237638/1‑150 (MQ=255)
                                                                              gAGCTGGAAGTGGAACTGACCCCACAGGGCACCCTCGCCGAAAAAATGCGCGCAGGTGGTGCCGGCATTCCGGCGTTCTACACCGCCACCGGCTATGGCACCCCGGTTGCCGATGGCAAGGAAGTCcgc                                                                             >  1:24715/1‑129 (MQ=255)
                                                                              gAGCTGGAAGTGGAACTGACCCCACAGGGCACCCTCGCCGAAAAAATGCGCGCAGGTGGTGCCGGCATTCCGGCGTTCTACACCGCCACCGGCTATGGCACCCCGGTTGCCGATGGCAAGGAAGTCcgc                                                                             <  2:24715/129‑1 (MQ=255)
                                                                                   ggAAGTGGAACTGACCCCACAGGGCACCCTCGCCGAAAAAATGCGCGCAGGTGGTGCCGGCATTCCGGCGTTCTACACCGCCACCGGCTATGGCACCCCGGTTGCCGATGGCAAGGAAGTCCGCGAGTTCAAGGGCCGCAAGTACATCCt                                                   <  1:257143/150‑1 (MQ=255)
                                                                                         ggAACTGACCCCACAGGGCACCCTCGCCGAAAAAATGCGCGCAGGTGGTGCCGGCATTCCGGCGTTCTACACCGCCACCGGCTATGGCACCCCGGTTGCCGATGGCAAGGAAGTCCGCGAGTTCAAGGGCCGCAAGTACATCCTGgaaga                                             >  1:224617/1‑150 (MQ=255)
                                                                                                  cccACAGGGCACCCTCGCCGAAAAAATGCGCGCAGGTGGTGCCGGCATTCCGGCGTTCTACACCGCCACCGGCTATGGCACCCCGGTTGCCGATGGCAAGGAAGTCCGCGAGTTCAAGGGCCGCAAGTACATCCTg                                                  <  1:59151/136‑1 (MQ=255)
                                                                                                  cccACAGGGCACCCTCGCCGAAAAAATGCGCGCAGGTGGTGCCGGCATTCCGGCGTTCTACACCGCCACCGGCTATGGCACCCCGGTTGCCGATGGCAAGGAAGTCCGCGAGTTCAAGGGCCGCAAGTACATCCTg                                                  >  2:59151/1‑136 (MQ=255)
                                                                                                   ccACAGGTCACCCTCGCCGAAAAAATGCGCGCAGGTGGTGCCGGCATTCCGGCGTTCTACACCGCCACCGGCTATGGCACCCCGGTTGCCGATGGCAAGGAAGTCCGCGAGTTCAAGGGCCGCAAGTACATCCTGGAAGAATCCATCAcc                                   >  2:93143/1‑150 (MQ=255)
                                                                                                    cacaGGGCACCCTCGCCGAAAAAATGCGCGCAGGTGGTGCCGGCATTCCGGCGTTCTACACCGCCACCGGCTATGGCACCCCGGTTGCCGATGGCAAGGAAGTCCGCGAGTTCAAGGGCCGCAAGTACATCCTGGAAGAATCCATCACCg                                  <  2:237638/150‑1 (MQ=255)
                                                                                                     acaGGGCACCCTCGCCGAAAAAATGCGCGCAGGTGGTGCCGGCATTCCGGCGTTCTACACCGCCACCGGCTATGGCACCCCGGTTGCCGATGGCAAGGAAGTCCGCGAGTTCAAGGGCCGCAAGTACATCCTGGAAGAATCCATCACCgg                                 <  1:302097/150‑1 (MQ=255)
                                                                                                                tCGCCGAAAAAATGCGCGCAGGTGGTGCCGGCATTCCGGCGTTCTACACCGCCACCGGCTATGGCACCCCGGTTGCCGATGGCAAGGAAGTCCGCGAGTTCAAGGGCCGCAAGTACATCCTGGAAGAATCCATCACCGGTGACTTCGCCa                      >  1:163499/1‑150 (MQ=255)
                                                                                                                 cGCCGAAAAAATGCGCGCAGGTGGTGCCGGCATTCCGGCGTTCTACACCGCCACCGGCTATGGCACCCCGGTTGCCGATGGCAAGGAAGTCCGCGAGTTCAAGGGCCGCAAGTACATCCTGGAAGAATCCATCACCGGTGACTTCGCCAt                     <  1:93143/150‑1 (MQ=255)
                                                                                                                  gCCGAAAAAATGCGCGCAGGTGGTGCCGGCATTCCGGCGTTCTACACCGCCACCGGCTATGGCACCCCGGTTGCCGATGGCAAGGAAGTCCGCGAGTTCAAGGGCCGCAAGTACATCCTGGAAGAATCCATCACCGGTGACTTCGCCATc                    <  1:356659/150‑1 (MQ=255)
                                                                                                                                    ggtggtGCCGGCATTCCGGCGTTCTACACCGCCACCGGCTATGGCACCCCGGTTGCCGATGGCAAGGAAGTCCGCGAGTTCAAGGGCCGCAAGTACATCCTGGAAGAATCCATCACCGGTGACTTCGCCATCGTCAAGGGCTGGAAGGcc  <  2:163499/150‑1 (MQ=255)
                                                                                                                                      gaacgCCGGAATGCCGGCGTCCTACACCGCCACCGGCTATGGCACCCCGGTTGCCGATGGCAAGGAAGTCCGCGAGTTCAAGGGCCGCAAGTACATCCTGgaa                                               >  2:206479/5‑103 (MQ=255)
                                                                                                                                      gaacgCCGGAATGCCGGCGTCCTACACCGCCACCGGCTATGGCACCCCGGTTGCCGATGGCAAGGAAGTCCGCGAGTTCAAGGGCCGCAAGTACATCCTGgaa                                               <  1:206479/99‑1 (MQ=255)
                                                                                                                                       aacgCCGGAATGCCGGCGTTCTACACCGCCACCGGCTATGGCACCCCGGTTGCCGATGGCAAGGAAGTCCGCGAGTTCAAGGGCCGCAAGTACATCCTGGAAGAATCCATCACCGGt                                >  1:346248/4‑117 (MQ=255)
                                                                                                                                       aacgCCGGAATGCCGGCGTTCTACACCGCCACCGGCTATGGCACCCCGGTTGCCGATGGCAAGGAAGTCCGCGAGTTCAAGGGCCGCAAGTACATCCTGGAAGAATCCATCACCGGt                                <  2:346248/114‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                       
GAAGACCGGCAGATTCGCAAGATGATCGCCTCCTACGTGGGCGAAAACGCCGAGTTCGAACGCCAGTTGCTCAGTGGCGAGCTGGAAGTGGAACTGACCCCACAGGGCACCCTCGCCGAAAAAATGCGCGCAGGTGGTGCCGGCATTCCGGCGTTCTACACCGCCACCGGCTATGGCACCCCGGTTGCCGATGGCAAGGAAGTCCGCGAGTTCAAGGGCCGCAAGTACATCCTGGAAGAATCCATCACCGGTGACTTCGCCATCGTCAAGGGCTGGAAGGCC  >  NC_002947/3533570‑3533851

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: