Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 710,940 A→G 8.7% V186V (GTA→GTG)  thiO → Glycine oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT2440710,9400AG 8.7% 135.3 / 3.4 46V186V (GTA→GTG) thiOGlycine oxidase
Reads supporting (aligned to +/- strand):  ref base A (16/26);  new base G (2/2);  total (18/28)
Fisher's exact test for biased strand distribution p-value = 6.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

TGGCCAATGTGCGCAACCCACGCTTGGTCAAGTCGCTGAAGGCGGCGTTGCTGGCGCTGCCGAACGTGAGCTTGCGCGAACACTGCCAGATCACCGGCTTCGTGCAGGAAAAGGGGCGCGTGACTGGTGTACAGACGGCTGATGGCGTACTGGCGGCGGATGAGGTTGTGCTCAGCGCAGGCGCCTGGAGTGGC‑GACCTGTTGCGTACCCTTGGCCTTGAACTGCCGGTGGAGCCCGTGAAGGGCCAGATGATTCTGTTCA  >  KT2440/710810‑711070
                                                                                                                                  |                                                                                                                                   
tgGCCAATGTGCGCAACCCACGCTTGGTCAAGTCGCTGAAGGCGGCGTTGCTGGCGCTGCCGAACGTGAGCTTGCGCGAACACTGCCAGATCACCGGCTTCGTGCAGGAAAAGGGGCGCGTGACTGGTGTGCAGa                                                                                                                                 >  6:183207/1‑135 (MQ=255)
tgGCCAATGTGCGCAACCCACGCTTGGTCAAGTCGCTGAAGGCGGCGTTGCTGGCGCTGCCGAACGTGAGCTTGCGCGAACACTGCCAGATCACCGGCTTCGTGCAGGAAAAGGGGCGCGTGACTGGTGTACAgc                                                                                                                                 >  5:26608/1‑134 (MQ=255)
tgGCCAATGTGCGCAACCCACGCTTGGTCAAGTCGCTGAAGGCGGCGTTGCTGGCGCTGCCGAACGTGAGCTTGCGCGAACACTGCCAGATCACCGGCTTCGTGCAGGAAAAGGGGCGCGTGACTGGTGTACAGa                                                                                                                                 >  4:146700/1‑135 (MQ=255)
tgGCCAATGTGCGCAACCCACGCTTGGTCAAGTCGCTGAAGGCGGCGTTGCTGGCGCTGCCGAACGTGAGCTTGCGCGAACACTGCCAGATCACCGGCTTCGGGCAGGAAAAGGGGCGCGTGACTGGTGTGCAGa                                                                                                                                 >  6:183209/1‑135 (MQ=255)
tgGCCAATGTGCGCAACCCACGCTTGGTCAAGTCGCTGAAGGCGGCGTTGCTGGCGCTGCCGAACGAGAGCTTGCGCGAACACTGCCAGATCCCCGGCTTCGTGCAGGAAAAGGGGCGCCTCACTGGTGTACAGa                                                                                                                                 >  3:194893/1‑135 (MQ=255)
     aaTGTGCGCAACCCACGCTTGGTCAAGTCGCTGAAGGCGGCGTTGCTGGCGCTGCCGAACGTGAGCTTGCGCGAACACTGCCAGATCACCGGCTTCGTGCAGGAAAAGGGGCGCGTGACTGGTGTACAGACGGc                                                                                                                             >  3:128882/1‑134 (MQ=255)
                     gCTTGGTCAAGTCGCTGAAGGCGGCGTTGCTGGCGCTGCCGAACGTGAGCTTGCGCGAACACTGCCAGATCACCGGCTTCGTTCAGGAAAAGGGGCGCGTGACTGGTGTACAGACGGCTGATGGCGTACTggcggc                                                                                                           <  8:87178/136‑1 (MQ=255)
                      cTTGGTCAAGTCGCTGAAGGCGGCGTTGCTGGCGCTGCCGAACGTGAGCTTGCGCGAACACTGCCAGATCACCGGCTTCGTGCAGGAAAAGGGGCGCGTGACTGGTGTACAGACGGCTGATGGCGTACTggcggc                                                                                                           >  6:97463/1‑135 (MQ=255)
                      cTTGGTCAAGTCGCTGAAGGCGGCGTTGCTGGCGCTGCCGAACGTGAGCTTGCGCGAACACTGCCAGATCACCGGCTTCGTGCAGGAAAAGGGGCGCGTGACTGGTGTACAGACGGCTGATGGCGTACTggcggc                                                                                                           >  6:102114/1‑135 (MQ=255)
                            cAAGTCGCTGAAGGCGGCGTTGCTGGCGCTGCCGAACGTGAGCTTGCGCGAACACTGCCAGATCACCGGCTTCGTGCAGGAAAAGGGGCGCGTGACTGGTGTACAGACGGCTGATGGCGTACTGGCGGCGGATGAg                                                                                                    <  7:99332/136‑1 (MQ=255)
                             aaGTCGCTGAAGGCGGCGTTGCTGGCGCTGCCGAACGTGAGCTTGCGCGAACACTGCCAGATCACCGGCTTCGTGCAGGAAAAGGGGCGCGTGACTGGTGTACAGACGGCTGATGGCGTACTGGCGGCGGATGAg                                                                                                    >  2:176667/1‑135 (MQ=255)
                             aaGTCGCTGAAGGCGGCGTTGCTGGCGCTGCCGAACGTGAGCTTGCGCGAACACTGCCAGATCACCGGCTTCGTGCAGGAAAAGGGGCGCGTGACTGGTGTACAGACGGCTGATGGCGTACTGGCGGCGGATGAg                                                                                                    >  2:59589/1‑135 (MQ=255)
                             aaGTCGCGGAAGGCGGCGTTGCTGGCGCTGCCGAACGTGAGCTTGCGCGAACACTGCCAGATCACCGGCTTCGTGCAGGAAAAGGGGCGCGTGACTGGTGTACAGACGGCTGATGGCGTACTGGCGGCGGATGAg                                                                                                    >  2:182657/1‑135 (MQ=255)
                                  gCTGAAGGCGGCGTTGCTGGCGCTGCCGAACGTGAGCTTGCGCGAACACTGCCAGATCACCGGCTTCGTGCAGGAAAAGGGGCGCGTGACTGGTGTACAGACGGCTGATGGCGTACTGGCGGCGGATGAGGTtgt                                                                                               >  1:189626/1‑135 (MQ=255)
                                         gcggcgTTGCTGGCGCTGCCGAACGTGAGCTTGCGCGAACACTGCCAGATCACCGGCTTCGTGCAGGAAAAGGGGCGCGTGACTGGTGTACAGACGGCTGATGGCGTACTGGCGGCGGATGAGGTTGTGCTCAgc                                                                                        <  1:42785/135‑1 (MQ=255)
                                                      cgcTGCCGAACGTGAGCTTGCGCGAACACTGCCAGATCACCGGCTTCGTGCAGGAAAAGGGGCGCGTGACTGGTGTACAGACGGCTGATGGCGTACTGGCGGCGGATGAGGTTGTGCTCAGCGCAGGCGCCTGGa                                                                           <  4:90071/135‑1 (MQ=255)
                                                           ccGAACGTGAGCTTGCGCGAACACTGCCAGATCACCGGCTTCGTGCAGGAAAAGGGGCGCGTGACTGGTGTACAGACGGCTGATGGCGTACTGGCGGCGGATGAGGTTGTGCTCAGCGCAGGCGCCTGGAGTGGc                                                                      <  5:85156/135‑1 (MQ=255)
                                                            cGAACGTGAGCTTGCGCGAACACTGCCAGATCACCGGCTTCGTGCAGGAAAAGGGGCGCGTGACTGGTGTACAGACGGCTGATGGCGTACTGGCGGCGGATGAGGTTGTGCTCAGCGCAGGCGCCTGGAGTGGc                                                                      <  3:146700/134‑1 (MQ=255)
                                                            cGAACGTGAGCTTGCGCGAACACTGCCAGATCACCGGCTTCGTGCAGGAAAAGGGGCGCGTGACTGGTGTACAGACGGCTGATGGCGTACTGGCGGCGGATGAGGTTGTGCTCAGCGCAGGCGCCTGGAGTGGC‑Ga                                                                   <  8:12024/136‑1 (MQ=255)
                                                             gAACGTGAGCTTGCGCGAACACTGCCAGATCCCCGGCTTCGTGCAGGAAAAGGGGCGCGTGACGGGTGTACAGACGGCTGATGGCGTACTGGCGGCGGATGAGGTTGTGCTCAgc                                                                                        <  1:241213/115‑1 (MQ=255)
                                                             gAACGTGAGCTTGCGCGAACACTGCCAGATCACCGGCTTCGTGCAGGAAAAGGGGCGCGTGACTGGTGTACAGACGGCTGATGGCGTACTGGCGGCGGATGAGGTTGTGCTCAgc                                                                                        >  2:241213/1‑115 (MQ=255)
                                                              aaCGTGAGCTTGCGCGAACACTGCCAGATCACCGGCTTCGTGCAGGAAAAGGGGCGCGTGACTGGTGTACAGACGGCTGATGGCGTACTGGCGGCGGATGAGGTTGTGCTCAGCGCAGGCGCCTGGAGTGGC‑ggc                                                                  >  6:260666/1‑133 (MQ=255)
                                                                   gAGCTTGCGCGAACACTGCCAGATCACCGGCTTCGTGCAGGAAAAGGGGCGCGTGACTGGTGTACAGACGGCTGATGGCGTACTGGCGGCGGATGAGGTTGTGCTCAGCGCAGGCGCCTGGAGTGGC‑GACCTGtt                                                             <  6:225602/135‑1 (MQ=255)
                                                                   gAGCTTGCGCGAACACTGCCAGATCACCGGCTTCGTGCAGGAAAAGGGGCGCGTGACTGGTGTACAGACGGCTGATGGCGTACTGGCGGCGGATGAGGTTGTGCTCAGCGCAGGCGCCTGGAGTGGC‑GACCTGtt                                                             <  1:182657/135‑1 (MQ=255)
                                                                   gAGCTTGCGCGAACACTGCCAGATCACCGGCTTCGTGCAGGAAAAGGGGCGCGTGACTGGTGTACAGACGGCTGATGGCGTACTGGCGGCGGATGAGGTTGTGCTCAGCGCAGGCGCCTGGAGTGGC‑GACCTGtt                                                             <  1:150043/135‑1 (MQ=255)
                                                                   gAGCTTGCGCGAACACTGCCAGATCACCGGCTTCGTGCAGGAAAAGGGGCGCGTGACTGGTGTACAGACGGCTGATGGCGTACTGGCGGCGGATGAGGTTGTGCTCAGCGCAGGCGCCTGGAGTGGC‑GACCTGtt                                                             <  3:274106/135‑1 (MQ=255)
                                                                        tGCGCGAACACTGCCAGATCACCGGCTTCGTGCAGGAAAAGGGGCGCGTGACTGGTGTACAGACGGCTGATGGCGTACTGGCGGCGGATGAGGTTGTGCTCAGCGCAGGCGCCTGGAGTGGC‑GACCTGTTg                                                            >  4:232809/1‑131 (MQ=255)
                                                                        tGCGCGAACACTGCCAGATCACCGGCTTCGTGCAGGAAAAGGGGCGCGTGACTGGTGTACAGACGGCTGATGGCGTACTGGCGGCGGATGAGGTTGTGCTCAGCGCAGGCGCCTGGAGTGGC‑GACCTGTTg                                                            <  3:232809/131‑1 (MQ=255)
                                                                          cgcgAACACTGCCAGATCACCGGCTTCGTGCAGGAAAAGGGGCGCGTGACTGGTGTGCAGACGGCTGATGGCGTACTGGCGGCGGATGAGGTTGTGCTCAGCGCAGGCGCCTGGAGTGGC‑GACCTGTTGCGTAc                                                       <  5:183207/134‑1 (MQ=255)
                                                                          cgcgAACACTGCCAGATCACCGGCTTCGTGCAGGAAAAGGGGCGCGTGACTGGTGTACAGACGGCTGATGGCGTACTGGCGGCGGATGAGGTTGTGCTCAGCGCAGGCGCCTGGAGTGGC‑GACCTGTTGCGTAc                                                       <  1:59589/134‑1 (MQ=255)
                                                                          cgcgAACACTGCCAGATCACCGGCTTCGTGCAGGAAAAGGGGCGCGTGACTGGTGTACAGACGGCTGATGGCGTACTGGCGGCGGATGAGGTTGTGCTCAGCGCAGGCGCCTGGAGTGGC‑GACCTGTTGCGTAc                                                       <  2:189626/134‑1 (MQ=255)
                                                                          cgcgAACACTGCCAGATCACCGGATGCGTGCAGGAAACGGGGCGCGTGACTGGTGTGCAGACGGCTGATGGCGTACTGGCGGCGGATGAGGTTGTGCTCAGCGCAGGCGCCTGGAGTGGC‑GACCTGTTGCGTAc                                                       <  5:183209/134‑1 (MQ=255)
                                                                                  cTGCCAGATCACCGGCTTCGTGCAGGAAAAGGGGCGCGTGACTGGTGTACAGACGGCTGATGGCGTACTGGCGGCGGATGAGGTTGTGCTCAGCGCAGGCGCCTGGAGTGGC‑GACCTGTTGCGTACCCTTGGCCt                                              <  3:11643/135‑1 (MQ=255)
                                                                                       aGATCTCCGGCTTCGTGCAGGAAAAGGGGCTCGTGAGTGGTGTACTGACGGCTGATGGCGTACTGGCGGAGGATGAGGTTGTGCTCAGCGCAGGCGGCTGGAGTGGC‑GAGCTGTTGCGTACCCTTGGCCTTGAAc                                         <  4:194893/135‑1 (MQ=255)
                                                                                                         aGGAACAGGGGCGCGGGACTGGTGTACAGACGGCTGATGGCGTACTGGCGGCGGATGAGGTTGTGCTCAGCGCAGGCGCCTGGAGTGGC‑GACCTGTTGCGTACCCTTGGCCTTGAACTGCCGGTGGAGCCCGTGaa                      <  7:237009/136‑1 (MQ=255)
                                                                                                         aGGAAACGGGGCGCGTGACTGGTGTACAGACGGCTGATGGCGTACTGGCGGCGGATGAGGTGGTGCTCAGCGCAGGCGCCTGGAGTGGC‑GACCTGTTGCGTACCCTTGGCCTTGAACTGCCGGTGGAGCCCGTGaa                      <  8:102350/136‑1 (MQ=255)
                                                                                                              aaGGGGCGCGTGACTGGTGTACAGACGGCTGATGGCGTACTGGCGGCGGATGAGGTTGTGCTCAGCGCAGGCGCCTGGAGTGGC‑GACCTGTTGCGTACCCTTGGCCTTGAACTGCCGGTGGAGCCCGTGAAggg                   >  2:196074/1‑134 (MQ=255)
                                                                                                                  ggCGCGTGACTGGTGTACAGACGGCTGATGGCGTACTGGCGGCGGATGAGGTTGTGCTCAGCGCAGGCGCCTGGAGTGGC‑GACCTGTTGCGTACCCTTGGCCTTGAACTGCCGGTGGAGCCCGTGAAGGGCCAga              <  5:115831/135‑1 (MQ=255)
                                                                                                                  ggCGCGTGACTGGTGTACAGACGGCTGATGGCGTACTGGCGGCGGATGAGGTTGTGCTCAGCGCAGGCGCCTGGAGTGGC‑GACCTGTTGCGTACCCTTGGCCTTGAACTGCCGGTGGAGCCCGTGAAGGGCCAga              <  4:250968/135‑1 (MQ=255)
                                                                                                                  ggCGCGTGACTGGTGTACAGACGGCTGATGGCGTACTGGCGGCGGATGAGGTTGTGCTCAGCGCAGGCGCCTGGAGTGGC‑GACCTGTTGCGTACCCTTGGCCTTAAACTGCCGGTGGAGCCCGTGAAGGGCCAga              <  1:29012/135‑1 (MQ=255)
                                                                                                                   gcgcGTGACTGGTGTACAGACGGCTGATGGCGTACTGGCGGCGGATGAGGTTGTGCTCAGCGCAGGCGCCTGGAGTGGCGGACGTGTTGCGTACCCTTGGCCTTGAAATGCCGGTGGGGCCCGGGAAGGGCCAga              >  1:228330/1‑135 (MQ=255)
                                                                                                                      cgTGACTGGTGGACAGACGGCTGATGGCGTACTGGCGGCGGATGAGGTTGTGCTCAGCGCAGGCGCCTGGAGTGGC‑GACCTGTTGCGTACCCTTGGCCGTGAACTGCCGGTGGAGCCCGTGAAGTGCCAGATGAtt         <  7:126485/136‑1 (MQ=255)
                                                                                                                         gACTGGTGTACAGACGGCTGATGGCGTACTGGCGGCGGATGAGGTTGTGCTCAGCGCAGGCGCCTGGAGTGGC‑GACCTGTTGCGTACCCTTGGCCTTGAACTGCCGGTGGAGCCCGTGAAGGGCCAGATGATTCTg      >  7:50053/1‑136 (MQ=255)
                                                                                                                            tGGTGTACAGACGGCTGATGGCGTACTGGCGGCGGATGAGGTTGTGCTCAGCGCAGGCGCCTGGAGTGGC‑GACCTGTTGCGTACCCTTGGCCTTGAACTGCCGGTGGAGCCCGTGAAGGGCCAGATGATTCTGtt    <  6:268477/135‑1 (MQ=255)
                                                                                                                             ggTGTACAGACGGCTGATGGCGTACTGGCGGCGGATGAGGGTGTGCTCAGCGCAGGCGCCTGGAGTGGC‑GACCTGTTGCGTACCCTTGGCCTTGAACTGCCGGTGGAGCCCGTGAAGGGCCAGATGATTCTGTTCa  <  8:137872/136‑1 (MQ=255)
                                                                                                                               tgtAGATACGGCTGTTTGCGTACTGGCGGGGGATGAGGTTGTGGTCAGCGCAGGCGGCTGTAGTGGC‑GACCTGTTGGGTTGCCTTGGCCTTGAACTGCGGGTGGATCCCGTTATGGGACAGATGATTCTGTTCa  <  2:30432/134‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                   
TGGCCAATGTGCGCAACCCACGCTTGGTCAAGTCGCTGAAGGCGGCGTTGCTGGCGCTGCCGAACGTGAGCTTGCGCGAACACTGCCAGATCACCGGCTTCGTGCAGGAAAAGGGGCGCGTGACTGGTGTACAGACGGCTGATGGCGTACTGGCGGCGGATGAGGTTGTGCTCAGCGCAGGCGCCTGGAGTGGC‑GACCTGTTGCGTACCCTTGGCCTTGAACTGCCGGTGGAGCCCGTGAAGGGCCAGATGATTCTGTTCA  >  KT2440/710810‑711070

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: ‑