Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 388,722 T→A 9.1% L198Q (CTG→CAG)  soxB → sarcosine oxidase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947388,7220TA 9.1% 155.3 / 6.1 55L198Q (CTG→CAG) soxBsarcosine oxidase subunit beta
Reads supporting (aligned to +/- strand):  ref base T (25/25);  new base A (3/2);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.69e-01

CAGGTTGCGGCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGCACATGGCCAAGC  >  NC_002947/388583‑388868
                                                                                                                                           |                                                                                                                                                  
cAGGTTGCGGCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGa                                                                                                                                           <  3:344118/149‑1 (MQ=255)
cAGGTTGCGGCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGa                                                                                                                                           <  7:171891/149‑1 (MQ=255)
  ggTTGCGGCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCAGGGCGTCggca                                                                                                                                         <  6:46997/149‑4 (MQ=255)
   gTTGCGGCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCt                                                                                                                                        <  2:187857/149‑1 (MQ=255)
            gAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGc                                                                                                                                >  7:242272/1‑148 (MQ=255)
                tCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAAcc                                                                                                                           <  8:242272/149‑1 (MQ=255)
                     tACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGt                                                                                                                      >  8:79028/1‑149 (MQ=255)
                        cTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGAt                                                                                                                   <  7:321012/149‑1 (MQ=255)
                           gACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCgg                                                                                                                >  3:133344/1‑149 (MQ=255)
                                cTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTcc                                                                                                           >  7:66673/1‑149 (MQ=255)
                                cTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTcc                                                                                                           >  3:184955/1‑149 (MQ=255)
                                 tCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCg                                                                                                          >  2:25347/1‑149 (MQ=255)
                                  cGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGTCGACGCCCAGGGCGTCggcagcg                                                                                                                                      <  5:261498/120‑7 (MQ=255)
                                  cGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGTCGACGCCCAGGGCGTCggcagcg                                                                                                                                      >  6:261498/1‑114 (MQ=255)
                                         cacCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGaa                                                                                                   >  7:189992/1‑148 (MQ=255)
                                             cGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACgg                                                                                              >  6:221891/1‑149 (MQ=255)
                                                    cGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTc                                                                                                            <  6:310626/128‑1 (MQ=255)
                                                    cGATCCTCGGCGCCACCGTACAGCGCCGTGGTAGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGaa                                                                                                                        >  1:345870/1‑116 (MQ=255)
                                                    cGATCCTCGGCGCCACCGTACAGCGCCGTGGTAGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGaa                                                                                                                        <  2:345870/116‑1 (MQ=255)
                                                        tcTCGGCGCCACCGTACCGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCg                                                                                   <  3:234361/148‑1 (MQ=255)
                                                          tCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGc                                                                                                         <  1:315052/125‑1 (MQ=255)
                                                          tCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGc                                                                                                         >  2:315052/1‑125 (MQ=255)
                                                          tCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGgt                                                                                 >  4:211109/1‑149 (MQ=255)
                                                          tCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGgt                                                                                 >  5:48197/1‑149 (MQ=255)
                                                          tCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCAGATCCAGCAAACCGAAGTCAACGGCTACCGCAAGGAAAACGGAGCGGACATCg                                                                                   >  4:234361/1‑147 (MQ=255)
                                                                  aCCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAAc                                                                         <  6:264835/149‑1 (MQ=255)
                                                                            gcCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGc                                                               <  3:211109/149‑1 (MQ=255)
                                                                            gcCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGc                                                               >  1:60358/1‑149 (MQ=255)
                                                                               gtggtgGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTc                                                            <  8:189992/149‑1 (MQ=255)
                                                                                tggtggCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCa                                                           <  7:79028/149‑1 (MQ=255)
                                                                                  gtggCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATc                                                         <  1:25347/149‑1 (MQ=255)
                                                                                     gCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGgc                                                      <  6:48197/149‑1 (MQ=255)
                                                                                       gTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGgcgc                                                                                           >  8:176543/1‑110 (MQ=255)
                                                                                       gTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGgcgc                                                                                           <  7:176543/110‑1 (MQ=255)
                                                                                                  cGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCg                                         <  2:60358/149‑1 (MQ=255)
                                                                                                        cGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTgg                                   >  3:331928/1‑149 (MQ=255)
                                                                                                               tGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACaaa                                                                  >  7:85946/1‑111 (MQ=255)
                                                                                                               tGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACaaa                                                                  <  8:85946/111‑1 (MQ=255)
                                                                                                                   gCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGt                        <  5:221891/149‑1 (MQ=255)
                                                                                                                     ttCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAAc                                                                         <  7:39323/98‑1 (MQ=255)
                                                                                                                     ttCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAAc                                                                         >  8:39323/1‑98 (MQ=255)
                                                                                                                       cGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAAc                     <  6:5937/148‑1 (MQ=255)
                                                                                                                        gCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTc                   >  7:330018/1‑149 (MQ=255)
                                                                                                                        gCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTc                   >  8:231925/1‑149 (MQ=255)
                                                                                                                           ccaccTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCg                 <  1:305684/144‑1 (MQ=255)
                                                                                                                            gcgcTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCggg               <  4:331928/149‑1 (MQ=255)
                                                                                                                            gcgcTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCggg               <  7:231925/149‑1 (MQ=255)
                                                                                                                             cgcTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTctg                 >  2:305684/1‑144 (MQ=255)
                                                                                                                                 gtcGACGCCCAGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGCACAt          >  8:79553/3‑149 (MQ=255)
                                                                                                                                  ccGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGCACATg         >  7:59028/1‑149 (MQ=255)
                                                                                                                                    gACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGCACATGGc       >  8:209611/1‑149 (MQ=255)
                                                                                                                                    gACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGCACATGGc       >  3:143625/1‑149 (MQ=255)
                                                                                                                                     aCGCCCAGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGCACATGGc       >  7:140094/1‑148 (MQ=255)
                                                                                                                                       gCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGCACATGGCCaa    <  4:184955/149‑1 (MQ=255)
                                                                                                                                         ccTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGCACATGGCCAAGc  >  6:89594/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                  
CAGGTTGCGGCCGAAATCCCGTACCTGGACTGCTCGAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTCCGGGCACATGGCCAAGC  >  NC_002947/388583‑388868

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