Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	331,300	A→G	Y613C (TAC→TGC)	betT →	choline:H(+) symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	331,300	0	A	G	100.0%	178.7 / NA	51	Y613C (TAC→TGC)	betT	choline:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/25);  total (26/25)

GAGCTAAAAAGCCTGCCACCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCG  >  NC_000913/331170‑331429                                                                                                                                   |                                                                                                                                  gAGCTAAAAAGCCTGCCACCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTGCCGcgcacg                                                                                                                          <  1:90150/140‑1 (MQ=255)      aaaaaGCCTGCCACCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAATAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTGCCGCGCACGCAGCg                                                                                                                     <  1:201284/140‑1 (MQ=255)      aaaaaGCCTGCCACCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTGCCGCGCACGCAGCg                                                                                                                     >  2:713269/1‑140 (MQ=255)        aaaGCCTGCCACCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTGCCGCGCACGCAGCGGt                                                                                                                   >  1:195942/1‑140 (MQ=255)            ccTGCCACCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTGCCGCGCACGCAGCGGTAAAt                                                                                                               <  2:788687/140‑1 (MQ=255)            ccTGCCACCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTGCCGCGCACGCAGCGGTAAAt                                                                                                               >  2:161298/1‑140 (MQ=255)            ccTGCCACCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTGCCGCGCACGCAGCGGTAAAt                                                                                                               >  2:175475/1‑140 (MQ=255)            ccTGCCACCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTGCCGCGCACGCAGCGGTAAAt                                                                                                               <  1:611321/140‑1 (MQ=255)            ccTGCCACCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTGCCGCGCACGCAGCGGTAAAt                                                                                                               <  1:161298/140‑1 (MQ=255)            ccTGCCACCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTGCCGCGCACGCAGCGGTAAAt                                                                                                               <  2:750031/140‑1 (MQ=255)               gCCACCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTGCCGcgcacgc                                                                                                                         >  1:503008/1‑127 (MQ=255)               gCCACCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTGCCGcgcacgc                                                                                                                         <  2:503008/127‑1 (MQ=255)                 cACCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTGCCGCGCACGCAGCGGTAAATCGAcc                                                                                                          <  2:701542/140‑1 (MQ=255)                      gAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTGCCGCGCACGCAGCGGTAAATc                                                                                                              >  2:57178/1‑131 (MQ=255)                      gAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTGCCGCGCACGCAGCGGTAAATc                                                                                                              <  1:57178/131‑1 (MQ=255)                                agcagTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTGCCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGaa                                                                                           >  2:337537/1‑140 (MQ=255)                                 gcagTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTGCCGCGCACGCAGCGGTAAATCGACCTACTACCgg                                                                                                 >  2:689100/1‑133 (MQ=255)                                 gcagTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTGCCGCGCACGCAGCGGTAAATCGACCTACTACCgg                                                                                                 <  1:689100/133‑1 (MQ=255)                                   agTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTGCCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAAcc                                                                                        <  1:570445/140‑1 (MQ=255)                                                ggATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTGCCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAg                                                                           <  1:14611/140‑1 (MQ=255)                                                  aTTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTGCCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTc                                                                                     <  1:217806/128‑1 (MQ=255)                                                  aTTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTGCCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTc                                                                                     >  2:217806/1‑128 (MQ=255)                                                      gttgGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTGCCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGcc                                                                     >  2:107047/1‑140 (MQ=255)                                                       ttgGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTGCCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCa                                                                    <  2:773254/140‑1 (MQ=255)                                                           tGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTGCCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGc                                                                >  2:737529/1‑140 (MQ=255)                                                                    gCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTGCCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTg                                                       >  1:356541/1‑140 (MQ=255)                                                                     cGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTGCCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGa                                                      >  2:342956/1‑140 (MQ=255)                                                                       aaGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTGCCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATg                                                    >  2:408509/1‑140 (MQ=255)                                                                         gagCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTGCCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGa                                                  >  1:229426/1‑140 (MQ=255)                                                                         gagCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTGCCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGa                                                  <  2:195942/140‑1 (MQ=255)                                                                         gagCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTGCCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGa                                                  <  2:496472/140‑1 (MQ=255)                                                                          agCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTGCCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGAc                                                 >  2:297093/1‑140 (MQ=255)                                                                            cAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTGCCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTa                                               >  2:156437/1‑140 (MQ=255)                                                                            cAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTGCCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTa                                               >  1:52973/1‑140 (MQ=255)                                                                              aaaCTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTGCCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACa                                             <  2:281111/140‑1 (MQ=255)                                                                              aaaCTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTGCCGCGCACGCAACGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACa                                             >  1:725030/1‑140 (MQ=255)                                                                                 cTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTGCCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCa                                          <  2:356541/140‑1 (MQ=255)                                                                                    tGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTGCCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAag                                       >  2:541631/1‑140 (MQ=255)                                                                                       cTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTGCCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGc                                    <  2:409858/140‑1 (MQ=255)                                                                                       cTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTGCCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGc                                    <  2:229426/140‑1 (MQ=255)                                                                                              aTTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTGCCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGAt                             <  1:337537/140‑1 (MQ=255)                                                                                              aTTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTGCCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGAt                             <  2:725030/140‑1 (MQ=255)                                                                                                        agcaATATTCGGTGCCGGGCTTTACCTGCCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGa                       <  2:544703/136‑1 (MQ=255)                                                                                                        agcaATATTCGGTGCCGGGCTTTACCTGCCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGa                       >  1:544703/1‑136 (MQ=255)                                                                                                              tttCGGTGCCGGGCTTTACCTGCCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGAc             <  2:685285/139‑1 (MQ=255)                                                                                                              atTCGGTGCCGGGCTTTACCTGCCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGAc             <  2:180565/140‑1 (MQ=255)                                                                                                                    tGCCGGGCTTTACCTGCCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGt         >  2:159279/1‑138 (MQ=255)                                                                                                                    tGCCGGGCTTTACCTGCCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGt         <  1:159279/138‑1 (MQ=255)                                                                                                                       cGGGCTTTACCTGCCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCGCCGATATTCTTGACCAGTACGAg    >  2:145922/1‑140 (MQ=255)                                                                                                                       cGGGCTTTACCTGCCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAg    >  2:72623/1‑140 (MQ=255)                                                                                                                         ggCTTTACCTGCCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCg  >  2:530654/1‑140 (MQ=255)                                                                                                                         ggCTTTACCTGCCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCg  >  1:132264/1‑140 (MQ=255)                                                                                                                                   |                                                                                                                                  GAGCTAAAAAGCCTGCCACCGGAAGAGGGACAGCAGTTGGGTCATCTGGATTTGTTGGTGCATATGGGCGAAGAGCAAAACTTTGTCTATCAGATTTGGCCGCAGCAATATTCGGTGCCGGGCTTTACCTACCGCGCACGCAGCGGTAAATCGACCTACTACCGGCTGGAAACCTTCCTGTTAGAAGGCAGCCAGGGCAACGACCTGATGGACTACAGCAAAGAGCAGGTGATCACCGATATTCTTGACCAGTACGAGCG  >  NC_000913/331170‑331429

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