Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	328,682	C→T	A18T (GCC→ACC)	betB ←	betaine aldehyde dehydrogenase, NAD‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	328,682	0	C	T	100.0%	165.3 / NA	51	A18T (GCC→ACC)	betB	betaine aldehyde dehydrogenase, NAD‑dependent
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/29);  total (22/29)

CCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTT  >  NC_000913/328554‑328806                                                                                                                                 |                                                                                                                             ccAGATTTTTCGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGTGGAGGTATaa                                                                                                                    <  2:1234639/139‑1 (MQ=255)   aGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAGGGTGCGACCGCTGGTGGTGGAGGTATAacc                                                                                                                  <  2:191777/139‑1 (MQ=255)     aTTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGTGGAGGTATAaccac                                                                                                                <  2:433781/139‑1 (MQ=255)      ttttttGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGTGGAGGTATAaccacc                                                                                                               <  1:1267263/139‑1 (MQ=255)       tttttGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGTGGAGGTATAAccacca                                                                                                              <  2:130518/139‑1 (MQ=255)        ttttGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGTGGAGGTATAACCACCAt                                                                                                             <  1:1208649/139‑1 (MQ=255)        ttttGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGTGGAGGTATAACCACCAt                                                                                                             <  2:48428/139‑1 (MQ=255)        ttttGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGTGGAGGTATAACCACCAt                                                                                                             <  1:289229/139‑1 (MQ=255)        ttttGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGTGGAGGTATAACCACCAt                                                                                                             <  2:1311278/139‑1 (MQ=255)           tGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGTGGAGGTATAaccacc                                                                                                               >  2:783855/1‑134 (MQ=255)           tGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGTGGAGGTATAaccacc                                                                                                               <  1:783855/134‑1 (MQ=255)               ccTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGTGGAGGTATAacc                                                                                                                  >  2:634154/1‑127 (MQ=255)               ccTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGTGGAGGTATAacc                                                                                                                  <  1:634154/127‑1 (MQ=255)                           ttttCACGGCGCGATCGACATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGTGGa                                                                                                                           >  1:1013742/1‑106 (MQ=255)                           ttttCACGGCGCGATCGACATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGTGGa                                                                                                                           <  2:1013742/106‑1 (MQ=255)                             ttCACGGCGCGATCGACATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGTGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTg                                                                                        <  2:330826/139‑1 (MQ=255)                                             cATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGTGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCAc                                                                        <  2:814159/139‑1 (MQ=255)                                              aTCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGTGGAGGTATAAccacca                                                                                                              <  1:1120496/100‑1 (MQ=255)                                              aTCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGTGGAGGTATAAccacca                                                                                                              >  2:1120496/1‑100 (MQ=255)                                              aTCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGTGGAGGTATAACCACCAt                                                                                                             <  2:1580092/101‑1 (MQ=255)                                              aTCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGTGGAGGTATAACCACCAt                                                                                                             >  1:1580092/1‑101 (MQ=255)                                                   cgcgCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGTGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTctcc                                                                  >  2:942302/1‑139 (MQ=255)                                                     cgcCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGTGGAGGtat                                                                                                                      <  2:477655/85‑1 (MQ=255)                                                     cgcCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGTGGAGGtat                                                                                                                      >  1:477655/1‑85 (MQ=255)                                                        ccGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGTGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGt                                                             <  1:305417/139‑1 (MQ=255)                                                          ggcggcCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGTGGAGGTATAaccacc                                                                                                               <  1:260598/87‑1 (MQ=255)                                                          ggcggcCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGTGGAGGTATAaccacc                                                                                                               >  2:260598/1‑87 (MQ=255)                                                           gcggcCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGTGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCa                                                                         <  1:1311654/124‑1 (MQ=255)                                                           gcggcCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGTGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCa                                                                         >  2:1311654/1‑124 (MQ=255)                                                                        tCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGTGGAGGTATAacc                                                                                                                  <  1:844878/70‑1 (MQ=255)                                                                        tCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGTGGAGGTATAacc                                                                                                                  >  2:844878/1‑70 (MQ=255)                                                                                            gCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGTGGAGGTATAAccacca                                                                                                              >  2:875121/1‑54 (MQ=255)                                                                                            gCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGTGGAGGTATAAccacca                                                                                                              <  1:875121/54‑1 (MQ=255)                                                                                            gCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGTGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGAc                                                                             <  1:1551526/87‑1 (MQ=255)                                                                                            gCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGTGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGAc                                                                             >  2:1551526/1‑87 (MQ=255)                                                                                            gCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGTGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCt                                          >  1:1016618/1‑122 (MQ=255)                                                                                            gCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGTGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCt                                          <  2:1016618/122‑1 (MQ=255)                                                                                                ggTTAATGGTCTCGAAGGTGCGACCGCTGGTGGTGGAGGTATAACCACCAt                                                                                                             <  2:1361770/51‑1 (MQ=255)                                                                                                ggTTAATGGTCTCGAAGGTGCGACCGCTGGTGGTGGAGGTATAACCACCAt                                                                                                             >  1:1361770/1‑51 (MQ=255)                                                                                                 gTTAATGGTCTCGAAGGTGCGACCGCTGGTGGTGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCgg                                                      <  2:551643/105‑1 (MQ=255)                                                                                                 gTTAATGGTCTCGAAGGTGCGACCGCTGGTGGTGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCgg                                                      >  1:551643/1‑105 (MQ=255)                                                                                                           tcGAAGGTGCGACCGCTGGTGGTGGAGGTATAAccacca                                                                                                              <  1:84299/39‑1 (MQ=255)                                                                                                           tcGAAGGTGCGACCGCTGGTGGTGGAGGTATAAccacca                                                                                                              >  2:84299/1‑39 (MQ=255)                                                                                                           tcGAAGGTGCGACCGCTGGTGGTGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAAt                                                         >  2:1570014/1‑92 (MQ=255)                                                                                                           tcGAAGGTGCGACCGCTGGTGGTGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAAt                                                         <  1:1570014/92‑1 (MQ=255)                                                                                                           tcGAAGGTGCGACCGCTGGTGGTGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAgcg          >  2:788273/1‑139 (MQ=255)                                                                                                           tcGAAGGTGCGACCGCTGGTGGTGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAgcg          >  2:126342/1‑139 (MQ=255)                                                                                                           tcGAAGGTGCGACCGCTGGTGGTGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAgcg          >  2:914562/1‑139 (MQ=255)                                                                                                           tcGAAGGTGCGACCGCTGGTGGTGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAgcg          <  1:942302/139‑1 (MQ=255)                                                                                                                  tGCGACCGCTGGTGGTGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGgtt   >  1:1612417/1‑139 (MQ=255)                                                                                                                   gCGACCGCTGGTGGTGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGtttt  >  2:214523/1‑139 (MQ=255)                                                                                                                                 |                                                                                                                             CCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTT  >  NC_000913/328554‑328806

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