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

CATCCACTGGCCGCGCCGTTCTCCATGGACCGTTTCTACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTACATGTTCTACCGTGACAACCCGCGCGGTATTCACCATGA  >  ACB122/382344‑382611
                                                                                                                                        |                                                                                                                                     
cATCCACTGGCCGCGCCGTTCTCCATGGACCGTTTCTACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGAGCt                                                                                                                                       <  8:51270/136‑1 (MQ=255)
cATCCACTGGCCGCGCCGTTCTCCATGGACCGTTTCTACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGAGCt                                                                                                                                       <  7:249083/136‑1 (MQ=255)
         gCCGCGCCTTTCTCCATGGACCGTTTCTACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCg                                                                                                                               <  8:332178/135‑1 (MQ=255)
          ccGCGCCGTTCTCCATGGACCGTTGCTACAACGGCGCACTGATCGACTAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCgaa                                                                                                                             <  7:3815/136‑1 (MQ=255)
             cgcCGTTCTCCATGGACCGTTTCTACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGCCACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGGGCGGCGCGCCgca                                                                                                                             >  5:391427/1‑131 (MQ=255)
              gcCGTTCTCCATGGACCGTTTCTACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCgaa                                                                                                                             >  2:110826/1‑132 (MQ=255)
                cGTTCTCCATGGACCGTTTCTACAACGGCGCACTGATCGACGAACCCGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGATGCATATTTTCTGTCCCCACTGCGGCGCGCAGCGCGCCccagaaca                                                                                                                        >  6:17421/1‑133 (MQ=255)
                      ccATGGACCGTTTCTACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCGAAGAAGAGtt                                                                                                                     <  7:270275/132‑1 (MQ=255)
                      ccATGGACCGTTTCTACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGAGCTGCGCGCCGAAGAAGAGtt                                                                                                                     >  8:270275/1‑132 (MQ=255)
                                tttCTACACCGGCGAACTGACCGACGAACCCGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGGCCCCACTGCGGCGAGCTGCGCTCCGAAGAAGAGTTCCACGCCTCTGGcc                                                                                                       <  7:46374/136‑1 (MQ=255)
                                   cTACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCTTCGCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTACCACCTTCGGCGGGCAGCGCCCCGCAGAAAACTTCAACGCCTCTGGcccag                                                                                                    >  8:235507/1‑133 (MQ=255)
                                    tACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGGGCTGCGCGCCGAAGAAGCGGTGCACGCCTCAGGGCCGgc                                                                                                   >  3:295228/1‑136 (MQ=255)
                                    tACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGATCGGCGGTCCGCAGAAGAGTTACACGCCTCAGGCCCGgc                                                                                                   >  3:268874/1‑136 (MQ=255)
                                           gcgcACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCCCTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGCGCTGCGCTCCCAATAAGAGTTGCACGCCTCGGGCCCggtctcgt                                                                                              >  6:70683/1‑128 (MQ=255)
                                               aCTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGCGCGGGGCGCCGCAGAAGAGTTCCACCCCGCGGGCCCGGCGCGCATTgcca                                                                                        >  7:367795/1‑135 (MQ=255)
                                                   aTCGACGAACACGGCGCCGCCGCCGGCGCCCACTAACCGGAGACACCGTCATGTTGAATATTTTCTGGCCCCACTGCGGCGAGCGGCGCTCCGAAGAAGCGTGCCACGCCTCTGGCCAGGCGCACATTgcccgccc                                                                                    <  4:339941/136‑1 (MQ=255)
                                                    tccccgaACCCGGCGCCGCCGCCTTCGCCCACTAACCGGAGACCCCGTCATGTTGAATATTTTCTGTCCCCACTGCGGCGCGCTGCGCTCCGAAGAAGAGTGCCACGCCTCTGGCCAGGCGCACATTgcccgccc                                                                                    <  7:197541/131‑1 (MQ=255)
                                                      gacgaACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCGCACATTGCccgcccgc                                                                                  >  2:78248/1‑135 (MQ=255)
                                                      gacgaACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGAGCGGGGCTCCGAAGCAGAGTTCCCCGCCGCTGGCCCGGCGCCCCTTGCccgcccgc                                                                                  >  6:189015/1‑135 (MQ=255)
                                                          accccGGCGCCGCCGCCGTCGCCCACTAACCGCAGACACCGTCATGTTCCATATTTTCTGCCCCCACTGCGGCGAGCGGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGAc                                                                             <  7:92167/132‑1 (MQ=255)
                                                          aaCACGGCGCCGCCGCCGTCGCCCAGTAACCGGAGACACCGTCATGTTGCATTTTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCTAAGAAGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGAc                                                                             <  3:186867/136‑1 (MQ=255)
                                                            cacGGCGCCGCCGCCGTCGCCCAATAACCGGACACACCGTCTTGTTGCAAATTTTCTGTCCCCACTGCGGCGAGCTGCGGTCCGAACACGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGAc                                                                             <  6:164254/134‑1 (MQ=255)
                                                               ggCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCTTGTTGCATATTTTCTGTCCCCACTGCGGCGAGCTGCGCGCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTaa                                                                        <  8:241409/136‑1 (MQ=255)
                                                                gcgcCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAAc                                                                       >  8:308363/1‑136 (MQ=255)
                                                                     gccgccgTCGCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGAGCGGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCt                                                                   <  1:78248/135‑1 (MQ=255)
                                                                     gccgccgTCGCCCACTAACCGGAGACACCGTAATGTTGCATATTTTCTGTCCCCACTGCGGAGAGCTGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCt                                                                   <  6:154890/135‑1 (MQ=255)
                                                                     gccgccgTCGCCAAATAACGGGAGACCCCGTCATTTTGATTATTTGCTTCCCCAACTGCGGCGCGCGGGGCTCGGATGATGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCt                                                                   <  2:325715/135‑1 (MQ=255)
                                                                         ccgGCGCCCACTACCCGGAGACCCGGTCATGTTGCATATTTTCTGCCCCCACTGCGGCGAGCGGCGCGCCGACGAAGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTcc                                                              <  7:2855/136‑1 (MQ=255)
                                                                           gacGCCAACAAACCGGAGACCCCGGAATGTTCAATATTTTCTGCCCCCACTGCGGCGAGCGGCGCTCCGAAGATTAGTGCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCcg                                                             <  2:102959/133‑1 (MQ=255)
                                                                           gTCCCCCACTAACCGGCGACACCGTCTTGTTGATATTTTTCTGCCCCCCCTGCGGCTGGCGGCGCCCCGAAGAAGAGTGCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCcg                                                             <  6:220279/135‑1 (MQ=255)
                                                                            tCGCCCACTACCCGGACACACCGTCAGTGTGAATATTTTCTGTCCCCCCTGCGGCGAGCGGCGCTCCGAAGCAGAGGTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCcgac                                                           <  7:135053/136‑1 (MQ=255)
                                                                              gCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCGAAGAAGAGtt                                                                                                                     >  8:257057/1‑76 (MQ=255)
                                                                              gCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCGAAGAAGAGtt                                                                                                                     <  7:257057/76‑1 (MQ=255)
                                                                                            gaCACCGTCATGTTGCATATTTTCTGCCCCCCCTGCGGCGAGCTGCGCGCCGAAGACGAGTCCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTAcc                                            <  6:64630/135‑1 (MQ=255)
                                                                                            cccacCGTCATGGTTCATATTTTCTGTCCCCACGGCGGCGAGCGGCGCTCCGAAGAAGAGTCCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTAcc                                            <  5:242980/133‑1 (MQ=255)
                                                                                              cacGGTCATGTGGAATATTTTCTGTCCCCTCTGCGGCGAGGTGCGCTCCGAAGAAGAGTACCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTAc                                        <  4:174261/136‑1 (MQ=255)
                                                                                                ccGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGCCCCAAACGCCTGCTCCGCCGAGGAGTGGGGTACCTACAt                                      >  8:208238/1‑136 (MQ=255)
                                                                                                    cATGTTGCATATTTTCTGCCCCCACTGCGGCGCGCGGCGCTCCGAACAAGAGTGCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTACATGtt                                   <  2:220277/135‑1 (MQ=255)
                                                                                                           cATATTTTCTGGCCCCCCTGCGGCGAGCTGCGCTCCGAAGAAGATTGCCACGCCTCTGGCCAGGCGCACTTTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGTGTACCTACATGTTCTACCGt                            <  8:284009/135‑1 (MQ=255)
                                                                                                            atatTTTCTGTCCCCACTGCGGCGCGCGGCGCGCCGAAGAAGAGTTCCACGCCTCTCGCCACGCGCACATTGCCCCCCCGCGGGCCCCGAACGCCGGCTCCGACGCGGAGAGGGGTACCGACATGTTCTACCGTg                           >  5:311532/1‑135 (MQ=255)
                                                                                                              atTTTCTGTCCCCACTGCGGCGCGCTGCCCGCCCAACAACACTTACACCCCGCAGGCCCCGCGCACAATGCCCGCCCCCAGGCCACAAACGCCAGCGCCCCCGAGGAGGGGGTGACCAACCATGTCTACCCTGACa                        >  4:86646/1‑136 (MQ=255)
                                                                                                              atTTTCTGTCCCCACTGCGGCGAGCAGCGCGCCGAAGAAGAGTTCCACCCCTCTGGCCAGGCGCCCATTTCCCCCCCCCGGGCCCCGAACGCCGGCTCCCACGAGGCGTTGGGTACCTACATTTTCTAACGTgg                          >  2:341299/1‑133 (MQ=255)
                                                                                                                    cTGACCCCACTGCGGAGAGCTGGGCTCCGAAGAAGAGTCCCAGGGCTGTGGCCAGGCGCACATTGGCCGCCCGCTGGACCCTAACGCCTGCTCCGTCGAGGAGTGGGGTACCTACATGTTCTACCGTGACAAccg                    <  3:155062/135‑2 (MQ=255)
                                                                                                                     tgcccccccGGCGGCGTGCGGCCCTCGGAACACGCGTACCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTACATGTTCTACCGTGACAAccc                    <  1:341299/131‑1 (MQ=255)
                                                                                                                     tgcccccACTGCGGCGAGCGGCGCGCCGAATAAGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTACATGTTCTACCGTGACAAccc                    <  1:110826/131‑1 (MQ=255)
                                                                                                                     tGTCCCCACTGCGGCGAGCTGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTACATGTTCTACCGTGACAAccc                    <  1:40575/134‑1 (MQ=255)
                                                                                                                      gTCCCCACTGCGGCGGGCTGCGCTCCGAAGAAGAGTTGCACGCCTCTGGCCAGGCGCCCATTGCCCGCCCGCTGGGCCCGAACGCCTGCTCCGGCGAGGGGTGGGGTAGCTACATGTTGTACCGTGGCACCCcgc                  >  2:91052/1‑135 (MQ=255)
                                                                                                                         ccgcgtGCGGCTCGCGGCGCTCCGACTTAGTGCGCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTACATGTTCTACCGTGACATCCCGCGCgg              <  7:337293/131‑1 (MQ=255)
                                                                                                                         ccccCTGCGGCGAGCTGCGCTCCGAAGAAGCGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTACATGTTCTACCGTGACAACCCGCGCgg              <  8:67275/136‑1 (MQ=255)
                                                                                                                          ccACTGCGGCGAGCGGCGCTCCGACGACGCGCCCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTACATGTTCTACCGTGACAACCCGCGCgg              <  5:259354/135‑1 (MQ=255)
                                                                                                                            aCTGCGGCGCGCGGCGCTCCGTAGAAGACTTCCCCGCCTCTGGCCAGGCGCACATTGCCCGCCCGCGGGCCCCTAACCCCTGCTCCGACGAGGAGTGGGGTTCCTACATGTTCTACCGTGACCACCCGCGCggga            >  1:50847/1‑133 (MQ=255)
                                                                                                                                 ggcgAGCTGCGCTCCGAAGAAGAGTGCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTACATGTTCTACCGTGACAACCCGCGCGGTATTCAcc      <  7:208238/136‑1 (MQ=255)
                                                                                                                                      gCTGCGCTCCGAAGAAGAGTTACACGCCTCTGGCCCGGCGCACATTTCCCGCCCGCTGGGCCCTAACGCCTGCTCCGCCGAGGAGTGGGGTACCGACATGTTCTACCGTGCCAACCCGCGCGGGATTCACCATGa  >  8:254453/1‑135 (MQ=255)
                                                                                                                                        |                                                                                                                                     
CATCCACTGGCCGCGCCGTTCTCCATGGACCGTTTCTACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCGCACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTACATGTTCTACCGTGACAACCCGCGCGGTATTCACCATGA  >  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: