Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,774,913 C→G 9.2% I260M (ATC→ATG)  betB → NAD‑dependent betaine aldehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,774,9130CG 9.2% 148.9 / 3.7 52I260M (ATC→ATG) betBNAD‑dependent betaine aldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (24/22);  new base G (2/3);  total (27/25)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCACCTGGCTGACCGAGCACCCGCGCATCGAAAAAGTGTCC‑TTCACCGGCGGCACCACCACCGGCAAGAAAGTCATGGCCAGCGCCT‑CCAGCTCTTCGCTGAAGGAAGTGACCATGGAGCTGGGCGGCAAGTCGCCACTGATCATCTGTGCCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCATGATGGCCAACTTCTACAGCTCGGGCCAGGTGTGCACCAACGGCACCCGCGTGTTCATCCCGGCTGAAATGAAGGCGGCCTTCGAAGCCAAGATCGCCGAG  >  NC_002947/5774770‑5775057
                                                                                                                                                 |                                                                                                                                                 
ggCACCTGGCTGACCGAGCACCCGCGCATCGAAAAAGTGTCC‑TTCACCTGCGGCACCACCACCGGCAAGAAAGTCATGGCCAGCGCCT‑CCAGCTCTTCGCTGAAGGAAGTGACCAAGGAGCTGGGCGGCAAGTGGGCACTGATAATCtt                                                                                                                                              >  8:124908/1‑148 (MQ=255)
ggCACCTGGCTGACCGAGCACCCGCGCATCGAAAAAGTGTCC‑TTCACCGGCGGCACCACCACCGGCAAGAAAGTCATGGCCAGCGCCT‑CCAGCTCTTCGCTGAAGGAAGTGACCATGGAGCTGGGCGGCAAGTCGCCACTGATCATCtg                                                                                                                                              >  7:170187/1‑149 (MQ=255)
ggCACCTGGCTGACCGAGCACCCGCGCATCGAAAAAGTGTCC‑TTCACCGGCGGCACCACCACCGGCAAGAAAGTCATGGCCAGCGCCT‑CCAGCTCTTCGCTGAAGGAAGTGACCATGGAGCTGGGCGGCAAGTCGCCACTGATCATCt                                                                                                                                               <  6:325790/148‑1 (MQ=255)
 gCACCTGGCTGACCGAGCACCCGCGCATCGAAAAAGTGTCC‑TTCACCGGCGGCACCACCACCGGCAAGAAAGTCATGGCCAGCGCCT‑CCAGCTCTTCGCTGAAGGAAGTGACCATGGAGCTGGGCGGCAAGTCGCCACTGATCATCtgt                                                                                                                                             <  5:31880/149‑1 (MQ=255)
       ggCTGACCGAGCACCCGCGCATCGAAAAAGTGTCC‑TTCACCGGCGGCACCACCACCGGCAAGAAAGTCATGGCCAGCGCCT‑CCAGCTCTTCGCTGAAGGAAGTGACCATGGAGCTGGGCGGCAAGTCGCCACTGatca                                                                                                                                                  <  3:299004/138‑1 (MQ=255)
       ggCTGACCGAGCACCCGCGCATCGAAAAAGTGTCC‑TTCACCGGCGGCACCACCACCGGCAAGAAAGTCATGGCCAGCGCCT‑CCAGCTCTTCGCTGAAGGAAGTGACCATGGAGCTGGGCGGCAAGTCGCCACTGatca                                                                                                                                                  >  4:299004/1‑138 (MQ=255)
        gCTGACCGAGCACCCGCGCATCGAAAAAGTGTCC‑TTCACCGGCGGCACCACCACCGGCAAGAAAGTCATGGCCAGCGCCT‑CCAGCTCTTCGCTGAAGGAAGTGACCATGGAGCTGGGCGGCAAGTCGCCACTGATCATCTGTgccgacg                                                                                                                                      >  3:163390/1‑149 (MQ=255)
         cTGACCGAGCACCCGCGCATCGAAAAAGTGTCC‑TTCACCGGCGGCACCACCACCGGCAAGAAAGTCATGGCCAGCGCCT‑CCAGCTCTTCGCTGAAGGAAGTGACCATGTAGCTGGGCGGCAAGTCGCCACTGATCATCTGTgccgacgc                                                                                                                                     <  8:170187/149‑1 (MQ=255)
          tGACCGAGCACCCGCGCATCGAAAAAGTGTCCTTTCACCGGCGGCACCACCACCGGCAAGAAAGTCATGGCCAGCGCCT‑CCAGCTCTTCGCTGAAGGAAGTGACCATGGAGCTGGGCGGCAAGTCGCCACTGATCATCTGTgccgatgc                                                                                                                                     >  1:47184/1‑149 (MQ=255)
                     ccGCGCATCGAAAAAGTGTCC‑TTCACCGGCGGCACCACCACCGGCAAGAAAGTCATGGCCAGCGCCT‑CCAGCTCTTCGCTGAAGGAAGTGACCATGGAGCTGGGCGGCAAGTCGCCACTGATCATCTGTGCCGATGCCGACCTGG‑ACaa                                                                                                                        <  2:47184/149‑1 (MQ=255)
                     ccGCGCATCGAAAAAGTGTCC‑TTCACCGGCGGCACCACCACCGGCAAGAAAGTCATGGCCAGCGCCT‑CCAGCTCTTCGCTGAAGGAAGTGACCATGGAGCTGGGCGGCAAGTCGCCACTGATCATCTGTGCCGACGCCGACCTGG‑Ac                                                                                                                          >  6:275801/1‑147 (MQ=255)
                       gcgcATCGAAAAAGTGTCC‑TTCACCGGCGGCACCACCACCGGCAAGAAAGTCATGGCCAGCGCCT‑CCAGCTCTTCGCTGAAGGAAGTGACCATGGAGCTGGGCGGCAAGTCGCCACTGATCATCTGTGCCGACGCCGACCTGG‑ACAAgg                                                                                                                      >  2:114428/1‑149 (MQ=255)
                               aaaaaGTGTCC‑TTCACCGGCGGCACCACCACCCGCAAGAAAGGCATGGCCAGCGCCT‑CCAGCTCTTCGCTGAAGGAAGGGACCATGGAGCTGGGCGGCAAGTCGCCACTGATCATCTGTGCCGACGCCGACCTTG‑ACAAGGCCGCTGAc                                                                                                              <  6:37120/149‑1 (MQ=255)
                                           ttCACCGGCGGCACCACCACCGGCAAGAAAGTCATGGCCAGCGCCT‑CCAGCTCTTCGCTGAAGGAAGTGACCATGGAGCTGGGCGGCAAGTCGCCACTGATCATCTGTGCCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCatgat                                                                                                   >  1:171992/1‑149 (MQ=255)
                                            tCACCGGCGGCACCACCACCGGCAAGAAAGTCATGGCCAGCGCCTGCCAGCTCTTCGCTGAAGGAAGTGACCATGGAGCTGGGCGGCAAGTCGCCACTGATCATCTGTGCCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCatga                                                                                                    >  1:245003/1‑148 (MQ=255)
                                             cACCGGCGGCACCACCACCGGCAAGAAAGTCATGGCCAGCGCCT‑CCAGCTCTTCGCTGAAGGAAGTGACCATGGAGCTGGGCGGCAAGTCGCCACTGATCATCTGTGCCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCatgatg                                                                                                  <  1:264341/148‑1 (MQ=255)
                                             cACCGGCGGCACCACCACCGGCAAGAAAGTCATGGCCAGCGCCT‑CCAGCTCTTCGCTGAAGGAAGTGACCATGGAGCTGGGCGGCAAGTCGCCACTGATCATCTGTGCCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCatgat                                                                                                   <  5:275801/147‑1 (MQ=255)
                                              aCCGGCGGCACCACCACCGGCAAGAAAGTCATGGCCAGCGCCT‑CCAGCTCTTCGCTGAAGGAAGTGACCATGGAGCTGGGCGGCAAGTCGCCACTGATCATCTGTGCCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCATGATGGc                                                                                                >  8:95795/1‑149 (MQ=255)
                                                   cggcACCACCCCCGGCAATTAAGTCATGGCCAGCGCCT‑CCTGGTCTTCGCTGTAGGAAGTGACCATTGAGCTGGGGGGGAAGTCGCCACTGATGATCTGTGCCGACGCCGACCTGGTAGAAGGCCGCTGACATCGCCATGATGGCCAAc                                                                                            <  1:211874/149‑1 (MQ=255)
                                                       accaccaccGGCAAGAAAGTCATGGCCAGCGCCT‑CCAGCTCTTCGCTGAAGGAAGTGACCATGGAGCTGGGCGGCAAGTCGCCACTGATCATCTGTGCCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCATGATGGCCAACTTCt                                                                                        >  5:190865/1‑148 (MQ=255)
                                                                 gCAAGAAAGTCATGGCCAGCGCCT‑CCAGCTCTTCGCTGAAGGAAGTGACCATGGAGCTGGGCGGCAAGTCGCCACTGATCATCTGTGCCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCATGATGGCCAACTTCTACAGCTCgg                                                                               <  8:259243/147‑1 (MQ=255)
                                                                 gCAAGAAAGTCATGGCCAGCGCCT‑CCAGCTCTTCGCTGAAGGAAGTGACCATGGAGCTGGGCGGCAAGTCGCCACTGATCATCTGTGCCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCATGATGGCCAACTTCTACAGCTCgg                                                                               >  7:259243/1‑147 (MQ=255)
                                                                    aGAAAGTCATGGCCAGCGCCT‑CCAGCTCTTCGCTGAAGGAAGTGACCATGGAGCTGGGCGGCAAGTCGCCACTGATCATCTGTGCCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCATGATGGCCAACTTCTACAGCTCGGGCCAg                                                                          <  4:163390/149‑1 (MQ=255)
                                                                     gAAAGTCATGGCCAGCGCCT‑CCAGCTCTTCGCTGAAGGAAGTGACCATGGAGCTGGGCGGCAAGTCGCCACTGATCATCTGTGCCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCATGATGGCCAACTTc                                                                                         >  7:160567/1‑133 (MQ=255)
                                                                     gAAAGTCATGGCCAGCGCCT‑CCAGCTCTTCGCTGAAGGAAGTGACCATGGAGCTGGGCGGCAAGTCGCCACTGATCATCTGTGCCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCATGATGGCCAACTTc                                                                                         <  8:160567/133‑1 (MQ=255)
                                                                       aaGTCATGGCCAGCGCCT‑CCAGCTCTTCGCTGAAGGAAGTGACCATGGAGCTGGGCGGCAAGTCGCCACTGATCATCTGTGCCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCATGATGGCCAACTTCTACAGCTCGGGCCAGgtg                                                                       >  5:307124/1‑149 (MQ=255)
                                                                       aaGTCATGGCCAGCGCCT‑CCAGCTCTTCGCTGAAGGAAGTGACCATGGAGCTGGGCGGCAAGTCGCCACTGATCATCTGTGCCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCATGATGGCCAACTTCTACAGCTCGGGCCAGgt                                                                        <  1:242311/148‑1 (MQ=255)
                                                                          tCATGGCCAGCGCCT‑CCAGCTCTTCGCTGAAGGAAGTGACCATGGAGCTGGGCGGCAAGTCGCCACTGATCATCTGTGCCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCATGATGGCCAACTTCTACAGCTCGGGCCAGGTGTGc                                                                    <  2:171992/149‑1 (MQ=255)
                                                                              ggCCAGCGCCT‑CCAGCTCTTCGATGAAGGAAGTGAACATGGAGCTGGGCGGCAAGTCGCCACTGATCATCTGTGCCGCCGCCGACCTGG‑ACAAGGCCGCTGACATCGCCATGATGGCCAACTTCTACAGCTCGGGCCAGGTGTGAAca                                                                 >  7:149401/1‑147 (MQ=255)
                                                                               gCCAGCGCCT‑CCAGCTCTTCGCTGAAGGAAGTGACCATGGAGCTGGGCGGCAAGTCGCCACTGATCATCTGTGCCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCATGATGGCCAACTTCTACAGCTCGGGCCAGGTGTGCACCaa                                                               >  4:77659/1‑149 (MQ=255)
                                                                               gCCAGCGCCT‑CCAGCTCTTCGCTGAAGGAAGTGACCATGGAGCTGGGCGGCAAGTCGCCACTGATCATCTGTGCCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCATGATGGCCAACTTCTACAGCTCGGGCCAGGTGTGCACCaa                                                               <  7:95795/149‑1 (MQ=255)
                                                                                ccAGCGCCT‑CCAGCTCTTCGCTGAAGGAAGTGACCATGGAGCTGGGCGGCAAGTCGCCACTGATCATCTGTGCCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCATGATGGCCAACTTCTACAGCTCGGGCCAGGTGTGCACCAAc                                                              <  6:307124/149‑1 (MQ=255)
                                                                                    cgcCT‑CCAGCTCTTCGCTGAAGGAAGTGACCATGGAGCTGGGCGGCAAGTCGCCACTGATCATCTGTGCCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCATGATGGCCAACTTCTACAGCTCGGGCCAGGTGTGCACCAACGGCa                                                          >  1:115859/1‑149 (MQ=255)
                                                                                      ccT‑CCAGCTCTTCGCTGAAGGAAGTGACCATGGAGCTGGGCGGCAAGTCGCCACTGATCATCTGTGCCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCATGATGGCCAACTTCTACAGCTCGGGCCAGGTGTGc                                                                    >  4:249560/1‑137 (MQ=255)
                                                                                      ccT‑CCAGCTCTTCGCTGAAGGAAGTGACCATGGAGCTGGGCGGCAAGTCGCCACTGATCATCTGTGCCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCATGATGGCCAACTTCTACAGCTCGGGCCAGGTGTGc                                                                    <  3:249560/137‑1 (MQ=255)
                                                                                           cAGCTCTTCGCTGAAGGAAGTGACCATGGAGCTGGGCGGCAAGTCGCCACTGATCATCTGTGCCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCATGATAGCCAACTTCTACAGCTCGGGCCAGGTGTGCACCAACGGCACCcgcg                                                    <  2:245003/149‑1 (MQ=255)
                                                                                                 ttCGCTGAAGGAAGTGACCATGGAGCTGGGCGGCAAGTCGCCACTGATCATCTGTGCCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCATGATGGCCAACTTCTACAGCTCGGGCCAGGTGTGCACCAACGGCACCCGCGTGTTCa                                              >  1:69091/1‑149 (MQ=255)
                                                                                                         aggaagTGACCATGGAGCTGGGCGGCAAGTCGCCACTGATCATCTGTGCCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCATGATGGCCAACTTCTACAGCTCGGGCCAGGTGTGCACCAACGGCACCCGCGTGTTCATCCCGGCt                                      <  1:114428/149‑1 (MQ=255)
                                                                                                           gaagTGACCATGGAGCTGGGCGGCAAGTCGCCACTGATCATCTGTGCCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCATGATGGCCAACTTCTACAGCTCGGGCCAGGTGTGCACCAACGGCACCCGCGTGTTCATCCCGGCTGa                                    <  2:115859/149‑1 (MQ=255)
                                                                                                            aagTGACCATGGAGCTGGGCGGCAAGTCGCCACTGATCATCTGTGCCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCATGATGGCCAACTTCTACAGCTCGGGCCAGGTGTGCACCAACGGCACCCGCGTGTTCATCCCGGCTGaa                                   <  3:77659/149‑1 (MQ=255)
                                                                                                                gACCATGGAGCTGGGCGGCAAGTCGCCACTGATCATCTGTGCCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCATGATGGCCAACTTCTACAGCTCGGGCCAGGTGTGCACCAACGGCACCCGCGTGTTCATCCCGGCTGAAATg                                >  8:296666/1‑148 (MQ=255)
                                                                                                                  ccATGGAGCTGGGCGGCAAGTCGCCACTGATCATCTGTGCCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCATGATGGCCAACTTCTACAGCTCGGGCCAGGTGTGCACCAACGGCACCCGCGTGTTCATCCCGGCTGAAATGAAg                             <  7:296666/149‑1 (MQ=255)
                                                                                                                   cATGGAGCTCGGCGGCAAGGCGCCACTGATCAGCTGTGCCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCATGATGGCCAACTTCTACAGCTCGGGCCAGGTGTGCACCAACGGCACCCGCGTGTTCATCCCGGCTGAAATGAAgg                            <  8:172613/149‑1 (MQ=255)
                                                                                                                         gCTGGGCGGCAAGTCGCCACTGATCATCTGTGCCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCATGATGGCCAACTTCTACAGCTCGGGCCAGGTGTGCACCAACGGCACCCGCGTGTTCATCCCGGCTGAAATGAAGGCGGCCt                      >  5:31981/1‑149 (MQ=255)
                                                                                                                              gcggcAAGTCGCCACTGATCATCTGTGCCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCATGATGGCCAACTTCTACAGCTCGGGCCAGGTGTGCACCAACCGCACCCGCGTGTTCATCCCGGCTGAAATGAAGGCGGCCTTCGaa                 <  6:31981/149‑1 (MQ=255)
                                                                                                                                      tCGCCACTGATCATCTGTGCCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCATGATGGCCAACTTCTACAGCTCGGGCCAGGTGTGCACCAACGGCACCCGCGTGTTCATCCCGGCTGAAATGAAGGCGGCCTTCGAAGcc              >  7:289718/1‑144 (MQ=255)
                                                                                                                                        ggcacaGATGATCAGTGGCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCATGATGGCCAACTTCTAc                                                                                      <  1:241611/68‑1 (MQ=37)
                                                                                                                                        ggcacaGATGATCAGTGGCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCATGATGGCCAACTTCTAc                                                                                      >  2:241611/3‑70 (MQ=37)
                                                                                                                                         ccACTGATCATCTGTGCCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCATGATGGCCAACTTCTACAGCTCGGGCCAGGTGTGCACCAACGGCACCCGCGTGTTCATCCCGGCTGAAATGAAGGCGGCCTTCGAAGCCAAGATCGc      >  6:157309/1‑149 (MQ=255)
                                                                                                                                          cacaGATGATCAGTGGCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCATGATGGCCAACTTCTACa                                                                                     >  1:83527/1‑69 (MQ=38)
                                                                                                                                          cacaGATGATCAGTGGCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCATGATGGCCAACTTCTACa                                                                                     <  2:83527/69‑1 (MQ=38)
                                                                                                                                          cACTGATCATCTGTGCCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCATGATGGCCAACTTCTACAGCTCGGGCCAGGTGTGCACCACCGGCACCCGCGTGTTCATCCCGGCTGAAATGAAGGCGGCCTTCGAAGCCAAGATCGcc     >  6:327954/1‑149 (MQ=255)
                                                                                                                                          cACTGATCATCTGTGCCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCATGATGGCCAACTTCTACAGCTCGGGCCAGGTGTGCACCAACGGCACCCGCGTGTTCATCCCGGCTGAAATGAAGGCGGCCTTCGAAGCCAAGATCGcc     >  6:332798/1‑149 (MQ=255)
                                                                                                                                             tGATCATCTGTGCCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCATGATGGCCAACTTCTACAGCTCGGGCCAGGTGTGCACCAACGGCACCCGCGTGTTCATCCCGGCTGAAATGAAGGCGGCCTTCGAAGCCAAGATCGCCGAg  >  5:82471/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GGCACCTGGCTGACCGAGCACCCGCGCATCGAAAAAGTGTCC‑TTCACCGGCGGCACCACCACCGGCAAGAAAGTCATGGCCAGCGCCT‑CCAGCTCTTCGCTGAAGGAAGTGACCATGGAGCTGGGCGGCAAGTCGCCACTGATCATCTGTGCCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCATGATGGCCAACTTCTACAGCTCGGGCCAGGTGTGCACCAACGGCACCCGCGTGTTCATCCCGGCTGAAATGAAGGCGGCCTTCGAAGCCAAGATCGCCGAG  >  NC_002947/5774770‑5775057

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