Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,136,737 G→T 25.9% P87T (CCC→ACC)  PP_3638 ← acyl‑CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,136,7370GT25.9% 29.7 / 4.0 15P87T (CCC→ACC) PP_3638acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (8/3);  new base T (2/2);  total (10/5)
Fisher's exact test for biased strand distribution p-value = 5.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGGCCACGCACCTGCAGCGCTTCGATACCGGACAGGAACTTCATGGCGTTGGACAGCCCGGTGGCACCGGCCAATTCACCGGTCAACAGCCGCGGCAGCAGGCGCTCACGCAGGGGCTGGTTGGGGCTGTGCAGCAGGTATTCGATGAAGGCACGCTGGCCCCAGCAGACGAAAGCCGCGGCCAATGAGCGGCTGGCAATGGCGGCGATGGCCTCGACCGCGTCGGTCACCTGCCCGCCCGTACCGCCCAGCGCTGGG  >  NC_002947/4136611‑4136870
                                                                                                                              |                                                                                                                                     
ggcggcCACGCACCTGCAGCGCTTCGATACCGGACAGGAACTTCATGGCGTTGGACAGCCCGGTGGCACCGGCCAATTCACCGGTCAACAGCCGCGGCAGCAGGCGCTCACGCAGGGGCTGGTTGGTGCTGTgcag                                                                                                                              >  1:12424/1‑136 (MQ=255)
ggcggcCACGCACCTGCAGCGCTTCGATACCGGACAGGAACTTCATGGCGTTGGACAGCCCGGTGGCACCGGCCAATTCACCGGTCAACAGCCGCGGCAGCAGGCGCTCACGCAGGGGCTGGTTGGTGCTGTgcag                                                                                                                              >  1:11096/1‑136 (MQ=255)
                  gcgcTTCGATACCGGACAGGAACTTCATGGCGTTGGACAGCCCGGTGGCACCGGCCAATTCACCGGTCAACAGCCGCGGCAGCAGGCGCTCACGCAGGGGCTGGTTGGGGCTGTGCAGCAGGTATTCGATGAAGGc                                                                                                            >  2:188503/1‑136 (MQ=255)
                        cGATACCGGACAGGAACTTCATGGCGTTGGACAGCCCGGTGGCACCGGCCAATTCACCGGTCAACAGCCGCGGCAGCAGGCGCTCACGCAGGGGCTGGTTGGGGCTGTGCAGCAGGTATTCGATGAAGGCACGCTg                                                                                                      >  2:205262/1‑136 (MQ=255)
                          aTACCGGACAGGAACTTCATGGCGTTGGACAGCCCGGTGGCACCGGCCAATTCACCGGTCAACAGCCGCGGCAGCAGGCGCTCACGCAGGGGCTGGTTGGGGCTGTGCAGCAGGTATTCGATGAAGGCACGCTGGc                                                                                                    >  2:239061/1‑136 (MQ=255)
                                   aGGAACTTCATGGCGTTGGACAGCCCGGTGGCACCGGCCAATTCACCGGTCAACAGCCGCGGCAGCAGGCGCTCACGCAGGGGCTGGTTGGGGCTGTGCAGCAGGTATTCGATGAAGGCACGCTGGCCCCAGCAGa                                                                                           >  1:256223/1‑136 (MQ=255)
                                         ttCATGGCGTTGGACAGCCCGGTGGCACCGGCCAATTCACCGGTCAACAGCCGCGGCAGCAGGCGCTCACGCAGGGGCTGGTTGGGGCTGTGCAGCAGGTATTCGATGAAGGCACGCTGGCCCCAGCAGACGAAAg                                                                                     <  1:31414/136‑1 (MQ=255)
                                                             ggTGGCACCGGCCAATTCACCGGTCAACAGCCGCGGCAGCAGGCGCTCACGCAGGGGCTGGTTGGGGCTGTGCAGCAGGTATTCGATGAAGGCACGCTGGCCCCAGCAGACGAAAGCCGCGGCCAATGAGCGGCTg                                                                 >  2:159862/1‑136 (MQ=255)
                                                             ggTGGCACCGGCCAATTCACCGGTCAACAGCCGCGGCAGCAGGCGCTCACGCAGGGGCTGGTTGGGGCTGTGCAGCAGGTATTCGATGAAGGCACGCTGGCCCCAGCAGACGAAAGCCGCGGCCAATGAGCGGCTg                                                                 >  2:142345/1‑136 (MQ=255)
                                                                       gCCAATTCACCGGTCAACAGCCGCGGCAGCAGGCGCTCACGCAGGGGCTGGTTGGGGCTGTGCAGCAGGTATTCGATGAAGGCACGCTGGCCCCAGCAGACGAAAGCCGCGGCCAATGAGCGGCTGGCAATggcgg                                                       >  1:164422/1‑136 (MQ=255)
                                                                                                 cagcagGCGCTCACGCAGGGGCTGGTTGGGGCTGTGCAGCAGGTATTCGATGAAGGCACGCTGGCCCCAGCAGACGAAAGCCGCGGCCAATGAGCGGCTGGCAATGGCGGCGATGGCCTCGACCGCGTCGGTCAcc                             >  2:250631/1‑136 (MQ=255)
                                                                                                               gCAGGGGCTGGTTGGGGCTGTGCAGCAGGTATTCGATGAAGGCACGCTGGCCCCAGCAGACGAAAGCCGCGGCCAATGAGCGGCTGGCAATGGCGGCGATGGCCTCGACCGCGTCGGTCACCTGCCCGCCCGTAcc               <  1:159862/136‑1 (MQ=255)
                                                                                                               gCAGGGGCTGGTTGGGGCTGTGCAGCAGGTATTCGATGAAGGCACGCTGGCCCCAGCAGACGAAAGCCGCGGCCAATGAGCGGCTGGCAATGGCGGCGATGGCCTCGACCGCGTCGGTCACCTGCCCGCCCGTAcc               <  1:142345/136‑1 (MQ=255)
                                                                                                                          ttGGTGCTGTGCAGCAGGTATTCGATGAAGGCACGCTGGCCCCAGCAGACGAAAGCCGCGGCCAATGAGCGGCTGGCAATGGCGGCGATGGCCTCGACCGCGTCGGTCACCTGCCCGCCCGTACCGCCCAGCGCTg    <  2:12424/136‑1 (MQ=255)
                                                                                                                          ttGGTGCTGTGCAGCAGGTATTCGATGAAGGCACGCTGGCCCCAGCAGACGAAAGCCGCGGCCAATGAGCGGCTGGCAATGGCGGCGATGGCCTCGACCGCGTCGGTCACCTGCCCGCCCGTACCGCCCAGCGCTg    <  2:11096/136‑1 (MQ=255)
                                                                                                                            ggggCTGTGCAGCAGGTATTCGATGAAGGCACGCTGGCCCCAGCAGACGAAAGCCGCGGCCAATGAGCGGCTGGCAATGGCGGCGATGGCCTCGACCGCGTCGGTCACCTGCCCGCCCGTACCGCCCAGCGCTggg  >  2:162068/1‑136 (MQ=255)
                                                                                                                            ggggCTGTGCAGCAGGTATTCGATGAAGGCACGCTGGCCCCAGCAGACGAAAGCCGCGGCCAATGAGCGGCTGGCAATGGCGGCGATGGCCTCGACCGCGTCGGTCACCTGCCCGCCCGTACCGCCCAGCGCTggg  >  1:144870/1‑136 (MQ=255)
                                                                                                                            ggggCTGTGCAGCAGGTATTCGATGAAGGCACGCTGGCCCCAGCAGACGAAAGCCGCGGCCAATGAGCGGCTGGCAATGGCGGCGATGGCCTCGACCGCGTCGCTCACCTGCCCGCCCGTACCGCCCAGCGCTggg  >  2:213251/1‑136 (MQ=255)
                                                                                                                              |                                                                                                                                     
GGCGGCCACGCACCTGCAGCGCTTCGATACCGGACAGGAACTTCATGGCGTTGGACAGCCCGGTGGCACCGGCCAATTCACCGGTCAACAGCCGCGGCAGCAGGCGCTCACGCAGGGGCTGGTTGGGGCTGTGCAGCAGGTATTCGATGAAGGCACGCTGGCCCCAGCAGACGAAAGCCGCGGCCAATGAGCGGCTGGCAATGGCGGCGATGGCCTCGACCGCGTCGGTCACCTGCCCGCCCGTACCGCCCAGCGCTGGG  >  NC_002947/4136611‑4136870

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: