Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 356,383 G→A 12.7% V50V (GTC→GTT)  lcdH_2 ← L‑carnitine dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT2440356,3830GA12.7% 103.9 / 4.8 39V50V (GTC→GTT) lcdH_2L‑carnitine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (12/22);  new base A (3/2);  total (15/24)
Fisher's exact test for biased strand distribution p-value = 3.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCTCTGGTGCGCTTTCCTGGATGAAATCGGCGTTGCGCACGCATTCCTCGATGGTTGCCACGAACTTCAGGCGGTCTTGCGATGCGCCAGGGGCCAGGCCTTGTTTTTCCAGTGCCGGCCAGGCGTTGGCGACACGTTTGCGCAGTGCCTGTTCGGCGCCGGGCGCCGGGTCCCAGGCCACCACATCCAGGCCATGGGCCAGGGCGCGGGCGACCCAGCCGCTGCCGATCACGCCGCTACCCAGGGCGGCGAAGGTTTTGATT  >  KT2440/356251‑356515
                                                                                                                                    |                                                                                                                                    
gcgcTCTGGTGCGCTTTCCTGGATGAAATCGGCGTTGCGCACGCATTCCTCGATGGTTGCCACGAACTTCAGGCGGTCTTGCGATGCGCCAGGGGCCAGGCCTTGTTTTTCCAGTGCCGGCCAGGCGTTGGCGac                                                                                                                                    >  6:135643/1‑135 (MQ=255)
gcgcTCTGGTGCGCTTTCCTGGATGAAATCGGCGTTGCGCACGCATTCCTCGATGGTTGCCACGAACTTCAGGCGGTCTTGCGATGCGCCAGGGGCCAGGCCTTGTTTTTCCAGTGCCGGCCAGGCGTTGGCGac                                                                                                                                    >  6:211563/1‑135 (MQ=255)
 cgcTCTGGTGCGCTTTCCTGGATGAAATCGGCGTTGCGCACGCATTCCTCGATGGTTGCCACGAACTTCAGGCGGTCTTGCGATGCGCCAGGGGCCAGGCCTTGTTTTTCCAGTGCCGGCCAGGCGTTGGCGaca                                                                                                                                   <  2:296207/135‑1 (MQ=255)
         tGCGCTTTCCTGGATGAAATCGGCGTTGCGCACGCATTCCTCGATGGTTGCCACGAACTTCAGGCGGTCTTGCGATGCGCCAGGGGCCAGGCCTTGTTTTTCCAGTGCCGGCCAGGCGTTGGCGACACGTTTg                                                                                                                             <  7:25139/133‑1 (MQ=255)
            gcTTTCCTGGATGAAATCGGCGTTGCGCACGCATTCCTCGATGGTTGCCACGAACTTCAGGCGGTCTTGCGATGCGCCAGGGGCCAGGCCTTGTTTTTCCAGTGCCGGCCAGGCGTTGGCGACACGTTTGCGCAg                                                                                                                        <  7:19262/135‑1 (MQ=255)
                   tGGATGAAATCGGCGTTGCGCACGCATTCCTCGATGGTTGCCACGAACTTCAGGCGGTCTTGCGATGCGCCAGGGGCCAGGCCTTGTTTTTCCAGTGCCGGCCAGGCGTTGGCGACACGTTTGCGCAGTGCCTGt                                                                                                                 >  7:156272/1‑135 (MQ=255)
                   tGGATGAAATCAGCGTTGCGCACGCATTCCTCGATGGTTGCCACGAACTTCAGGCGGTCTTGCGATGCGCCAGGGGCCAGGCCTTGTTTTTCCAGTGCCGGCCAGGCGTTGGCGACACGTTTGCGCAGTGCCTGt                                                                                                                 <  1:223288/135‑1 (MQ=255)
                             cGGCGTTGCGCACGCATTCCTCGATGGTTGCCACGAACTTCAGGCGGTCTTGCGATGCGCCAGGGGCCAGGCCTTGTTTTTCCAGTGCCGGCCAGGCGTTGGCGACACGTTTGCGCAGTGCCt                                                                                                                   <  5:193244/123‑1 (MQ=255)
                             cGGCGTTGCGCACGCATTCCTCGATGGTGGCCACGAACTTCAGGCGGTCTTGCGATGCGCCAGGGGCCAGGCCTTGTTTTTCCAGTGCCGGCCAGGCGTTGGCGACACGTTTGCGCAGTGCCt                                                                                                                   >  6:193244/1‑123 (MQ=255)
                                cGTTGCGCACGCATTCCTCGATGGTTGCCACGAACTTCAGGCGGTCTTGCGATGCGCCAGGGGCCAGGCCTTGTTTTTCCAGTGCCGGCCAGGCGTTGGCAACACGTTTGCGCAGTGCCTGTTCggcgccgggcg                                                                                                    <  1:269429/135‑1 (MQ=255)
                                  ttGCGCACGCATTCCTCGATGGTTGACACGAACTTCAGGCGGTCTTGCGATGCGCCAGGGGCCAGGCCTTGTTTTTCCAGTGCCGGCCAGGCGTTGGCGACACGTTTGCGCAGTGCCTGTTCggcgccgggcgcc                                                                                                  >  1:107783/1‑135 (MQ=255)
                                     cgcacgcaGTCCTCGCAGGTTGCCACGAACTTCAGGCGGTCTTGCGAGCGGCCAGGGGCCAGGCCTTGTTTTTCCAGTGCGGGCCAGCCGTGGGCGACACGTCTGCGCAGTGCCTGTGCGGcgccgggcgccggg                                                                                               <  6:183120/135‑1 (MQ=255)
                                        acgcaTTCCTCGATGGTTGCCACGAACTTCAGGCGGTCTTGCGATGCGCCAGGGGCCAGGCCTTGTTTTTCCAGTGCCGGCCAGGCGTTGGCGACACGTTTGCGCAGTGCCTGTTCGGCGCCGGGCGCCGGGTcc                                                                                            >  8:192027/1‑135 (MQ=255)
                                         cgcaTTCCTCGATGGGTGCCACGAACTTCAGGCGGTCTTGCGATGCGCCAGGGGCCAGGCCTTGTTTTTCCAGTGCCGGCCAGGCGTTGGCGACACGTTTGCGCAGTGCCTGTTCGGCGCCGGGCGCCGGGTccc                                                                                           <  7:60876/135‑1 (MQ=255)
                                             ttCCTCGATGGTTGCCACGAACTTCAGGCGGTCTTGCGATGCGCCAGGGGCCAGGCCTTGTTTTTCCAGTGCCGGCCAGGCGTTGGCGACACGTTTGCGCAGTGCCTGTTCGGCGCCGGGCGCCGGGTCCCAGGc                                                                                       <  7:259896/135‑1 (MQ=255)
                                             ttCCTCGATGGTTGCCACGAACTTCAGGCGGTCTTGCGATGCGCCAGGGGCCAGGCCTTGTTTATCCAGTGCCGGCCAGGCGTTGGCGACACGTTTGCGCAGTGCCTGTTCGGCGCCGGGCGCCGGGTCCCAGGc                                                                                       <  3:284878/135‑1 (MQ=255)
                                              tCCTCGATGGTTGCCACGAACTTCAGGCGGTCTTGCGATGCGCCAGGGGCCAGGCCTTGTTTTTCCAGTGCCGGCCAGGCGTTGGCGACACGTTTGCGCAGTGCCTGTTCGGCGCCGGGCGCCGGGTCCCAGGcc                                                                                      <  6:106903/135‑1 (MQ=255)
                                                    aTGGTTGCCACGAACTTCAGGCGGTCTTGCGATGCGCCAGGGGCCAGGCCTTGTTTTTCCAGTGCCGGCCAGGCGTTGGCGACACGTTTGCGCAGTGCCTGTTCGGCGCCGGGCGCCGGGTCCCAGGCcaccac                                                                                 <  1:186011/134‑1 (MQ=255)
                                                     tGGTTGCCACGAACTTCAGGCGGTCTTGCGATGCGCCAGGGGCCAGGCCTTGTTTTTCCAGTGCCGGCCAGGCGTTGGCGACACGTTTGCGCAGTGCCTGTTCGGCGCCGGGCGCCGGGTCCCAGGCCACCACAt                                                                               <  1:162637/135‑1 (MQ=255)
                                                     tGGTTGCCACGAACTTCAGGCGGTCTTGCGATGCGCCAGGGGCCAGGCCTTGTTTTTCCAGTGCCGGCCAGGCGTTGGCGACACGTTTGCGCAGTGCCTGTTCGGCGCCGGGCGCCGGGTCCCAGGCCACCACAt                                                                               <  2:107783/135‑1 (MQ=255)
                                                         tGCCACGAACTTCAGGCGGTCTTGCGATGCGCCAGGGGCCAGGCCTTGTTTTTCCAGTGCCGGCCAGGCGTTGGCGACACGTTTGCGCAGTGCCTGTTCGGCGCCGGGCGCCGGGTCCCAGGCCACCACATCCAg                                                                           <  7:50176/135‑1 (MQ=255)
                                                               gAACTTCAGGCGGGCTTGCGATGCGCCAGGGGCCAGGCCTTGTTTTTCCAGTGCCGGCCAGGCGTTGGCGACACGTTTGCGCAGTGCCTGTTCGGCGCCGGGCGCCGGGTCCCAGGCCACCACATCCAGGCCATg                                                                     <  8:46530/135‑1 (MQ=255)
                                                                                         cAGGGGCCAGGCCTTGTTTTTCCAGTGCCGGCCAGGCGTTGGCGACACGTTTGCGCAGTGCCTGTTCGGCGCCGGGCGCCGGGTCCCAGGCCACCACATCCAGGCCATGGGCCAGGGCGCGGGCGACCCAGCCGc                                           >  2:77100/1‑135 (MQ=255)
                                                                                         cAGGGGCCAGGCCTTGTTTTTCCAGTGCCGGCCAGGCGTTGGCAACACGTTTGCGCAGTGCCTGTTCGGCGCCGGGCGCCGGGTCCCAGGCCACCACATCCAGGCCATGGGCCAGGGCGCGGGCGACCCAGCCGc                                           >  6:2716/1‑135 (MQ=255)
                                                                                               ccAGGCCTTGTTTTTCCAGTGCCGGCCAGGCGTGGGCGACACGTTTGCGCAGTGCCTGTGCGGCGCCGGGCGCCGGGTCCCAGGCCACCACATCCAGGCCATGGGCCAGGGCGCGGGCGACCCAGCCGCTGccta                                     <  5:187423/135‑3 (MQ=255)
                                                                                               ccAGGCCTTGTTTTGCCCGGGCCGGCCTGGGGTGGGCGACAGGTTGGCGCTGGCCGTGTGCGGCGCCGGGCGCCGGGTCCCCGGCCACCACATCCCGGCCAGGGGCCAGGGCGCGGGCGACCCAGCCGCTGCCGa                                     <  8:156272/135‑1 (MQ=255)
                                                                                               ccAGGCCTTGTTTTGCCAGTGCCGGCCAGGCGTGGGCGACCCGTTTGCGCAGTGCCTGTGCGGCGCCGGGCGCCGGGTCCCAGGCCACCACATCCAGGCCATGGGCCAGGGCGCGGGCGACCCAGCCGCGGCCGa                                     <  4:53402/135‑1 (MQ=255)
                                                                                                   gCCTTGTTTTTCCAGTGCCGGCCAGGCGTTGGCGACACGTTTGCGCAGTGCCTGTTCggcgccgggc                                                                                                     <  8:283275/67‑1 (MQ=255)
                                                                                                   gCCTTGTTTTTCCAGTGCCGGCCAGGCGTTGGCGACACGTTTGCGCAGTGCCTGTTCggcgccgggc                                                                                                     >  7:283275/1‑67 (MQ=255)
                                                                                                    cgttGTTTTTCCAGTGCCGGCGTGGGGTTGGCGACAGGTTTGCGCAGTGCCTGTTGGGCGCCGGGCGCCGGGTCCCAGGCCACCACATCCAGGCCATGGGCCAGGGCGGGGGCGACCCTGCCGCTGCCGATCAcg                                <  6:221473/133‑1 (MQ=255)
                                                                                                    ccTTGTTTTTCCAGTGCCGGCCAGGCGTTGGCGACACGTTTGCGCAGTGCCTGTTCGGCGCCGGGCGCCGGGTCCCAGGCCACCACATCCAGGCCATGGGCCAGGGCGCGGGCGACCCAGCCCCTGCCGATCAcg                                >  5:221473/1‑135 (MQ=255)
                                                                                                    ccTTGTTTTTCCAGTGCCGGCCAGGCGTTGGCAACACGTTTGCGCAGTGCCTGTGCGGCGCCGGGCGCCGGGTCCCAGGCCACCACATCCAGGCCATGGGCCAGGGCGCGGGCGACCCAGCCGCTGCCGATCAcg                                <  5:53543/135‑1 (MQ=255)
                                                                                                      ttGTTTTTCCAGTGCCGGCCAGGCGTTGGCAACACGTTTGCGCAGTGCCTGTTCGGCGCCGGGCGCCGGGTCCCAGGCCACCACATCCAGGCCATGGGCCAGGGCGCGGGCGACCCAGCCGCTGCCGATCAcgcc                              >  5:25814/1‑135 (MQ=255)
                                                                                                      ttGTTTTTCCAGTGCCGGCCAGGCGTTGGCAACACGTTTGCGCAGTGACTGTTCGGCGCCGGGCGCCGGGTCCCAGGCCACCACATCCAGGCCATGGGCCAGGGCGCGGGCGACCCAGCCGCTGCCGCACAcgcc                              >  5:174201/1‑135 (MQ=255)
                                                                                                                aGTGCCGGCCAGGCGTTGGCGACACGTTTGCGCAGTGCCTGTTCGGCGCCGGGCGCCGGGTCCCAGGCCACCACATCCAGGCCATGGGCCAGGGCGCGGGCGCCCCAGCCGCTGCCGATCACGCCGCGACCCAgg                    >  1:256889/1‑135 (MQ=255)
                                                                                                                     cGGCCAGGCGTTGGCGACACGTTTGCGCAGTGCCTGTTCGGCGCCGGGCGCCGGGTCCCAGGCCACCACATCCAGGCCATGGGCCAGGGCGCGGGCGACCCAGCCGCTGCCGATCACGCCGCTACCCAGggcgg                >  2:268309/1‑134 (MQ=255)
                                                                                                                               ttGGCGACACGTTTGCGCAGTGCCTGTTCGGCGCCGGGCGCCGGGTCCCAGGCCACCACATCCAGGCCATGGGCCAGGGCGCGGGCGACCCAGCCGCTGCCGATCACGCCGCTACCCAGGGCGGCGGAGGTTTTg     >  5:57288/1‑135 (MQ=255)
                                                                                                                                  gCGACTTGTTTGCGCAGGGCCGGTGCGGCGCCGGGCGCCGGGTCCCAGGCCACCAAATCCCGGCCATGGGCCAGGGCGCGGGCTACCCAGCGGCTGCCGATCACGCCGCTACCCAGGGCGGCGAAGGTTTTGAtt  <  5:211563/135‑1 (MQ=255)
                                                                                                                                   cGACAAGTTTGCGCAGTGCCTGTGCGGCGCCGGGCGCCGGGTCCCAGGCCACCACAGCCAGGCCATGGGCCAGGGCGCGGGCGACCCAGCCGCTGCCGATCACGCCGCTACCCAGGGCGGCGAAGGTTTTGAtt  <  8:247967/134‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                    
GCGCTCTGGTGCGCTTTCCTGGATGAAATCGGCGTTGCGCACGCATTCCTCGATGGTTGCCACGAACTTCAGGCGGTCTTGCGATGCGCCAGGGGCCAGGCCTTGTTTTTCCAGTGCCGGCCAGGCGTTGGCGACACGTTTGCGCAGTGCCTGTTCGGCGCCGGGCGCCGGGTCCCAGGCCACCACATCCAGGCCATGGGCCAGGGCGCGGGCGACCCAGCCGCTGCCGATCACGCCGCTACCCAGGGCGGCGAAGGTTTTGATT  >  KT2440/356251‑356515

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