Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 4,135,607 A→C 48.9% Q32P (CAG→CCG)  tgnC → (Z)‑2‑((N‑methylformamido)methylene)‑5‑hydroxybu tyrolactone dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1224,135,6070AC48.9% ‑0.8 / 15.4 26Q32P (CAG→CCG) tgnC(Z)‑2‑((N‑methylformamido)methylene)‑5‑hydroxybu tyrolactone dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (5/7);  new base C (9/4);  total (14/12)
Fisher's exact test for biased strand distribution p-value = 2.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

GCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAGCTGCTGCGTCGTCTG  >  ACB122/4135479‑4135740
                                                                                                                                |                                                                                                                                     
gCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCACGCCTg                                                                                                                                 >  4:373517/1‑135 (MQ=255)
gCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGAACCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAAACCGGCCACCCCCCAGGCCTg                                                                                                                                 >  7:190459/1‑135 (MQ=255)
    ccTCACCGAATCCGCCAAGGGGATTGTTGCATGACCATCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGATCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCCCCCCGGCCAGGgcg                                                                                                                             >  5:185915/1‑135 (MQ=255)
     cTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCCCCCCCGCCAGGgcgc                                                                                                                            >  7:89872/1‑135 (MQ=255)
     cTCACCGAATACGCCAAGGGGATTGTAGCATGACCCTCGTTCGTTTCCAGATGTGCATCGAAGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCCCCCCCCAGGCCTGGgcg                                                                                                                             >  2:47632/1‑134 (MQ=255)
          cGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACCGCCTCAACCCGGCCACCGCCCCGGCCGGGGGGCCAcgg                                                                                                                       >  6:220273/1‑133 (MQ=255)
              tCCGCCAAGGGTTTTGTTGCAGGACCCTGGTTCGTTTACAGATGTGCATCGACGGCCAATGGCGCGATGCACAGTGCTGCAAGACCTTCGACAGCCTAAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCg                                                                                                                   <  8:177514/135‑1 (MQ=255)
                 gCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCCGCCTGGGCGCCACTGCccgccg                                                                                                                >  2:237323/1‑135 (MQ=255)
                  ccAAGGGTTTTGTTGCTTGACCCTCGTTCTTTTCAAGATGTCCATCGACGGCAAATGGCGCGATGCCCAGAGCGGAAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGc                                                                                                               <  2:94703/135‑1 (MQ=255)
                   cAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCAGGGCGCCACTGCCCgccgcc                                                                                                              >  6:24521/1‑135 (MQ=255)
                   cAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCCCCCCGGCCAGCGCGCAACGGCCCCAccc                                                                                                               >  4:374049/1‑132 (MQ=255)
                          tttgttgCATCACCCGCGTTCGTTTCCAGATGTGCAGCGACGCCCAAGGGCGCGCTGCCCAGAGCGGAAACACCTTCGCCAGCCTCAACCCGCCCACCGCCCAGGCCGGGGCGCAACTGCCCGACGCCGATGAAg                                                                                                       <  8:241233/134‑1 (MQ=255)
                             gttgCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCGGGGCGCCCCAGCCCGACGCCGCTGAAGc                                                                                                      >  5:342207/1‑133 (MQ=255)
                                                           gCATCCCCGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGc                                                                      <  4:165782/135‑1 (MQ=255)
                                                             aTCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCACGCCCGCGCGCCGCGGGCCGACGCCGAGCCAGCCGAGGTGGGCCTGGCCGGTCACGCCGCCCCGcca                                                                    >  5:324754/1‑133 (MQ=255)
                                                               cGACGGCCAATGGCGCGATGCCCAGAGCGGCGAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCTGGGCGCCACTGCCCGCCGCCGCTGCAGCCGCTGTGGAGCTGGCCGTTCAGGCCGCCCCGCGTg                                                                   >  8:130886/1‑134 (MQ=255)
                                                               cGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGGAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGc                                                                  >  2:285739/1‑135 (MQ=255)
                                                                     ccAAGGGCGCGATCCCCGTAGCGCCCAGACCTGCGCCCCCCCCACCCCGCCCACCGCCCTGGCGTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGa                                                            <  6:337504/135‑1 (MQ=255)
                                                                                 tccccAGCGCGCCAACCCCTCCGACCCCCCCACCCCGCCCCCCCCCCCGGCCGGGCCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTg                                                <  7:30252/133‑1 (MQ=255)
                                                                                 cgCCCAGCGCGCCAACCCCTTCGACCGCCCCCCCCCGCCCCCCCCCCCGCCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTg                                                <  5:135354/134‑1 (MQ=255)
                                                                                     ccgagCGGAAACACCTTCGCCCGCCACAACCCCCCCCCCGCCCAGCCCTGGCCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGgcgc                                            <  3:373517/133‑1 (MQ=255)
                                                                                            gCAAGACCTTCTACAGCCTCAACCCGGCCACCGCCCAGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCa                                     <  8:190459/135‑1 (MQ=255)
                                                                                                gACCTTCGACAGCCTCAACCCGGCCACCCCCCCCGCCGGGGGGCCCCTCGCCGCCGCCCAGGAACGCGCTCTTGAGCAGGCCCGTCAGGCCGCCCCGCGTGCCGTCGCCAGCAAGGCGGGGCGCCGCCTccaccc                                 >  2:175752/1‑130 (MQ=255)
                                                                                                                     gccaaCCCCCCCGCCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTGAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAgctgct            <  5:34651/131‑1 (MQ=255)
                                                                                                                        cACCGCCCCGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAGCtgctgct         <  1:101055/135‑2 (MQ=255)
                                                                                                                               cAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAGCTGCTGCGTCGTCTg  <  4:304553/135‑1 (MQ=255)
                                                                                                                                |                                                                                                                                     
GCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAGCTGCTGCGTCGTCTG  >  ACB122/4135479‑4135740

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