Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	331,199	A→G	G579G (GGA→GGG)	betT →	choline:H(+) symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	331,199	0	A	G	100.0%	92.1 / NA	28	G579G (GGA→GGG)	betT	choline:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (13/15);  total (13/15)

ATGAATTATCCGGGCACGCGTTACACTAAACAGATGATGGAGACGGTCTGTTACCCGGCAATGGAAGAAGTGGCGCAGGAGTTGCGGTTGCGCGGCGCGTACGTGGAGCTAAAAAGCCTGCCACCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACC  >  NC_000913/331065‑331332                                                                                                                                       |                                                                                                                                      atgaATTATCCGGGCACGCGTTACACTAAACAGATGATGGAGACGGTCTGTTACCCGGCAATGGAAGAAGTGGCGCAGGAGTTGCGGTTGCGCGGCGCGTACGTGGAGCTAAAAAGCCTGCCACCGGAAGAGGGgcagca                                                                                                                                  >  2:370793/1‑140 (MQ=255)         tCCGGGCACGCGTTACACTAAACAGATGATGGAGACGGTCTGTTACCCGGCAATGGAAGAAGTGGCGCAGGAGTTGCGGTTGCGCGGCGCGTACGTGGAGCTAAAAAGCCTGCCACCGGAAGAGGGGCAGCAGTTGGGTc                                                                                                                          >  2:48127/1‑140 (MQ=255)           cGGGCACGCGTTACACTAAACAGATGATGGAGACGGTCTGTTACCCGGCAATGGAAGAAGTGGCGCAGGAGTTGCGGTTGCGCGGCGCGTACGTGGAGCTAAAAAGCCTGCCACCGGAAGAGGGGCAGCAGTTGGGTCAt                                                                                                                        >  2:596965/1‑140 (MQ=255)                              aCAGATGATGGAGACGGTCTGTTACCCGGCAATGGAAGAAGTGGCGCAGGAGTTGCGGTTGCGCGGCGCGTACGTGGAGCTAAAAAGCCTGCCACCGGAAGAGGGGCAGCAGTTGGGTCATCTGGATTTGTTGGTGCata                                                                                                     >  2:612774/1‑140 (MQ=255)                              aCAGATGATGGAGACGGTCTGTTACCCGGCAATGGAAGAAGTGGCGCAGGAGTTGCGGTTGCGCGGCGCGTACGTGGAGCTAAAAAGCCTGCCACCGGAACAGGGGCACCAGTTGGGTCATCTGGATTTGTTGGTGCata                                                                                                     >  2:274376/1‑140 (MQ=255)                                  atgatgGAGACGGTCTGTTACCCGGCAATGGAAGAAGTGGCGCAGGAGTTGCGGTTGCGCGGCGCGTACGTGGAGCTAAAAAGCCTGCCACCGGAAGAGGGGCAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATggg                                                                                                 >  2:478869/1‑140 (MQ=255)                                    gatgGAGACGGTCTGTTACCCGGCAATGGAAGAAGTGGCGCAGGAGTTGCGGTTGCGCGGCGCGTACGTGGAGCTAAAAAGCCTGCCACCGGAAGAGGGGCAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCg                                                                                               <  2:31578/140‑1 (MQ=255)                                    gatgGAGACGGTCTGTTACCCGGCAATGGAAGAAGTGGCGCAGGAGTTGCGGTTGCGCGGCGCGTACGTGGAGCTAAAAAGCCTGCCACCGGAAGAGGGGCAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCg                                                                                               <  2:635747/140‑1 (MQ=255)                                      tgGAGACGGTCTGTTACCCGGCAATGGAAGACGTGCCGCCGGAGTTGCTGTTGCGCGGCGCGTACGTGGAGCTAAAAAGCCTGCCACCGGAAGAGGGGCAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGaa                                                                                             <  2:184231/140‑1 (MQ=255)                                           aCGGTCTGTTACCCGGCAATGGAAGAAGTGGCGCAGGAGTTGCGGTTGCGCGGCGCGTACGTGGAGCTAAAAAGCCTGCCACCGGAAGAGGGGCAGCAGTGGGGCCATCTAGATGTGGTGGTGCATATGGGCGAAGAGCa                                                                                        >  1:185022/1‑140 (MQ=255)                                               tCTGTTACCCGGCAATGGAAGAAGTGGCGCAGGAGTTGCGGTTGCGCGGCGCGTACGTGGAGCTAAAAAGCCTGCCACCGGAAGAGGGGCAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAAc                                                                                    >  2:363556/1‑140 (MQ=255)                                                   ttACCCGGCAATGGAAGAAGTGGCGCAGGAGTTGCGGTTGCGCGGCGCGTACGTGGAGCTAAAAAGCCTGCCACCGGAAGAGGGGCAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTg                                                                                <  1:80950/140‑1 (MQ=255)                                                        cGGCAATGGAAGAAGTGGCGCAGGAGTTGCGGTTGCGCGGCGCGTACGTGGAGCTAAAAAGCCTGCCACCGGAAGAGGGGCAGCAGTTGGGTCATCTGGATTCGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTAt                                                                           <  2:42982/140‑1 (MQ=255)                                                         ggCAATGGAAGAAGTGGCGCAGGAGTTGCGGTTGCGCGGCGCGTACGTGGAGCTAAAAAGCCTGCCACCGGAAGAGGGGCAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATc                                                                          >  2:589685/1‑140 (MQ=255)                                                         ggCAATGGAAGAAGTGGCGCAGGAGTTGCGGTTGCGCGGCGCGTACGTGGAGCTAAAAAGCCTGCCACCAGAAGAGGGGCAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACttt                                                                                 <  2:560126/133‑1 (MQ=255)                                                         ggCAATGGAAGAAGTGGCGCAGGAGTTGCGGTTGCGCGGCGCGTACGTGGAGCTAAAAAGCCTGCCACCAGAAGAGGGGCAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACttt                                                                                 >  1:560126/1‑133 (MQ=255)                                                           cAATGGAAGAAGTGGCGCAGGAGTTGCGGTTGCGCGGCGCGTACGTGGAGCTAAAAAGCCTGCCACCGGAAGAGGGGCAGCAGTTGGGTCATCTGGATTTGTTGGTGCATCTGGGCGGAGAGCAAAACTTTGTCTATCAg                                                                        >  1:321567/1‑140 (MQ=255)                                                                  agaagTGGCGCAGGAGTTGCGGTTGCGCGGCGCGTACGTGGAGCTAAAAAGCCTGCCACCGGAAGAGGGGCAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGc                                                                 <  1:363556/140‑1 (MQ=255)                                                                    aagTGGCGCAGGAGTTGGGGTTGCGCGGCGCCTACTTGGTGCTAAAAAGCCTGCCACCGGAAGAGGGGCAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCg                                                               <  2:253947/140‑1 (MQ=255)                                                                               gAGTTGCGGTTGCGCGGCGCGTACGTGGAGCTAAAAAGCCTGCCACCGGAAGAGGGGCAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTc                                                    >  1:563301/1‑140 (MQ=255)                                                                                     cggttgcgCGGCGCGTACGTGGAGCTAAAAAGCCTGCCACCGGAAGAGGGGCAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGcc                                              <  1:344010/140‑1 (MQ=255)                                                                                         tgcgCGGCGCGTACGTGGAGCTAAAAAGCCTGCCACCGGAAGAGGGGCAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGcc                                              <  2:570767/136‑1 (MQ=255)                                                                                         tgcgCGGCGCGTACGTGGAGCTAAAAAGCCTGCCACCGGAAGAGGGGCAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGcc                                              >  1:570767/1‑136 (MQ=255)                                                                                         tgcgCGGCGCGTACGTGGAGCTAAAAAGCCTGCCACCGGAAGAGGGGCAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGc                                          <  1:370793/140‑1 (MQ=255)                                                                                                      cGTGGAGCTAAAAAGCCTGCCACCGGAAGAGGGGCAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACcgcg                             <  1:478869/140‑1 (MQ=255)                                                                                                             cTAAAAAGCCTGCCACCGGAAGAGGGGCAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAg                      >  2:348979/1‑140 (MQ=255)                                                                                                              tAAAAAGCCTGCCACCGGAAGAGGGGCAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGc                     <  2:522136/140‑1 (MQ=255)                                                                                                                  aaGCCTGCCACCGGAAGAGGGGCAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTa                 <  2:563301/140‑1 (MQ=255)                                                                                                                                 agagGGGCAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTAcc  <  1:348979/140‑1 (MQ=255)                                                                                                                                       |                                                                                                                                      ATGAATTATCCGGGCACGCGTTACACTAAACAGATGATGGAGACGGTCTGTTACCCGGCAATGGAAGAAGTGGCGCAGGAGTTGCGGTTGCGCGGCGCGTACGTGGAGCTAAAAAGCCTGCCACCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACC  >  NC_000913/331065‑331332

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