Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 382,484 T→G 28.2% S14A (TCC→GCC)  soxD → Sarcosine oxidase subunit delta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT2440382,4840TG28.2% 44.5 / 2.4 56S14A (TCC→GCC) soxDSarcosine oxidase subunit delta
Reads supporting (aligned to +/- strand):  ref base T (10/27);  new base G (5/12);  total (16/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGCCGCGCCGTTCTCCATGGACCGTTTCTACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCG‑CACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTACATGTTCTACCGTGACAACCCGCGCGGTATTCACCATGAACTGTG  >  KT2440/382350‑382617
                                                                                                                                      |                                                                                                                                      
cTGGCCGCGCCGTTCTCCATGGACCGTTTCTACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGAGCTGCGCt                                                                                                                                        <  6:143440/135‑1 (MQ=255)
    ccGCGCCGTTCTCCATGGACCGTTTCTACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCTAGCTGCGCTCCgaa                                                                                                                                   <  8:308088/136‑1 (MQ=255)
     cgcgCCGTTCTCCATGGACCGTTTCTACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCgaa                                                                                                                                   <  2:9251/135‑1 (MQ=255)
      gcgcCGTTCTCCATGGACCGTTTCTACAACGGCGCACTGATCGACGCACACGGCGCCGCCGCCGTCGCCCACTACCCGGAGCCACCGGCATGTTGCATAGTTTCTGTCCCCACTGCGGCGGGCTGCGGGCCgaa                                                                                                                                   >  3:122924/1‑134 (MQ=255)
      gcgcCGTTCTCCATGGACCGTTTCTACAACGGCGCACTGATCGACGAACCCGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGCGCTGCGCGCCgca                                                                                                                                   >  4:349914/1‑132 (MQ=255)
             tctcCTTGGCCGGTTTATACACGGGCGCACTGGTCCCCCACCCCGGCGCCGCCGCCGTCGCCCAAAAACCGGACACACGGTATTGTTGCATATTTTCTGCCCCCACTGCGGCGAGCTGCGCTCCGAAGAAGAGtt                                                                                                                           <  4:28843/135‑1 (MQ=255)
             tctcCAGGGACGGTTGCTCCACCGGGGAAGTGAGCCACCACCACGGCGCCGCCGCCGTCGCCCACTAACGGGAGACACCGTCATGTTGAATTTTTTCTGCCCCCAGTGCGTCGAGCTGCGCTCCGAATAAGAttt                                                                                                                           <  2:322520/135‑1 (MQ=255)
                  aTGCACCGTTTCTACACCGGCGCACTGAGCGACCAACCCGGCGCCGCCGCCGGCGCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCGAAGAAGAGTTCCACg                                                                                                                      <  4:164230/135‑1 (MQ=255)
                       ccGTTTCTACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGAGCTGCGCGCCGAAGAAGAGTTCCACGCctct                                                                                                                 <  5:234257/135‑1 (MQ=255)
                         gTTTCTACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGGGCTGCGCTCCGAAGAAGAGTTCCCCCCCGCTgg                                                                                                               >  6:323609/1‑135 (MQ=255)
                         gTTTCTACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGCGCGGCGCTCCGCAGAAGAGTTCCCCCCCGCTGGc                                                                                                              >  7:15283/1‑136 (MQ=255)
                         gTTTCTACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCgaa                                                                                                                                   >  7:70821/1‑115 (MQ=255)
                         gTTTCTACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCgaa                                                                                                                                   <  8:70821/115‑1 (MQ=255)
                            tCTACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGGGCTGCGCTCCCAACAATAGTTACACCCCGCTGGccc                                                                                                            >  2:189253/1‑134 (MQ=255)
                            tCTACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTAGCATATTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCa                                                                                                            >  4:60932/1‑135 (MQ=255)
                                        cACTGATCGACGAACACGGAGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACCGCGGCGCGCTGCGCGCCGAAGAAGAGTTACACGCCACTGGCCCGGCG‑CCCAATgc                                                                                               >  4:379484/1‑135 (MQ=255)
                                              tCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCGAAGAATAGTTCCACGCCTCTGGCCAGGCG‑CACATTGcccgcccg                                                                                        <  7:81685/136‑1 (MQ=255)
                                              tCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTAATGTTGCATATTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCGAAGAAGAGTGCCACGCCTCTGGCCAGGCG‑CACATTgcccgccc                                                                                         <  5:284888/135‑1 (MQ=255)
                                              tCGACGAACACGGCGCCGCCGCCGGCGCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGCGCTGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCG‑CACATTGcccgcccg                                                                                        <  8:435962/136‑1 (MQ=255)
                                                gacgaACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGAGCGGCGCGCCGAAGCACAGTTGCCCGCCGCTGGCCAGGCG‑CACATTgcccgccccc                                                                                       >  1:367116/1‑133 (MQ=255)
                                                        cGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGAGCGGCGCTCCGCAGAACAGTTCCACGCCGCTGGCCACGCG‑GACATTGCCCGCCCGCGGGcccc                                                                                >  6:234388/1‑134 (MQ=255)
                                                            gccgccgccGTCGCCCACTAACCGGAGACCCCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGAGGGGCGCGCCCAAGAAGAGTTGCACGCCACGGGCCCGGGGCCCCATTG‑CCGCCCGCTGGACCCTAACg                                                                           >  6:189567/1‑135 (MQ=255)
                                                             ccgccgcggTCGCCCACTACCCGGCGACACGGTCATGTTCCATATTTTCTGACCCCTGTGGGGCGAGCTGCGCGCCTACGAAGTGTTCCACGCCTCTGGCCAGGCG‑CACATTGCCCGCCCGCTGGACCCTAACGc                                                                          <  2:56580/135‑1 (MQ=255)
                                                              cgccgcctGCGCCCACTAACCGGAGACACCGTCATGTTGCATTTTTGCTGCCCCCACTGCGGCGAGCGGCGCTCCGAACACGCGTGCCACGCCTCTGGCCAGGCG‑CACATTGCCCGCCCGCTGGACCCTAACGCCt                                                                        <  8:434991/136‑1 (MQ=255)
                                                                   ccgTCGCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCG‑CACATTGCCCGCCCGCTGGACCCTAACGCCTGCTcc                                                                   <  8:95157/136‑1 (MQ=255)
                                                                    cggcGCCCACTCCCCGGACCCACGGTATTGTTTATTATTTTCTGCCCCCCCGCCGGCGCGCGGCGCTCCCAAGAGGCTTCCCACGCCTCTGGCCAGGCG‑CACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCcg                                                                  <  7:404980/133‑1 (MQ=255)
                                                                       cGCCCACTAACCGGACACCCCGTCATGTTGAATATTTTCTGCCCCCCCTGCGGCGAGCGGCGCGCCGAAGACGCGTCCCACGCCTCTGGCCAGGCG‑CACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCcgac                                                                <  4:376835/135‑1 (MQ=255)
                                                                          ccACTACCCGGAGACCCCGTCATGTTCAATATTTTCTGCCCCCCCCGCGGCGAGCTGCGCTCCGACGAAGAGTTCCACGCCTCTGGCCAGGCG‑CACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACgag                                                             <  3:379484/135‑1 (MQ=255)
                                                                          ccACTACCCGGACACACCGTCATGTTGCATATTTTCGGCCCCCCCGGCGGCGAGCTGCGCTCCGACTAAGAGTCCCACGCCTCTGGCCAGGCG‑CACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACgag                                                             <  2:367116/135‑1 (MQ=255)
                                                                          ccACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCGAAGAAGAGTCCCACGCCTCTGGCCAGGCG‑CACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACgag                                                             <  3:221480/135‑1 (MQ=255)
                                                                           cACTAACCGGAGACAGCGTCATGTTTCATATTTTCTGTCCCCACTGCGGCGAGCGGCGCTCCTAAGAAGAGTCCCACGCCTCTGGCCAGGCG‑CACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACgagg                                                            <  3:58232/135‑1 (MQ=255)
                                                                            ccTAACCGGAGACACCGTCATGTTGAATATTTTCTGGCCCCACGGCGGCGAGCTGCGCGCCGAAGAAGAGGTCCACGCCTCTGGCCAGGCG‑CACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACgagga                                                           <  4:355164/134‑1 (MQ=255)
                                                                               gccgggAGCCACGGTCTTGTTGCATATTTTCTCCCCCCCCGGCGGCGAGCGGCGCCCCGAACATGAGTCCCACGCCTCTGGCCAGGCG‑CACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTg                                                        <  1:284903/131‑1 (MQ=255)
                                                                               aaCCGCAGACACCGTCATGTTCCATATTTTCTGCCCCCACTGCGGCGAGCGGCGCTCCGAAGAAGAGTGCCACGCCTCTGGCCAGGCG‑CACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTg                                                        <  3:183178/135‑1 (MQ=255)
                                                                                      caccccGTCTTGTTGCATATTTGCTCCCCCCCCGGCGGCGCGCGGCGCGCCTAAGACGAGTCCCACGCCTCTGGCCAGGCG‑CACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTAcc                                                 <  5:223361/131‑1 (MQ=255)
                                                                                       acacCGTCTGTTTTCATATTTTGTGTCCGCAGTGGGTCTCTGTCCGATCCTAATAAGAGTGCCACGCCTCTGGCCAGGCG‑CACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTAcc                                                 <  6:357604/134‑1 (MQ=255)
                                                                                         acCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCG‑CACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTa                                               >  3:85617/1‑134 (MQ=255)
                                                                                          ccGTCATGTTGCATATTTTCGGCCCCCACTGCGGCGAGCTGCGCGCCGAAGACGAGGTCCACGCCTCTGGCCAGGCG‑CACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTACAt                                            <  8:184784/136‑1 (MQ=255)
                                                                                              cATGTTGCATATTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCGAAGAAGAGTTCCACGCCAATGGGCAGGCG‑CCCATTGCCCGCCCGCTGGTCCCAAACGCCTGCGCCGACGAGGAGTGGGGTACCTACATggt                                         >  5:83684/1‑133 (MQ=255)
                                                                                              cATGTTGCATATTTTCTGTCCCCACTGCGGCGAGCTGCGCACCGAACAACAGTTACACGCCACTGGCCACGCG‑CACATTGCCCGCCCGCTGGCCCCGAAACCCTGCTCCCACCAGGCGGGGGGGTCCGACATGtt                                         >  1:250249/1‑135 (MQ=255)
                                                                                              cATGTTGCATATTTTCTGTCCCCACTGCGGCGAGCGGCGGTCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCG‑CACCTTTCCCGCCCGCTGGGCCCTAACGCCTGCTACGACCAGGAGTGGGGTACCTACATGtc                                         >  2:284903/1‑134 (MQ=255)
                                                                                                        atTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCCGGCG‑CCCATTGCCCGCCCGCTGGCCCCTAACGCCTGCTCCGCCGCGGAGTGGGGTACCTACATGTTCTACCGTgcc                               >  3:392507/1‑133 (MQ=255)
                                                                                                              tgcccccccTGGGGGGCGCTGCGGGCCGAAGACGGGTCCCCCGCGTCTGGCCAGGCG‑CACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTACATGTTCTACCGTGACAAccc                          <  5:67757/131‑1 (MQ=255)
                                                                                                              tgcccccACTGCGGCGAGCCGCGCTCCGAAGAACAGTCCCACGCCTCTGGCCAGGCG‑CACATTGACCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTACATGTTCTACCGTGACAAccc                          <  3:82045/131‑1 (MQ=255)
                                                                                                              tgaccccccTGCGGCGAGCGGCGCGCCGAAGCCGAGTCCCACGCCTCTGGCCAGGCG‑CACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTACATGTTCTACCGTGACAAccc                          <  3:45818/131‑1 (MQ=255)
                                                                                                              tGTCCCCACTGCGGCGCGCTGCGCTCCGTAGAAGCGTTCCACGCCTCTGGCCAGGCG‑CACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTACATGTTCTACCGTGACAAccc                          <  2:50334/134‑1 (MQ=255)
                                                                                                              tGTCCCCACTGCGGCGAGCTGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCG‑CACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTACATGTTCTACCGTGACAAccc                          <  4:85617/134‑1 (MQ=255)
                                                                                                                 ccccACTGCGGCGTGGTGCGCGCGGATCAACAGGCCCACGCCTCTGGCCAGGCG‑CACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTACATGTTCTACCGTGACAACCcgc                        <  3:60932/133‑1 (MQ=255)
                                                                                                                  ccccCTGCGTCGAGCTGCTCTCCGAATAAGAGTTCCACGCCTCTGGCGAGGCG‑CACATTGCGCGCCCGCTGGACCCTATCGGCTGCTCCGACGAGGAGTGGGGTACGTACTTGTTCTACCGTGAGATCCCGAGCgg                    <  8:245108/136‑1 (MQ=255)
                                                                                                                    cccgggggcgCGCTGCCCGCCTAACATGCGTCCCACGCCTCTGGCCAGGCG‑CACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTACATGTTCTACCGTGACAACCCGCGCGGt                   <  5:355824/129‑1 (MQ=255)
                                                                                                                    cACTGGGGCGCGGTGCGCTCCTAATAAGAGTCCCACGCCTCTGGCCAGGCG‑CACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTACATGTTCTACCGTGACAACCCGCGCGGt                   <  5:170764/135‑1 (MQ=255)
                                                                                                                          gggtcgCTGCGCTCGTACCAAGAGTTCGACGCCTCTGGCCAGGCG‑CAGATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGAGGAGTTGTGTGGTACCTACATGTTCTAGCGTGACAACGCGCGCGGTATTCAcc            <  7:377588/131‑1 (MQ=255)
                                                                                                                          ggcgAGCGGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCG‑CACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTACATGTTCTACCGTGACAACCCGCGCGGTATTCAcc            <  7:284038/136‑1 (MQ=255)
                                                                                                                           gcgcgcGGCGCGCCGAACAAGTGTCCCACGCCTCTGGCCAGGCG‑CACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTACATGTTCTACCGTGACAACCCGCGCGGTATTCAcc            <  1:208882/135‑1 (MQ=255)
                                                                                                                             gcgCTGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCG‑CACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTACATGTTCTACCGTGACAACCCGCGCGGTATTCACCAt          <  4:166706/133‑1 (MQ=255)
                                                                                                                             gAGCGGCGCGCCTAAGAAGAGTTCCACGCCTCTGGCCAGGCG‑CACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTACATGTTCTACCGTGACAACCCGCGCGGTATTCACCAt          <  2:250249/135‑1 (MQ=255)
                                                                                                                              aGCTGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCG‑CACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTACATGTTCTACCGTGACAACCCGCGCGGTATTCACCATg         <  4:417936/135‑1 (MQ=255)
                                                                                                                                      tCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCG‑CACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTACATGTTCTACCGTGACAACCCGCGCGGTATTCACCATGAACtgtg  >  3:99797/1‑134 (MQ=255)
                                                                                                                                      |                                                                                                                                      
CTGGCCGCGCCGTTCTCCATGGACCGTTTCTACAACGGCGCACTGATCGACGAACACGGCGCCGCCGCCGTCGCCCACTAACCGGAGACACCGTCATGTTGCATATTTTCTGTCCCCACTGCGGCGAGCTGCGCTCCGAAGAAGAGTTCCACGCCTCTGGCCAGGCG‑CACATTGCCCGCCCGCTGGACCCTAACGCCTGCTCCGACGAGGAGTGGGGTACCTACATGTTCTACCGTGACAACCCGCGCGGTATTCACCATGAACTGTG  >  KT2440/382350‑382617

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