Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,774,921	C→G	8.7%	A263G (GCC→GGC)	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,921	0	C	G	8.7%	142.6 / 3.9	46	A263G (GCC→GGC)	betB	NAD‑dependent betaine aldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (22/20);  new base G (2/2);  total (24/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.97e-01

GCTGACCGAGCACCCGCGCATCGAAAAAGTGTCC‑TTCACCGGCGGCACCACCACCGGCAAGAAAGTCATGGCCAGCGCCT‑CCAGCTCTTCGCTGAAGGAAGTGACCATGGAGCTGGGCGGCAAGTCGCCACTGATCATCTGTGCCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCATGATGGCCAACTTCTACAGCTCGGGCCAGGTGTGCACCAACGGCACCCGCGTGTTCATCCCGGCTGAAATGAAGGCGGCCTTCGAAGCCAAGATCGCCGAGCGCGTTGCC  >  NC_002947/5774778‑5775066                                                                                                                                                  |                                                                                                                                                   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                                                                                                  <  8:160567/133‑1 (MQ=255)                                                              gAAAGTCATGGCCAGCGCCT‑CCAGCTCTTCGCTGAAGGAAGTGACCATGGAGCTGGGCGGCAAGTCGCCACTGATCATCTGTGCCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCATGATGGCCAACTTc                                                                                                  >  7:160567/1‑133 (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                                                                             <  3:249560/137‑1 (MQ=255)                                                                               ccT‑CCAGCTCTTCGCTGAAGGAAGTGACCATGGAGCTGGGCGGCAAGTCGCCACTGATCATCTGTGCCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCATGATGGCCAACTTCTACAGCTCGGGCCAGGTGTGc                                                                             >  4:249560/1‑137 (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)                                                                                                                                                  ccgacgccgacCTGG‑ACAAGGCCGCTGACATCGCCATGATGGCCAACTTCTACAGCTCGGGCCAGGTGTGCACCAACGGCACCCGCGTGTTCATCCCGGCTGAAATGAAGGCGGCCTTCGAAGCCAAGATCGCCGAGCGCGTTGcc  <  6:211060/146‑1 (MQ=255)                                                                                                                                                  ccgacgccgacCTGG‑ACAAGGCCGCTGACATCGCCATGATGGCCAACTTCTACAGCTCGGGCCAGGTGTGCACCAACGGCACCCGCGTGTTCATCCCGGCTGAAATGAAGGCGGCCTTCGAAGCCAAGATCGCCGAGCGCGTTGcc  >  5:211060/1‑146 (MQ=255)                                                                                                                                                  |                                                                                                                                                   GCTGACCGAGCACCCGCGCATCGAAAAAGTGTCC‑TTCACCGGCGGCACCACCACCGGCAAGAAAGTCATGGCCAGCGCCT‑CCAGCTCTTCGCTGAAGGAAGTGACCATGGAGCTGGGCGGCAAGTCGCCACTGATCATCTGTGCCGACGCCGACCTGG‑ACAAGGCCGCTGACATCGCCATGATGGCCAACTTCTACAGCTCGGGCCAGGTGTGCACCAACGGCACCCGCGTGTTCATCCCGGCTGAAATGAAGGCGGCCTTCGAAGCCAAGATCGCCGAGCGCGTTGCC  >  NC_002947/5774778‑5775066

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