Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,310,151 T→G 16.1% L408R (CTG→CGG) ‡ lhgO → L‑2‑hydroxyglutarate oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,310,1510TG16.1% 48.9 / 5.6 25L408R (CTG→CGG) ‡lhgOL‑2‑hydroxyglutarate oxidase
Reads supporting (aligned to +/- strand):  ref base T (12/9);  new base G (2/2);  total (14/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 = 7.84e-01

TCGTCAACACCGCACGCAGTGTCAACGTGTGCAACGCACCGTCGCCAGCTGCTACCTCGGCCATCCCGATCGGCGCCTACATCGTCGAGAAGGTGTGCGAGCAGGTCGGCCTGCAGGGCGGCAGCTTCCCCAAGGCCGACCTGTCGGCCAACCAGAGGGCTGCCAGCTGACAGCCGGGTAACACCTTCTGCGGTTTCGTTTTACTCAACCCGTACTGCCCTGCTCACCATTCAGGGCAGGGCGGGCTTTTTCGCGCCTGCGTCGGTGACGTGCAGCCGGCGGC  >  NC_002947/3310010‑3310292
                                                                                                                                             |                                                                                                                                             
tCGTCAACACCGCACGCAGTGTCAACGTGTGCAACGCACCGTCGCCAGCTGCTACCTCGGCCATCCCGATCGGCGCCTACATCGTCGAGAAGGTGTGCGAGCAGGTCGGCCTGCAGGGCGGCAGCTTCCCCAAGGCCGACAGGTCGGcct                                                                                                                                       <  1:231573/150‑2 (MQ=255)
  gTCAACACCGCACGCAGTGTCAACGTGTGCAACGCACCGTCGCCAGCTGCTACCTCGGCCATCCCGATCGGCGCCTACATCGTCGAGAAGGTGTGCGAGCAGGTCGGCCTGCAGGGCGGCAGCTTCCCCAAGGCCGACAGGTCGGcct                                                                                                                                       >  2:76724/1‑147 (MQ=255)
  gTCAACACCGCACGCAGTGTCAACGTGTGCAACGCACCGTCGCCAGCTGCTACCTCGGCCATCCCGATCGGCGCCTACATCGTCGAGAAGGTGTGCGAGCAGGTCGGCCTGCAGGGCGGCAGCTTCCCCAAGGCCGACAGGTCGGcct                                                                                                                                       <  1:76724/148‑2 (MQ=255)
   tCAACACCGCACGCAGTGTCAACGTGTGCAACGCACCGTCGCCAGCTGCTACCTCGGCCATCCCGATCGGCGCCTACATCGTCGAGAAGGTGTGCGAGCAGGTCGGCCTGCAGGGCGGCAGCTTCCCCAAGGCCGACCTGTCGGCCAAcc                                                                                                                                    >  1:417068/1‑150 (MQ=255)
   tCAACACCGCACGCAGTGTCAACGTGTGCAACGCACCGTCGCCAGCTGCTACCTCGGCCATCCCGATCGGCGCCTACATCGTCGAGAAGGTGTGCGAGCAGGTCGGCCTGCAGGGCGGCAGCTTCCCCAAGGCCGACCTGTCGGCCAAc                                                                                                                                     <  2:1000/149‑1 (MQ=255)
   tCAACACCGCACGCAGTGTCAACGTGTGCAACGCACCGTCGCCAGCTGCTACCTCGGCCATCCCGATCGGCGCCTACATCGTCGAGAAGGTGTGCGAGCAGGTCGGCCTGCAGGGCGGCAGCTTCCCCAAGGCCGACCTGTCGGCCAAc                                                                                                                                     >  1:1000/1‑149 (MQ=255)
                                 aCGCACCGTCGCCAGCTGCTACCTCGGCCATCCCGATCGGCGCCTACATCGTCGAGAAGGTGTGCGAGCAGGTCGGCCTGCAGGGCGGCAGCTTCCCCAAGGCCGACCTGTCGGCCAACCAGAGGGCTGCCAGCTGACAGCCGGGTAaca                                                                                                      >  2:730579/1‑150 (MQ=255)
                                   gCACCGTCGCCAGCTGCTACCTCGGCCATCCCGATCGGCGCCTACATCGTCGAGAAGGTGTGCGAGCAGGTCGGCCTGCAGGGCGGCAGCTTCCCCAAGGCCGACCTGTCGGCCAACCAGAGGGCTGCCAg                                                                                                                       <  2:731019/131‑1 (MQ=255)
                                   gCACCGTCGCCAGCTGCTACCTCGGCCATCCCGATCGGCGCCTACATCGTCGAGAAGGTGTGCGAGCAGGTCGGCCTGCAGGGCGGCAGCTTCCCCAAGGACGACCTGTCGGCCAACCAGAGGGCTGCCAg                                                                                                                       >  1:731019/1‑131 (MQ=255)
                                         tCGCCAGCTGCTACCTCGGCCATCCCGATCGGCGCCTACATCGTCGAGAAGGTGTGCGAGCAGGTCGGCCTGCAGGGCGGCAGCTTCCCCAAGGCCGACCTGTCGGCCAACCAGAGGGCTGCCAGCTGACAGCCGGGTAACACCTTCTGc                                                                                              >  2:105664/1‑150 (MQ=255)
                                            ccAGCTGCTACCTCGGCCATCCCGATCGGCGCCTACATCGTCGAGAAGGTGTGCGAGCAGGTCGGCCTGCAGGGCGGCAGCTTCCCCAAGGCCGACCTGTCGGCCAACCAGAGGGCTGCCAGCTGACAGCCGGGTAACACCTTCTGCGGt                                                                                           >  1:434772/1‑150 (MQ=255)
                                                     aCCTCGGCCATCCCGATCGGCGCCTACATCGTCGAGAAGGTGTGCGAGCAGGTCGGCCTGCAGGGCGGCAGCTTCCCCAAGGCCGACCTGTCGGCCAACCAGAGGGCTGCCAGCTGACAGCCGGGTAACACCTTCTGCGGTTTCGTTTTa                                                                                  >  1:566277/1‑150 (MQ=255)
                                                                 ccGATCGGCGCCTACATCGTCGAGAAGGTGTGCGAGCAGGTCGGCCTGCAGGGCGGCAGCTTCCCCAAGGCCGACCTGTCGGCCAACCAGAGTGCTGTCATATGACAGCCGGGTAACACCTTCTGCGGTTTCGTTTTACTCAACCCGTAc                                                                      <  2:195446/150‑1 (MQ=255)
                                                                  cGATCGGCGCCTACATCGTCGAGAAGGTGTGCGAGCAGGTCGGCCTGCAGGCCGGCAGCTTCCCCAAGGCCGACCTGTCGGCCAACCAGAGGGCTGCCAGCTGACAGCCGGGTAACACCTTCTGCGGTTTCGTTTTACTCAACCCGTACt                                                                     >  2:359564/1‑150 (MQ=255)
                                                                        gcgcCTACATCGTCGAGAAGGTGTGCGAGCAGGTCGGCCTGCAGGGCGGCAGCTTCCCCAAGGCCGACCTGTCGGCCAACCAGAGGGCTGCCAGCTGACAGCCGGGTAac                                                                                                       <  1:618289/110‑1 (MQ=255)
                                                                        gcgcCTACATCGTCGAGAAGGTGTGCGAGCAGGTCGGCCTGCAGGGCGGCAGCTTCCCCAAGGCCGACCTGTCGGCCAACCAGAGGGCTGCCAGCTGACAGCCGGGTAac                                                                                                       >  2:618289/1‑110 (MQ=255)
                                                                                              gtgCGAGCAGGTCGGCCTGCAGGGCGGCAGCTTCCCCAAGGCCGACCTGTCGGCCAACCAGAGGGCTGCCAGCTGACAGCCGGGTAACACCTTCTGCGGTTTCGTTTTACTCAACCCGTACTGCCCTGCTCACCATTCAGGGCAgggcgg                                         >  2:54451/1‑150 (MQ=255)
                                                                                                  gAGCAGGTCGGCCTGCAGGGCGGCAGCTTCCCCAAGGCCGAACTATAGGCAAATCAGAGGGCTGCCAGCTGACAGACGGGTAACACCTTCTGCGGTTTCGTTTTACTCAACCCGTACTGCCCTGCTCACCATTCAGGGCAGGGCGGGCtt                                     <  1:776613/150‑1 (MQ=255)
                                                                                                         tCGGCCTGCAGGGCGGCAGCTTCCCCAAGGCCGACCTGTCGGCCAACCAGAGGGCTGCCAGCTGACAGCCGGGTAACACCTTCTGCGGTTTCGTTTTACTCAACCCGTACTGCCCTGCTCACCATTCAGGGCAGGGCGGGCTTTTTcgcg                              >  1:614199/1‑150 (MQ=255)
                                                                                                          cGGCCTGCAGGGCGGCAGCTTCCCCAAGGCCGACCTGTCGGCCAACCAGAGGGCTGCCAGCTGACAGCCGGGTAACACCTTCTGCGGTTTCGTTTTACTCAACCCGTACTGCCCTGCTCACCATTCAGGGCAGGGCGGGCTTTTTCgcgc                             <  1:140762/150‑1 (MQ=255)
                                                                                                                  aGGGCGGCAGCTTCCCCAAGGCCGACCTGTCGGCCAACCAGAGGGCTGCCAGCTGACAGCCGGGTAACACCtt                                                                                                  <  2:171780/73‑1 (MQ=255)
                                                                                                                  aGGGCGGCAGCTTCCCCAAGGCCGACCTGTCGGCCAACCAGAGGGCTGCCAGCTGACAGCCGGGTAACACCtt                                                                                                  >  1:171780/1‑73 (MQ=255)
                                                                                                                            cTTCCCCAAGGCCGACCTGTCGGCCAACCAGAGGGCTGCCAGCTGACAGCCGGGTAACACCTTCTGCGGTTTCGTTTTACTCAACCCGTACTGCCCTGCTCACCATTCAGGGCAGGGCGGGCTTTTTCGCGCCTGCGTCGGTGACGTGCa           <  2:742708/150‑1 (MQ=255)
                                                                                                                                  gttggCCGACAGGTCGGCCAACCAGAGGGCTGCCAGCTGACAGCCGGGTAACACCTTCTGCGGTTTCGTTTTACTCAACCCGTACTGCCCTGCTCACCATTCAGGGCAGGGCGGGCTTTTTCGCGCCTGCGTCGGTGACGTGCAGCcggc     >  2:561890/4‑150 (MQ=255)
                                                                                                                                     ggCCGACCTGTCGGCCAACCAGAGGGCTGCCAGCTGACAGCCGGGTAACACCTTCTGCGGTTTCGTTTTACTCAACCCGTACTGCCCTGCTCACCATTCAGGGCAGGGCGGGCTTTTTCGCGCCTGCGTCGGTGACGTGCAGCCggcggc  <  2:566277/150‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                             
TCGTCAACACCGCACGCAGTGTCAACGTGTGCAACGCACCGTCGCCAGCTGCTACCTCGGCCATCCCGATCGGCGCCTACATCGTCGAGAAGGTGTGCGAGCAGGTCGGCCTGCAGGGCGGCAGCTTCCCCAAGGCCGACCTGTCGGCCAACCAGAGGGCTGCCAGCTGACAGCCGGGTAACACCTTCTGCGGTTTCGTTTTACTCAACCCGTACTGCCCTGCTCACCATTCAGGGCAGGGCGGGCTTTTTCGCGCCTGCGTCGGTGACGTGCAGCCGGCGGC  >  NC_002947/3310010‑3310292

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