Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 382,479 T→G 36.3% L12R (CTG→CGG)  soxD → Sarcosine oxidase subunit delta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB122382,4790TG36.3% 36.6 / 5.7 51L12R (CTG→CGG) soxDSarcosine oxidase subunit delta
Reads supporting (aligned to +/- strand):  ref base T (11/17);  new base G (6/14);  total (20/31)
Fisher's exact test for biased strand distribution p-value = 5.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCCACTGGCCGCGCCGTTCTCCATGGACCGTTTCTACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATA‑TTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTA‑CATGTTCTACCGTGACAACCCGCGCGGTATTCACCATGA  >  ACB122/382344‑382611
                                                                                                                                        |                                                                                                                                     
cATCCACTGGCCGCGCCGTTCTCCATGGACCGTTTCTACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATA‑TTTTCTGTCCCCACTGCGGCGAGCt                                                                                                                                       <  8:51270/136‑1 (MQ=255)
cATCCACTGGCCGCGCCGTTCTCCATGGACCGTTTCTACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATA‑TTTTCTGTCCCCACTGCGGCGAGCt                                                                                                                                       <  7:249083/136‑1 (MQ=255)
         gCCGCGCCTTTCTCCATGGACCGTTTCTACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATA‑TTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCg                                                                                                                               <  8:332178/135‑1 (MQ=255)
          ccGCGCCGTTCTCCATGGACCGTTGCTACAACGGCGCACTGATCGACTAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATA‑TTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCgaa                                                                                                                             <  7:3815/136‑1 (MQ=255)
             cgcCGTTCTCCATGGACCGTTTCTACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGCCACCGTCATGTTGCATA‑TTTTCTGTCCCCACTGCGGCGGGCGGCGCGCCgca                                                                                                                             >  5:391427/1‑131 (MQ=255)
              gcCGTTCTCCATGGACCGTTTCTACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATA‑TTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCgaa                                                                                                                             >  2:110826/1‑132 (MQ=255)
                cGTTCTCCATGGACCGTTTCTACAACGGCGCACTGATCGACGAACCCGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGATGCATA‑TTTTCTGTCCCCACTGCGGCGCGCAGCGCGCCccagaaca                                                                                                                        >  6:17421/1‑133 (MQ=255)
                      ccATGGACCGTTTCTACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATA‑TTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCGAAGAAGAGtt                                                                                                                     <  7:270275/132‑1 (MQ=255)
                      ccATGGACCGTTTCTACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATA‑TTTTCTGTCCCCACTGCGGCGAGCTGCGCGCCGAAGAAGAGtt                                                                                                                     >  8:270275/1‑132 (MQ=255)
                                tttCTACACCGGCGAACTGACCGACGAACCCGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATA‑TTTTCTGGCCCCACTGCGGCGAGCTGCGCTCCGAAGAAGAGTTCCACGCCTCTGGcc                                                                                                       <  7:46374/136‑1 (MQ=255)
                                   cTACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCTTCGCCCACTAACCGGAGACACCGTCATGTTGCATA‑TTTTCTGTACCACCTTCGGCGGGCAGCGCCCCGCAGAAAACTTCAACGCCTCTGGcccag                                                                                                    >  8:235507/1‑133 (MQ=255)
                                    tACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATA‑TTTTCTGTCCCCACTGCGGCGGGCTGCGCGCCGAAGAAGCGGTGCACGCCTCAGGGCCGgc                                                                                                   >  3:295228/1‑136 (MQ=255)
                                    tACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATA‑TTTTCTGTCCCCACTGCGGCGATCGGCGGTCCGCAGAAGAGTTACACGCCTCAGGCCCGgc                                                                                                   >  3:268874/1‑136 (MQ=255)
                                           gcgcACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCCCTAACCGGAGACACCGTCATGTTGCATA‑TTTTCTGTCCCCACTGCGGCGCGCTGCGCTCCCAATAAGAGTTGCACGCCTCGGGCCCggtctcgt                                                                                              >  6:70683/1‑128 (MQ=255)
                                               aCTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATA‑TTTTCTGTCCCCACTGCGGCGCGCGGGGCGCCGCAGAAGAGTTCCACCCCGCGGGCCCGGCGCGCATTgcca                                                                                        >  7:367795/1‑135 (MQ=255)
                                                   aTCGACGAACACGGCGCCGCCGCCGGCGCCCACTAACCGGAGACACCGTCATGTTGAATA‑TTTTCTGGCCCCACTGCGGCGAGCGGCGCTCCGAAGAAGCGTGCCACGCCTCTGGCCAGGCGCACATTgcccgccc                                                                                    <  4:339941/136‑1 (MQ=255)
                                                    tccccgaACCCGGCGCCGCCGCCTTCGCCCACTAACCGGAGACCCCGTCATGTTGAATA‑TTTTCTGTCCCCACTGCGGCGCGCTGCGCTCCGAAGAAGAGTGCCACGCCTCTGGCCAGGCGCACATTgcccgccc                                                                                    <  7:197541/131‑1 (MQ=255)
                                                      gacgaACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATA‑TTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCGCACATTGCccgcccgc                                                                                  >  2:78248/1‑135 (MQ=255)
                                                      gacgaACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATA‑TTTTCTGTCCCCACTGCGGCGAGCGGGGCTCCGAAGCAGAGTTCCCCGCCGCTGGCCCGGCGCCCCTTGCccgcccgc                                                                                  >  6:189015/1‑135 (MQ=255)
                                                          accccGGCGCCGCCGCCGTCGCCCACTAACCGCAGACACCGTCATGTTCCATA‑TTTTCTGCCCCCACTGCGGCGAGCGGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGAc                                                                             <  7:92167/132‑1 (MQ=255)
                                                          aaCACGGCGCCGCCGCCGTCGCCCAGTAACCGGAGACACCGTCATGTTGCATT‑TTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCTAAGAAGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGAc                                                                             <  3:186867/136‑1 (MQ=255)
                                                            cacGGCGCCGCCGCCGTCGCCCAATAACCGGACACACCGTCTTGTTGCAAA‑TTTTCTGTCCCCACTGCGGCGAGCTGCGGTCCGAACACGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGAc                                                                             <  6:164254/134‑1 (MQ=255)
                                                               ggCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCTTGTTGCATA‑TTTTCTGTCCCCACTGCGGCGAGCTGCGCGCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTaa                                                                        <  8:241409/136‑1 (MQ=255)
                                                                gcgcCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATA‑TTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAAc                                                                       >  8:308363/1‑136 (MQ=255)
                                                                     gccgccgTCGCCCACTAACCGGAGACACCGTCATGTTGCATA‑TTTTCTGTCCCCACTGCGGCGAGCGGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCt                                                                   <  1:78248/135‑1 (MQ=255)
                                                                     gccgccgTCGCCCACTAACCGGAGACACCGTAATGTTGCATA‑TTTTCTGTCCCCACTGCGGAGAGCTGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCt                                                                   <  6:154890/135‑1 (MQ=255)
                                                                     gccgccgTCGCCAAATAACGGGAGACCCCGTCATTTTGATTA‑TTTGCTTCCCCAACTGCGGCGCGCGGGGCTCGGATGATGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCt                                                                   <  2:325715/135‑1 (MQ=255)
                                                                         ccgGCGCCCACTACCCGGAGACCCGGTCATGTTGCATA‑TTTTCTGCCCCCACTGCGGCGAGCGGCGCGCCGACGAAGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTcc                                                              <  7:2855/136‑1 (MQ=255)
                                                                           gacGCCAACAAACCGGAGACCCCGGAATGTTCAATA‑TTTTCTGCCCCCACTGCGGCGAGCGGCGCTCCGAAGATTAGTGCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCcg                                                             <  2:102959/133‑1 (MQ=255)
                                                                           gTCCCCCACTAACCGGCGACACCGTCTTGTTG‑ATATTTTTCTGCCCCCCCTGCGGCTGGCGGCGCCCCGAAGAAGAGTGCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCcg                                                             <  6:220279/135‑1 (MQ=255)
                                                                            tCGCCCACTACCCGGACACACCGTCAGTGTGAATA‑TTTTCTGTCCCCCCTGCGGCGAGCGGCGCTCCGAAGCAGAGGTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCcgac                                                           <  7:135053/136‑1 (MQ=255)
                                                                              gCCCACTAACCGGAGACACCGTCATGTTGCATA‑TTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCGAAGAAGAGtt                                                                                                                     >  8:257057/1‑76 (MQ=255)
                                                                              gCCCACTAACCGGAGACACCGTCATGTTGCATA‑TTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCGAAGAAGAGtt                                                                                                                     <  7:257057/76‑1 (MQ=255)
                                                                                            gaCACCGTCATGTTGCATA‑TTTTCTGCCCCCCCTGCGGCGAGCTGCGCGCCGAAGACGAGTCCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTAcc                                            <  6:64630/135‑1 (MQ=255)
                                                                                            cccacCGTCATGGTTCATA‑TTTTCTGTCCCCACGGCGGCGAGCGGCGCTCCGAAGAAGAGTCCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTAcc                                            <  5:242980/133‑1 (MQ=255)
                                                                                              cacGGTCATGTGGAATA‑TTTTCTGTCCCCTCTGCGGCGAGGTGCGCTCCGAAGAAGAGTACCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTA‑c                                        <  4:174261/136‑1 (MQ=255)
                                                                                                ccGTCATGTTGCATA‑TTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGCCCCAAACGCCTGCTCCGCCGAGGAGTGGGGTACCTA‑CAt                                      >  8:208238/1‑136 (MQ=255)
                                                                                                    cATGTTGCATA‑TTTTCTGCCCCCACTGCGGCGCGCGGCGCTCCGAACAAGAGTGCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTA‑CATGtt                                   <  2:220277/135‑1 (MQ=255)
                                                                                                           cATA‑TTTTCTGGCCCCCCTGCGGCGAGCTGCGCTCCGAAGAAGATTGCCACGCCTCTGGCCAGGCGCACTTTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGTGTACCTA‑CATGTTCTACCGt                            <  8:284009/135‑1 (MQ=255)
                                                                                                            ata‑tTTTCTGTCCCCACTGCGGCGCGCGGCGCGCCGAAGAAGAGTTCCACGCCTCTCGCCACGCGCACATTGCCCCCCCGCGGGCCCCGAACGCCGGCTCCGACGCGGAGAGGGGTACCGA‑CATGTTCTACCGTg                           >  5:311532/1‑135 (MQ=255)
                                                                                                              a‑tTTTCTGTCCCCACTGCGGCGCGCTGCCCGCCCAACAACACTTACACCCCGCAGGCCCCGCGCACAATGCCCGCCCCCAGGCCACAAACGCCAGCGCCCCCGAGGAGGGGGTGACCAACCATG‑TCTACCCTGACa                        >  4:86646/1‑136 (MQ=255)
                                                                                                              a‑tTTTCTGTCCCCACTGCGGCGAGCAGCGCGCCGAAGAAGAGTTCCACCCCTCTGGCCAGGCGCCCATTTCCCCCCCCCGGGCCCCGAACGCCGGCTCCCACGAGGCGTTGGGTACCTA‑CATTTTCTAACGTgg                          >  2:341299/1‑133 (MQ=255)
                                                                                                                    cTGACCCCACTGCGGAGAGCTGGGCTCCGAAGAAGAGTCCCAGGGCTGTGGCCAGGCGCACATTGGCCGCCCGCTGGACCCTAACGCCTGCTCCGTCGAGGAGTGGGGTACCTA‑CATGTTCTACCGTGACAAccg                    <  3:155062/135‑2 (MQ=255)
                                                                                                                     tgcccccccGGCGGCGTGCGGCCCTCGGAACACGCGTACCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTA‑CATGTTCTACCGTGACAAccc                    <  1:341299/131‑1 (MQ=255)
                                                                                                                     tgcccccACTGCGGCGAGCGGCGCGCCGAATAAGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTA‑CATGTTCTACCGTGACAAccc                    <  1:110826/131‑1 (MQ=255)
                                                                                                                     tGTCCCCACTGCGGCGAGCTGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTA‑CATGTTCTACCGTGACAAccc                    <  1:40575/134‑1 (MQ=255)
                                                                                                                      gTCCCCACTGCGGCGGGCTGCGCTCCGAAGAAGAGTTGCACGCCTCTGGCCAGGCGCCCATTGCCCGCCCGCTGGGCCCGAACGCCTGCTCCGGCGAGGGGTGGGGTAGCTA‑CATGTTGTACCGTGGCACCCcgc                  >  2:91052/1‑135 (MQ=255)
                                                                                                                         ccgcgtGCGGCTCGCGGCGCTCCGACTTAGTGCGCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTA‑CATGTTCTACCGTGACATCCCGCGCgg              <  7:337293/131‑1 (MQ=255)
                                                                                                                         ccccCTGCGGCGAGCTGCGCTCCGAAGAAGCGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTA‑CATGTTCTACCGTGACAACCCGCGCgg              <  8:67275/136‑1 (MQ=255)
                                                                                                                          ccACTGCGGCGAGCGGCGCTCCGACGACGCGCCCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTA‑CATGTTCTACCGTGACAACCCGCGCgg              <  5:259354/135‑1 (MQ=255)
                                                                                                                            aCTGCGGCGCGCGGCGCTCCGTAGAAGACTTCCCCGCCTCTGGCCAGGCGCACATTGCCCGCCCGCGGGCCCCTAACCCCTGCTCCGACGAGGAGTGGGGTTCCTA‑CATGTTCTACCGTGACCACCCGCGCggga            >  1:50847/1‑133 (MQ=255)
                                                                                                                                 ggcgAGCTGCGCTCCGAAGAAGAGTGCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTA‑CATGTTCTACCGTGACAACCCGCGCGGTATTCAcc      <  7:208238/136‑1 (MQ=255)
                                                                                                                                      gCTGCGCTCCGAAGAAGAGTTACACGCCTCTGGCCCGGCGCACATTTCCCGCCCGCTGGGCCCTAACGCCTGCTCCGCCGAGGAGTGGGGTACCGA‑CATGTTCTACCGTGCCAACCCGCGCGGGATTCACCATGa  >  8:254453/1‑135 (MQ=255)
                                                                                                                                        |                                                                                                                                     
CATCCACTGGCCGCGCCGTTCTCCATGGACCGTTTCTACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATA‑TTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTA‑CATGTTCTACCGTGACAACCCGCGCGGTATTCACCATGA  >  ACB122/382344‑382611

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