Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 4,294,789 G→T 7.5% Q371K (CAG→AAG)  soxA_2 ← Monomeric sarcosine oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1434,294,7890GT 7.5% 161.4 / 4.2 62Q371K (CAG→AAG) soxA_2Monomeric sarcosine oxidase
Reads supporting (aligned to +/- strand):  ref base G (30/27);  new base T (2/3);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.36e-01

TAACTCGTGGCGCTTTTTGGTTGTATTTCCAAAACACCATACTAGCCAGCACCGATTCCTGTTTCAGGATAGGTAGACTGTGACCCGCGCACGCTAGTGGGAAAACCGCTCCCGGCTGAAGGCAGGAACCTGTAAGGCCAACGACTGGCCCGCTAAAGCGCGGGCCAGGCAACTACCCAATGCCGGCGCAAACTTGAAACCATGGCCGGAACAGGCCGAGACATAGAATACACTTGGGAGCTCCTTGCTTTCCCCTATGATGAA  >  ACB143/4294658‑4294921
                                                                                                                                   |                                                                                                                                    
tAACTCGTGGCGCTTTTTGGTTGTATTTCCAAAACACCATACTAGCCAGCACCGAATCCTGTTTCAGGATAGGTAGACTGTGACCCGCGCACGCTAGTGGGAAAACCGCTCCCGGCTGAAGGCAGGAACCTGTAAg                                                                                                                                  <  8:17825/136‑1 (MQ=255)
 aaCTCGTGGCGCTTTTTGGTTGTATTTCCAAAACACCATACTAGCCAGCACCGATTCCTGTTTCAGGATAGGTAGACTGTGACCCGCGCACGCTAGTGGGAAAACCGCTCCCGGCTGAAGGCAGGAACCTGTAAg                                                                                                                                  <  4:333150/135‑1 (MQ=255)
         gcgcTTTTTGGTTGTATTTCCAAAACACCATACTAGCCAGCACCGATTCCTGTTTCAGGATAGGTAGACTGTGACCCGCGCACGCTAGTGGGAAAACCGCTCCCGGCTGAAGGCAGGAACCTGTAAGGCCAACGa                                                                                                                          <  6:374901/135‑1 (MQ=255)
           gcTTTTTGGTTGTATTTCCAAAACACCATACTAGCCAGCACCGATTCCTGTTTCAGGATAGGTAGACTGTGACCCGCGCACGCTAGTGGGAAAACCGCTCCCGGCTGAAGGCAGGAACCTGTAAGGCCAACGACTg                                                                                                                       >  1:313654/1‑136 (MQ=255)
               tttGGTTGTATTTCCAAAACACCATACTAGCCAGCACCGATTCCTGTTTCAGGATAGGTAGACTGTGACCCGCGCACGCTAGTGGGAAAACCGCTCCCGGCTGAAGGCAGGAACCTGTAAGGCCAACGACTGGc                                                                                                                     <  5:266687/134‑1 (MQ=255)
                      gTATTTCCAAAACACCATACTAGCCAGCACCGATTCCTGTTTCAGGATAGGTAGACTGTGACCCGCGCACGCTAGTGGGAAAACCGCTCCCGGCTGAAGGCAGGAACCTGTAAGGCCAACGACTGGCCCGCTaaa                                                                                                             <  4:6330/135‑1 (MQ=255)
                        aTTTCCAAAAAACCTTACTAGCAAGCACCGATTCATTTTTCTGGATAGGTAGACTGGGCCCGGCGCACGCTTGTGGTAAAACCGCTCCCGGCTGAAGGCAGTAACCTGTAAGGCCAACGACTGGCCCGATAAAgc                                                                                                           <  5:324563/135‑1 (MQ=255)
                          ttCCAAAACACCATACTAGCCAGCACCGATTCCTGTTTCAGGATAGGTAGACTGTGACCCGCGCACGCTAGTGGGAAAACCGCTCCCGGCTGAAGGCAGGAACCTGTAAGGCCAACGACTGGCCCGCTAAAGcgcg                                                                                                        >  2:342201/1‑136 (MQ=255)
                                aaCACCATACTAGCCAGCACCGATTCGTGTTTCTGGATAGGTCGACTGTGACCCGCGCACGCTAGTGGGAAAACCGCCCCCGGCTGAAGGCAGCAACCTGTAAGGCCAACGACTGGCCCGCTAAAGCGCGGGCCAg                                                                                                  <  4:234488/136‑1 (MQ=255)
                                 acacCATACTAGCCCGCACCGATTCCTGTTTAAGGATAGGTAGACTGTGACCCGCGCACGCTAGTGGGAAAACCGCTCCCGGCTGAAGGCAGGAACCTGTAAGGCCAACGACTGGCCCGCTAAAGCGCGGGCCAgg                                                                                                 <  3:297293/136‑1 (MQ=255)
                                    ccATACTAGCCAGCACCGATTCCTGTTTCAGGATAGGTAGACTGTGACCCGCGCACGCTAGTGGGAAAACCGCTCCCGGCTGAAGGCAGGAACCTGTAAGGCCAACGACTGGCCCGCTAAAGCGCGGGCCAGGCaa                                                                                              >  4:397064/1‑136 (MQ=255)
                                               aGCACCGATTCCTGTTTCAGGATAGGTAGACTGTGACCCGCGCACGCTAGTGGGAAAACCGCTCCCGGCTGAAGGCAGGAACCTGTAAGGCCAACGACTGGCCCGCTAAAGCGCGGGCCAGGCAACTACCCAATGc                                                                                   >  4:358607/1‑136 (MQ=255)
                                               aGCACCGATTCCTGTTTCAGGATAGGTAGACTGTGACCCGCGCACGCTAGTGGGAAAACCGCTCCCGGCTGAAGGCAGGAACCTGTAAGGCCAACGACTGGCCCGCTAAAGCGCGGGCCAGGCAACTACCCAATGc                                                                                   >  4:191704/1‑136 (MQ=255)
                                               aGCACCGATTCCTGTTTCAGGATAGGTAGACTGTGACCCGCGCACGCTAGTGGGAAAACCGCTCCCGGCTGAAGGCAGGAACCTGTAAGGCCAACGACTGGCCCGCTAAAGCGCGGGCCAGGCAACTACCCAATGc                                                                                   >  8:192630/1‑136 (MQ=255)
                                                  aCCGATTCCTGTTTCAGGAAAGGTAGACTGTGATCCGCGCACGCTAGTGGGAAAACCGCTCCCGGCTGAAGGCAGTAACCTGTAAGGCCAACGACTGGCCCGCTAAAGCGCGGGCCAGGCAACTACCCAATGCCgg                                                                                <  2:261610/136‑1 (MQ=255)
                                                    cGATTCCTGTTTCAGGATAGGTAGACTGTGACCCGCGCACGCTAGTGGGAAAACCGCTCCCGGCTGAAGGCAGGAACCTGTAAGGCCAACGACTGGCCCGCTAAAGCGCGGGCCAGGCAACTACCCAATGCCGgcg                                                                              >  5:75934/1‑136 (MQ=255)
                                                    cGATTCCTGTTTCAGGATAGGTAGACTGTGACCCGCGCACGCTAGTGGGAAAACCGCTCCCGGCTGAAGGCAGGAACCTGTAAGGCCAACGACTGGCCCGCTAAAGCGCGGGCCAGGCAACTACCCAATGCCGgcg                                                                              >  1:179778/1‑136 (MQ=255)
                                                    cGATTCCTGTTTCAGGATAGGTAGACTGTGACCCGCGCACGCTAGTGGGAAAACCGCTCCCGGCTGAAGGCAGGAACCTGTAAGGCCAACGACTGGCCCGCTAAAGCGCGGGCCAGGCAACTACCCAATGCCGgcg                                                                              >  5:284699/1‑136 (MQ=255)
                                                    cGATTCCTGTTTCAGGATAGGTAGACTGTGACCCGCGCACGCTAGTGGGAAAACCGCTCCCGGCTGAAGGCAGGAACCTGTAAGGCCAACGACTGGCCCGCTAAAGCGCGGGCCAGGCAACTACCCAATGCCGgcg                                                                              >  2:6221/1‑136 (MQ=255)
                                                         ccTGTTTCAGGATAGGTAGACTGTGACCCGCGCACGCTAGTGGGAAAACCGCTCCCGGCTGAAGGCAGGAACCTGTAAGGCCAACGACTGGCCCGCTAAAGCGCGGGCCAGGCAACTACCCAATGCCGGCGCaaa                                                                          >  4:95119/1‑135 (MQ=255)
                                                         ccTGTTTCAGGATAGGTAGACTGTGACCCGCGCACGCTAGTGGGAAAACCGCTCCCGGCTGAAGGCAGGAACCTGTAAGGCCAACGACTGGCCCGCTAAAGCGCGGGCCAGGCAACTACCCAATGCCGGCGCAAAc                                                                         >  4:272143/1‑136 (MQ=255)
                                                         ccTGTTTCAGGATAGGTAGACTGTGACCCGCGCACGCTAGTGGGAAAACCGCTCCCGGCTGAAGGCAGGAACCTGTAAGGCCAACGACTGGCCCGCTAAAGCGCGGGCCAGGCAACTACCCAAAGCCGGCGCAAAc                                                                         >  6:115776/1‑136 (MQ=255)
                                                         ccTGTTTCAGGATAGGTAGACTGTGACCAGCGCACGCTAGTGGGAAAACCGCTCCCGGCTGAAGTCAGGAACCTGTAAGGCCAACGACTGGCCCGCTAAAGCGCGGGCCAGGCAACTACCCAATGCCGGCGCaaa                                                                          >  1:133867/1‑135 (MQ=255)
                                                                        gTAGACTGTGACCCGCGCACGCTAGTGGGAAAACCGCTCCCGGCTGAAGGCAGGAACCTGTAAGGCCAACGACTGGCCCGCTAAAGCGCGGGCCAGGCAACTACCCAATGCCGGCGCAAACTTGAAACCATGGCCg                                                          <  5:157602/136‑1 (MQ=255)
                                                                        gTAGACGGTGACCCGCGCACGCTAGTGGTAAAACCGCTCCCGGCTGAAGGCAGGAACCTGTAAGGCCAACGACTGGCCCGCTAAAGCGCGGGCCAGGCAACTACCCAATGCCGGCGCAAACTTGAAACCATGGCCg                                                          <  7:37702/136‑1 (MQ=255)
                                                                             ctgtcccccGCCCCCGCTAGTGGCAAACCCGCCCCCGGCTGAAGCCAGCAACTTGTACGCCCAACCACTGGCCCGCAAAAGCGCGGGCCAGGCAACTACCCAATGCCGGCGCAAACTTGAAACCATGGCCGGAACa                                                     <  2:181765/130‑1 (MQ=255)
                                                                                tgACCCGCGCACTCTAGTGGTAAAACCGCTCCCGGCTGAAGGCAGGAACATGTAAGGCCAACTACTGGCCCGCTAAAGCGCGGGCCAGGCAACTACCCAATGCCGGCGCAAACTTGAAACCATGGCCGGAACAgg                                                   <  7:126603/135‑1 (MQ=255)
                                                                                 gACCCGCGCACGCTAGTGGTAAAACCGCTCCCGGCTGAAGGCAGGAACCTGTAAGGCCAACGACTGGCCCGCTAAAGCGCGGGCCAGGCAACTACCCAATGCCGGCGCAAACTTGAAACCATGGCCGGAACAGGcc                                                 <  1:12027/136‑1 (MQ=255)
                                                                                  aCCCGCGCACGCTAGTGGGAAAACCGCTCCCGGCTGAAGGCAGGAACCTGTAAGGCCAACGACTGGCCCGCTAAAGCGCGGGCCAGGCAACTACCCAATGCCGGCGCAAACTTGAAACCATGGCCGGAACAGGcc                                                 >  7:301231/1‑135 (MQ=255)
                                                                                  aCCCGCGCACGCTAGTGGGAAAACCGCTCCCGGCTGAAGGCAGGAACCTGTAAGGCCAACGACTGGCCCGCTAAAGCGCGGGCCAGGCAACTACCCAATGCCGGCGCAAACTTGAAACCATGGCCGGAACAGGcc                                                 >  7:301237/1‑135 (MQ=255)
                                                                                  aCCCGCGCACGCTAGGGGGAAAACCGCTCCCGGCTGAAGGCAGGAACCTGTAAGGCCAACGACTGGCCCGCTAAAGCGCGGGCCAGGCAACTACCCAATGCAGGCGCAAACTTGAAACCATGGCCGGAACAGGcc                                                 <  1:341563/135‑1 (MQ=255)
                                                                                   cacgcgCACGCTAGTGGGAAAACCGCTCCCGGCTGAAGGCAGGACCCTGTAAGGCCAACGACTGGCCCGCAAAAGCGCGGGCCAGTCACATACCCAATGCCGGAGCAAACTAGAAACCATGGCCGGAACAGGcc                                                 >  4:149029/3‑134 (MQ=255)
                                                                                    ccGCGCACGCTAGTGGTAAAACCGCTCCCGGCTGAAGGCAGGAACCTGTAAGGCCAACGACTGGCCCGCTAAAGCGCGGGCCAGGCAACTACCCAATGCCGGCGCAAACTTGAAACCATGGCCGGAACAGGCCgag                                              <  1:99471/136‑1 (MQ=255)
                                                                                    ccGCGCACGCTAGGGGCAAAACCGCCCCCGGCTGAAGGCAGGAACCTGTACGGCCAACGACTGGCCCGCTAAAGCGCGGGCCAGGCAACTACCCAATGCCGGCGCAAACTTGAAACCATGGCCGGAACAGGCCgag                                              <  5:149930/136‑1 (MQ=255)
                                                                                     cgcgCACGCTAGTGGGAAAACCGCTCCCGGCTGAAGGCAGGAACCTGTAAGGCCAACGACTGGCCCGCTAAAGCGCGGGCCAGGCAACTACCCAATGCCGGCGCAAACTTGAAACCATGGCCGGAACAGGCCgaga                                             >  2:205390/1‑136 (MQ=255)
                                                                                     cgcgCACGCTAGTGGGAAAACCGCTCCCGGCTGAAGGCAGGAACCTGTAAGGCCAACGACTGGCCCGCTAAAGCGCGGGCCAGGCAACTACCCAATGCCGGCGCAAACTTGAAACCATGGCCGGAACAGGCCgaga                                             >  1:94344/1‑136 (MQ=255)
                                                                                        gcACGCTAGTGGGAAAACCGCTCCCGGCTGAAGGCAGGAACCTGTAAGGCCAACGACTGCCCCGCTAACGCGCGGGCCAGGCAACTACCCAATGCCGGCGCAAACTTGAAACCATGGCCGGAACAGGCCGAGACAt                                          <  3:262870/136‑1 (MQ=255)
                                                                                          ccGCTTGTGGTAAAACCGCTCCCGGCTGAAGGCAGGAACCTGTAAGGCCAACGACTGGCCCGCTAACGCGCGGGCCAGGCAACTACCCAATGCCGGCGCAAACTTGAAACCATGGCCGGAACAGGCCGAGACATAg                                        <  2:289583/135‑1 (MQ=255)
                                                                                          aCGCTAGTGGGAAAACCGCTCCCGGCTGAAGGCAGGAACCTGTAAGGCCAACGACTGGCCCGCTAAAGCGCGGGCCAGGCAACTCCCCAATGCCGGCGCAAACTTGAAACCATGGCCGGAACAGGCCGAGACa                                           >  2:81188/1‑133 (MQ=255)
                                                                                          aCGCTAGTGGGAAAACCGCTCCCGGCTGAAGGCAGGAACCTGTAAGGCCAACGACTGGCCCGCTAAAGCGCGGGCCAGGCAACTACCCAATGCCGGCGCAAACTTGAAACCATGGCCGGAACAGGCCGAGACa                                           >  1:123164/1‑133 (MQ=255)
                                                                                          aCGCTAGTGGGAAAACCGCTCCCGGCTGAAGGCAGGAACCTGTAAGGCCAACGACTGGCCCGCTAAAGCGCGGGCCAGGCAACTACCCAATGCCGGCGCAAACTTGAAACCATGGCCGGAACAGGCCGAGACa                                           >  2:105814/1‑133 (MQ=255)
                                                                                          aCGCTAGTGGGAAAACCGCTCCCGGCTGAAGGCAGGAACCTGTAAGGCCAACGACTGGCCCGCTAAAGCGCGGGCCAGGCAACTACCCAATGCCGGCGCAAACTTGAAACCATGGCCGGAACAGGCCGAGACa                                           >  3:191870/1‑133 (MQ=255)
                                                                                                    taaaaCCGCTCCCGGCTCACGGCAGTACCCTTTACGCCCACCGACTGCCCCGATAAAGCGCGGCCCAGGCAACTACCCAATGCCGGCGCAAACTTGAAACCATGGCCGGAACAGGCCGAGACATAGAATACACTTg                              <  6:213610/135‑1 (MQ=255)
                                                                                                     caacccGGTCCCGGCTTCAGGCAGGAACCTGTAAGGCGACGGACTGGCCCGCTCAAGCGCGGGCCAGGCAACTACCCAATGCCGGCGGAAACTTGAAACCATGGCCGGAACAGGCCGAGACATAGAATACACTTgg                             <  6:70791/132‑1 (MQ=255)
                                                                                                      aaaCCGCTCCCGGCTGAAGGCAGGAACCTGTAAGGCCAACGACTGGCCCGCTAAAGCGCGGGCCAGGCAACTACCCAATGCCGGCGCAAACTTGAAACCATGGCCGGAACAGGCCGAGACATAGAATACACTTgg                             <  4:327300/135‑1 (MQ=255)
                                                                                                         ccGCTCCCGGCTGACGGCAGGAACCTGTAAGGCCAACGACTGGCCCGCTAAAGCGCGGGCCAGGCAACTACCCAATGCCGGCGCAAACTTGAAACCATGGCCGGAACAGGCCGAGACATAGAATACACTTGGGAGc                         <  6:326753/136‑1 (MQ=255)
                                                                                                         ccGCTCCCGGCTGAAGGCAGGAACCTTTAAGGCCAACGACTGGCCCGCTAAAGCGCggg                                                                                                      <  2:90126/59‑1 (MQ=255)
                                                                                                         ccGCTCCCGGCTGAAGGCAGGAACCTTTAAGGCCAACGACTGGCCCGCTAAAGCGCggg                                                                                                      >  1:90126/1‑59 (MQ=255)
                                                                                                         ccGCTCCCGGCTGAAGGCAGGAACCTTTAAGGCCAACGACTGGCCCGCTAAAGCGCggg                                                                                                      >  1:96141/1‑59 (MQ=255)
                                                                                                         ccGCTCCCGGCTGAAGGCAGGAACCTTTAAGGCCAACGACTGGCCCGCTAAAGCGCggg                                                                                                      <  2:96141/59‑1 (MQ=255)
                                                                                                         ccGCTCCCGGCTGAAGGCAGGAACCTGTAAGGCCAACGACTGGCCCGCTAAAGCGCGGGCCAGGCAACTACCCAATGCCGGCGCAAACTTGAAACCATGGCCGGAACAGGCCGAGACATAGAATACACTTGGGAGc                         >  1:298928/1‑136 (MQ=255)
                                                                                                         ccGCCCCCGGCTGAAGGCAGGAACCTGTACGGCCAACGACTGGCCCGATAACGCGCGGGCCAGGCAACTACCCAATGCCGGCGCAAACTTGAAACCATGGCCGGAACAGGCCGAGACATAGAATACACTTGGGAGc                         <  1:18276/136‑1 (MQ=255)
                                                                                                           gCTCCCGGCTGAAGGCAGGAACCTGTAAGGCCAACGACTGGCCCGCTAAAGCGCGGGCCAGGCAACTACCCAATGCCGGCGCAAACTTGAAACCATGGCCGGAACAGGCCGAGACATAGAATACACTTGGGAGCTc                       >  5:63439/1‑136 (MQ=255)
                                                                                                             tCCCGGCTGAAGGCAGGAACCTGTAAGGCCAACGACTGGCCCGCTAAAGCGCGGGCCAGGCAACTACCCAATGCCGGCACAAACTTGAAACCATGGCCGGAACAGGCCGAGACATAGAATACACTTGGGAGCTCCt                     >  7:62464/1‑136 (MQ=255)
                                                                                                              cccGGCTGAAGGCAGGAACCTGTAAGGCCAACGACTGGCCCGCTAAAGCGCGGGCCAGGCAACTACCCAATGCCGGCGCAAACTTGAAACCATGGCCGGAACAGGCCGAGACATAGAATACACTTGGGAGCTCCtt                    <  2:94344/136‑1 (MQ=255)
                                                                                                              cccGGCTGAAGGCAGGAACCTGTAAGGCCAACGACTGGCCCGCTAAAGCGCGGGCCAGGCAACTACCCAATGCCGGCGCAAACTTGAAACCATGGCCGGAACAGGCCGAGACATAGAATACACTTGGGAGCTCCtt                    <  4:47668/136‑1 (MQ=255)
                                                                                                                   cTGAAGGAAGCAACCTGAACGGCCAAGGACTGGCCCGCTAAAGCGCGGGCCTGTCAACTACCCAATGCCGGCGCAAACTTGAAACCATGGCCGGAACAGGCCGAGACATAGAATACACTTGGGAGCTCCTTGCttt               <  5:158915/136‑1 (MQ=255)
                                                                                                                         gCAGGAACCTGTAAGGCCAACGACTGGCCCGCTAAAGCGCGGGCCAGGCAACTACCCAATGCCGGCGCAAACTTGAAACCATGGCCGGAACAGGCCGAGACATAGAATACACTTGGGAGCTCCTTGCTTTCCCCTa         >  3:290660/1‑136 (MQ=255)
                                                                                                                                ccTGTAAGGCCAACGACTGGCCCGCTAAAGCGCGGGCCAGGCAACTACCCAATGCCGGCGCAAACTTGAAACCATGGCCGGAACAGGCCGAGACAt                                          <  7:373443/96‑1 (MQ=255)
                                                                                                                                ccTGTAAGGCCAACGACTGGCCCGCTAAAGCGCGGGCCAGGCAACTACCCAATGCCGGCGCAAACTTGAAACCATGGCCGGAACAGGCCGAGACAt                                          >  8:373443/1‑96 (MQ=255)
                                                                                                                                ccTGTAAGGCCAACGACTGGCCCGCTAAAGCGCGGGCCAGGCAACTACCCAATGCCGGCGCAAACTTGAAACCATGGCCGGAACAGGCCGAGACATAGAATACACTTGGGAGCTCCTTGCTTTCCCCTATGATGaa  >  3:314847/1‑136 (MQ=255)
                                                                                                                                ccTGTAAGGCCAACGACTGGCCCGCTAAAGCGCGGGCCAGGCAACTACCCAATGCCGGCGCAAACTTGAAACCATGGCCGGAACAGGCCGAGACATAGAATACACTTGGGAGCTCCTTGCTATCCCCTATGATGaa  >  1:64301/1‑136 (MQ=255)
                                                                                                                                   |                                                                                                                                    
TAACTCGTGGCGCTTTTTGGTTGTATTTCCAAAACACCATACTAGCCAGCACCGATTCCTGTTTCAGGATAGGTAGACTGTGACCCGCGCACGCTAGTGGGAAAACCGCTCCCGGCTGAAGGCAGGAACCTGTAAGGCCAACGACTGGCCCGCTAAAGCGCGGGCCAGGCAACTACCCAATGCCGGCGCAAACTTGAAACCATGGCCGGAACAGGCCGAGACATAGAATACACTTGGGAGCTCCTTGCTTTCCCCTATGATGAA  >  ACB143/4294658‑4294921

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