Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 3,784,646 T→G 35.2% D571A (GAC→GCC)  dmdC_4 ← 3‑methylmercaptopropionyl‑CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1433,784,6460TG35.2% 4.0 / 2.1 22D571A (GAC→GCC) dmdC_43‑methylmercaptopropionyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (7/7);  new base G (4/4);  total (11/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 = 8.90e-01

TCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGCCGCGTCGATGGCCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCTGCTTTTCCCGGCGCTGCGCCGCCCCCACTGGCAACGCC  >  ACB143/3784524‑3784750
                                                                                                                          |                                                                                                         
tCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGCCGCGGCGATGgccgcccc                                                                                             >  7:83114/1‑134 (MQ=255)
tCGGCGGACCGCGGCATGCGCACCCGCCCCTCGGCAGGCCAGCCGACGCCGGTGGTCATGCCGCCCTCGCAGGGGGGTTGCAGACCGATCCGGCAATGAATGGCACCGGCGTGGCGGCCGCGTCGATGGCCGCCAc                                                                                             <  4:235942/136‑1 (MQ=255)
       aCCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGCCGCGGGGGTGGGCGCCCCCAGcccg                                                                                      >  7:160381/1‑134 (MQ=255)
          gCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGCCGCGGCGATGGCCGCCCCCGGCCGGGCg                                                                                   >  8:304923/1‑136 (MQ=255)
                   gcTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATGCGGCAATGAATGGCAACGGCGTGGCGGCCCCGTGGATGGCCCCCCCCCGCCCGGCGCgggcgggc                                                                          >  1:121762/1‑136 (MQ=255)
                         gCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGCCGCGTCGATGGCCGCCCCCGGCCGGTCGAGGTCGGGCCGCAgg                                                                    >  6:167765/1‑136 (MQ=255)
                                ggCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGCCGCGGGGGTGGCCGCCACCTGCCGGTCGAGGGCGGGCCGCAGGtttccgg                                                             >  2:264888/1‑130 (MQ=255)
                                 gCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGTTTCGGAAATGAATGGCAACGGCGTGGCGGCCGCGTCTATGGCCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAg                                                            <  7:382270/136‑1 (MQ=255)
                                 gCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCATTGAATGGCAACGGCGGGGCGGCCGCGGCGATGGCCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAg                                                            <  7:304687/136‑1 (MQ=255)
                                  cACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGCCCCGTGGGATGGCCCCCACAGCCGGTCCATGTGGGGCACCAGGTATCCgcc                                                            >  7:372926/1‑133 (MQ=255)
                                        agccgcagccgGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGCCGCGTGGATGGCCGCCCCCAGCCCGGCGAGGGGGGGCCGCCGGGGGCCGAGGAAGAAc                                                     >  4:81914/1‑136 (MQ=255)
                                                     gttattGACGACATTGCACGTGGGTTTCATACCTATTCGGCATTTAATGGAAACGGGGTGGGGGCGGCGTCGATGGCCGACACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCt                                        <  8:322510/131‑1 (MQ=255)
                                                      ttCATTACGCCCTCGCACGTGGGTTTCAGAACTATTCGCCAATGAATGGCACGGGGGTGGCGGCCGCGTCGATGGCCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCt                                        <  1:323877/135‑1 (MQ=255)
                                                                 cTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGGCCCGTCCGAGGCCCCCCCCAGCCCGGCGAGGTCGGGGCCCAGGGTGCCCCGGCAGCAACCAACGCGTTGCAGCTTTTaccc                              >  1:377885/1‑133 (MQ=255)
                                                                   cGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGCCGCGTCGATGGCCGCCCCCAGCCGGTCGAGGTCGGGCCGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCGGCTTTTACCCgcgc                          >  2:127189/1‑136 (MQ=255)
                                                                             tttCATAACGATTCGGAAATTAAGGGCAAGGGGGGGGGGGCGGCGTCGATGGCCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCTGCTTTTCCCGGCGCTGCGCCGccc                <  2:280585/136‑1 (MQ=255)
                                                                                  taaCTATGCGGAAATGATGGGCACGGGGGGGGGGGCCGGGTCGAGGGCCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCTGCTTTTCCCGGCGCTGCGCCGCCCCCACt           <  2:330186/135‑1 (MQ=255)
                                                                                  caaCTATTCGGCAATCAATGGCACCGGCGGGGCGGCCGCTGCGATGGCCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCTGCTTTTCCCGGCGCTGCGCCGCCCCCACt           <  7:22864/135‑1 (MQ=255)
                                                                                    actaTTCGGCAATGAATGGCAACGGCGGGGCGGCCGCGTCGATGGCCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCTGCTTTTCCCGGCGCTGCGCCGCCCCCACTgg         <  1:48253/133‑1 (MQ=255)
                                                                                    aCGATTGGGCAATTATTGGCACCGGCGGGGGGGCGGGGGCGATGGCCCCCACCTGCCGGTCGCGGTCGGGCTGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCTTCTTTTCCCCGCGCTGCGCCGCCCCCACt           <  6:78307/134‑1 (MQ=255)
                                                                                    aCGATTCGGCAATGAATGGCAACGGCGGGGCGGCCGCGGCGTGGGCCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCTGCTTTTCCCGGCGCTGCGCCGCCCCCACTgg         <  7:19142/136‑1 (MQ=255)
                                                                                           ggCAATGAATGGCAACGGCGTGGCGGCCGCGTCGATGGCCGCCACCAGCCGGGCGAGGTGGGGCACCACGTAGCCCAGGAAGAACCCCACGCTTTACTGCGTTTAACCGCGCATCCCCCCCCCCACAGGCAAcccc  >  6:132581/1‑136 (MQ=255)
                                                                                                                          |                                                                                                         
TCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGCCGCGTCGATGGCCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCTGCTTTTCCCGGCGCTGCGCCGCCCCCACTGGCAACGCC  >  ACB143/3784524‑3784750

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: