Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,774,917	T→A	8.3%	C262S (TGT→AGT)	betB →	NAD‑dependent betaine aldehyde dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,774,917	0	T	A	8.3%	134.3 / 2.9	48	C262S (TGT→AGT)	betB	NAD‑dependent betaine aldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (23/21);  new base A (2/2);  total (25/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.85e-01

GGCACCTGGCTGACCGAGCACCCGCGCATCGAAAAAGTGTCC‑TTCACCGGCGGCACCACCACCGGCAAGAAAGTCATGGCCAGCGCCT‑CCAGCTCTTCGCTGAAGGAAGTGACCATGGAGCTGGGCGGCAAGTCGCCACTGATCATCTGTGCCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCATGATGGCCAACTTCTACAGCTCGGGCCAGGTGTGCACCAACGGCACCCGCGTGTTCATCCCGGCTGAAATGAAGGCGGCCTTCGAAGCCAAGATCGCCGAGCGCGTTGC  >  NC_002947/5774770‑5775065                                                                                                                                                      |                                                                                                                                                      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)         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                                                                                       >  7:259243/1‑147 (MQ=255)                                                                  gCAAGAAAGTCATGGCCAGCGCCT‑CCAGCTCTTCGCTGAAGGAAGTGACCATGGAGCTGGGCGGCAAGTCGCCACTGATCATCTGTGCCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCATGATGGCCAACTTCTACAGCTCgg                                                                                       <  8:259243/147‑1 (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)                                                                                                                                                      tgtgCCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCATGATGGCCAACTTCTACAGCTCGGGCCAGGTGTGCACCAACGGCACCCGCGTGTTCATCCCGGCTGAAATGAAGGCGGCCTTCGAAGCCAAGATCGCCGAGCGCGTTGc  >  2:241666/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                      GGCACCTGGCTGACCGAGCACCCGCGCATCGAAAAAGTGTCC‑TTCACCGGCGGCACCACCACCGGCAAGAAAGTCATGGCCAGCGCCT‑CCAGCTCTTCGCTGAAGGAAGTGACCATGGAGCTGGGCGGCAAGTCGCCACTGATCATCTGTGCCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCATGATGGCCAACTTCTACAGCTCGGGCCAGGTGTGCACCAACGGCACCCGCGTGTTCATCCCGGCTGAAATGAAGGCGGCCTTCGAAGCCAAGATCGCCGAGCGCGTTGC  >  NC_002947/5774770‑5775065

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