Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	2,750,562	T→C	V197V (GTA→GTG)	betA ←	choline dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	2,750,562	0	T	C	100.0%	191.4 / NA	52	V197V (GTA→GTG)	betA	choline dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/36);  total (16/36)

TTTTTCTTATACGTAACACCAGTTGCTCTTCTACCTTCATAATGAATTTCAGTTACAAAGGCACGTGTTTCAACGGTTAAGTTTTTACGCTTCATCGCTGGATGTAAATATGCTCTTGAAGCTGACATTCGGCGACCACGATGTACTTGACTATCGAACGGTCCAAAACCTTCTTGTCTAAATCCATTCACATCAGGTGTTTTATGATAGCCTGCTTCAACACCTGCATCAAAGAATGACTGGAATAAAGGATTCGTTGCTGGCCCTCGTTTTAACTTAATT  >  CP000730/2750419‑2750700                                                                                                                                                |                                                                                                                                           tttttCTTATACGTAACACCAGTTGCTCTTCTACCTTCATAATGAATTTCAGTTACAAAGGCACGTGTTTCAACGGTTAAGTTTTTACGCTTCATCGCTGGATGTAAATATGCTCTTGAAGCTGACATTCGGCGGCCACGATGCACtt                                                                                                                                        <  2:382554/148‑1 (MQ=255) tttttCTTATACGTAACACCAGTTGCTCTTCTACCTTCATAATGAATTTCAGTTACAAAGGCACGTGTTTCAACGGTTAAGTTTTTACGCTTCATCGCTGGATGTAAATATGCTCTTGAAGCTGACATTCGGCGACCACGATGCACtt                                                                                                                                        <  2:382548/148‑1 (MQ=255) tttttCTTATACGTAACACCAGTTGCTCTTCTACCTTCATAATGAATTTCAGTTACAAAGGCACGTGTTTCAACGGTTAAGTTTTTACGCTTCATCGCTGGATGTAAATATGCTCTTGAAGCTGACATTCGGCGACCACGATGCACtt                                                                                                                                        <  2:258236/148‑1 (MQ=255) tttttCTTATACGTAACACCAGTTGCTCTTCTACCTTCATAATGAATTTCAGTTACAAAGGCACGTGTTTCAACGGTTAAGTTTTTACGCTTCATCGCTGGATGTAAATATGCTCTTGAAGCTGACATTCGGCGACCACGATGCACtt                                                                                                                                        >  1:21095/1‑148 (MQ=255)  ttttCTTATACGTAACACCAGTTGCTCTTCTACCTTCATAATGAATTTCAGTTACAAAGGCACGTGTTTCAACGGTTAAGTTTTTACGCTTCATCGCTGGATGTAAATATGCTCTTGAAGCTGACATTCGGCGACCACGATGCACTTg                                                                                                                                       >  2:351563/1‑148 (MQ=255)  ttttCTTATACGTAACACCAGTTGCTCTTCTACCTTCATAATGAATTTCAGTTACAAAGGCACGTGTTTCAACGGTTAAGTTTTTACGCTTCATCGCTGGATGTAAATATGCTCTTGAAGCTGACATTCGGCGACCACGATGCACTTg                                                                                                                                       >  2:220668/1‑148 (MQ=255)         ataCGTAACACCAGTTGCTCTTCTACCTTCATAATGAATTTCAGTTACAAAGGCACGTGTTTCAACGGTTAAGTTTTTACGCTTCATCGCTGGATGTAAATATGCTCTTGAAGCTGACATTCGGCGACCACGATGCACTTGACTATCg                                                                                                                                >  1:43097/1‑148 (MQ=255)         ataCGTAACACCAGTTGCTCTTCTACCTTCATAATGAATTTCAGTTACAAAGGCACGTGTTTCAACGGTTAAGTTTTTACGCTTCATCGCTGGATGTAAATATGCTCTTGAAGCTGACATTCGGCGACCACGATGCACTTGACTATCg                                                                                                                                >  1:9186/1‑148 (MQ=255)            cGTAACACCAGTTGCTCTTCTACCTTCATAATGAATTTCAGTTACAAAGGCACGTGTTTCCACGGTTAAGTTTTTACGCTTCATCGCTGGATGTAAATATGCTCTTGAAGCTGACATTCGGCGACCACGATGCACTTGACTATCGAAc                                                                                                                             >  2:205894/1‑148 (MQ=255)              tAACACCAGTTGCTCTTCTACCTTCATAATGAATTTCAGTTACAAAGGCACGTGTTTCAACGGTTAAGTTTTTACGCTTCATCGCTGGATGTAAATATGCTCTTGAAGCTGACATTCGGCGACCACGATGCACTTGACTATCGAACgg                                                                                                                           >  1:164766/1‑148 (MQ=255)              tAACACCAGTTGCTCTTCTACCTTCATAATGAATTTCAGTTACAAAGGCACGTGTTTCAACGGTTAAGTTTTTACGCTTCATCGCTGGATGTAAATATGCTCTTGAAGCTGACATTCGGCGACCACGATGCACTTGACTATCGAACgg                                                                                                                           >  2:121607/1‑148 (MQ=255)                  acCAGTTGCTCTTCTACCTTCATAATGAATTTCAGTTACAAAGGCACGTGTTTCAACGGTTAAGTTTTTACGCTTCATCGCTGGATGTAAATATGCTCTTGAAGCTGACATTCGGCGACCACGATGCACTTGACTATCGAACGGTCCa                                                                                                                       <  2:354463/148‑1 (MQ=255)                  acCAGTTGCTCTTCTACCTTCATAATGAATTTCAGTTACAAAGGCACGTGTTTCAACGGTTAAGTTTTTACGCTTCATCGCTGGATGTAAATATGCTCTTGAAGCTGACATTCGGCGACCACGATGCACTTGACTATCGAACGGTCCa                                                                                                                       <  2:403074/148‑1 (MQ=255)                       ttGCTCTTCTACCTTCATAATGAATTTCAGTTACAAAGGCACGTGTTTCAACGGTTAAGTTTTTACGCTTCATCGCTGGATGTAAATATGCTCTTGAAGCTGGCATTCGGCGACCACGATGCACTTGACTATCGAACGGTCCAAAAcc                                                                                                                  <  1:305536/148‑1 (MQ=255)                        tGCTCTTCTACCTTCATAATGAATTTCAGTTACAAAGGCACGTGTTTCAACGGTTAAGTTTTTACGCTTCATCGCTGGATGTAAATATGCTCTTGAAGCTGACATTCGGCGACCACGATGCACTTGACTATCGAACGGTCCAAAACct                                                                                                                 <  2:83685/148‑1 (MQ=255)                              tctACCTTCATAATGAATTTCAGTTACAAAGGCACGTGTTTCAACGGTTAAGTTTTTACGCTTCATCGCTGGATGTAAATATGCTCTTGAAGCTGACATTCGGCGACCACGATGCACTTGACTATCGAACGGTCCAAAACCTTCTTGt                                                                                                           <  2:145914/148‑1 (MQ=255)                                            gAATTTCAGTTACAAAGGCACGTGTTTCAACGGTTAAGTTTTTACTCTTCATCGCTGGATGTAAATATGCTCTTGAAGCTGACATTCGGCGACCACGATGCACTTGACTATCGAACGGTCCAAAACCTTCTTGTCTAAATCCATTcac                                                                                             >  1:219698/1‑148 (MQ=255)                                            gAATTTCAGTTACAAAGGCACGTGTTTCAACGGTTAAGTTTTTACGCTTCATCGCTGGATGTAAATATGCTCTTGAAGCTGACATTCGGCGACCACGATGCACTTGACTATCGAACGGTCCAAAACCTTCTTGTCTAAATCCATTcac                                                                                             >  1:115003/1‑148 (MQ=255)                                            gAATTTCAGTTACAAAGGCACGTGTTTCAACGGTTAAGTTTTTACGCTTCATCGCTGGATGTAAATATGCTCTTGAAGCTGACATTCGGCGACCACGATGCACTTGACTATCGAACGGTCCAAAACCTTCTTGTCTAAATCCATTcac                                                                                             >  1:398978/1‑148 (MQ=255)                                            gAATTTCAGTTACAAAGGCACGTGTTTCAACGGTTAAGTTTTTACGCTTCATCGCTGGATGTAAATAAAATCTTGAAGCTGACATTCGGCGACCACGATGCACTTGACTATCCAACGGGCCAAAACCTTCTTGTCTAAATCCATTcac                                                                                             >  2:373780/1‑148 (MQ=255)                                              aTTTCAGTTACAAAGGCACGTGTTTCAACGGTTAAGTTTTTACGCTTCATCGCTGGATGTAAATATGCTCTTGAAGCTGACATTCGGCGACCACGATGCACTTGACTATCGAACGGTCCAAAACCTTCTTGTCTAAATCCATTCACAt                                                                                           >  1:253302/1‑148 (MQ=255)                                                     ttACAAAGGCACGTGTTTCAACGGTTAAGTTTTTACGCTTCATCGCTGGATGTAAATATGCTCTTGAAGCTGACATTCGGCGACCACGATGCACTTGACTATCGAACGGTCCAAAACCTTCTTGTCTAAATCCATTCACATCAGgtgt                                                                                    <  1:57346/148‑1 (MQ=255)                                                     ttACAAAGGCACGTGTTTCAACGGTTAAGTTTTTACGCTTCATCGCTGGATGTAAATATGCTCTTGAAGCTGACATTCGGCGACCACGATGCACTTGACTATCGAACGGTCCAAAACCTTCTTGTCTAAATCCATTCACATCAGgtgt                                                                                    <  1:193694/148‑1 (MQ=255)                                                           aGGCACGTGTTTCAACGGTTAAGTTTTTACGCTTCATCGCTGGATGTAAATATGCTCTTGAAGCTGACATTCGGCGACCACGATGCACTTGACTATCGAACGGTCCAAAACCTTCTTGTCTAAATCCATTCACATCAGGTGTTTTATg                                                                              <  1:373990/148‑1 (MQ=255)                                                             gCACGTGTTTCAACGGTTAAGTTTTTACGCTTCATCGCTGGATGTAAATATGCTCTTGAAGCTGACATTCGGCGACCACGATGCACTTGACTATCGAACGGTCCAAAACCTTCTTGTCTAAATCCATTCACATCAGGTGTTTTATGAt                                                                            >  2:350894/1‑148 (MQ=255)                                                                   gtTTCAACGGTTAAGTTTTTACGCTTCATCGCTGGATGTAAATATGCTCTTGAAGCTGACATTCGGCGACCACGATGCACTTGACTATCGAACGGTCCAAAACCTTCTTGTCTAAATCCATTCACATCAGGTGTTTTATGATAGCCTg                                                                      >  1:102003/1‑148 (MQ=255)                                                                          cGGTTAAGTTTTTACGCTTCATCGCTGGATGTAAATATGCTCTTGAAGCTGACATTCGGCGACCACGATGCACTTGACTATCGAACGGTCCAAAACCTTCTTGTCTAAATCCATTCACATCAGGTGTTTTATGATAGCCTGCTTCAac                                                               <  2:376260/148‑1 (MQ=255)                                                                             ttAAGTTTTTACGCTTCATCGCTGGATGTAAATATGCTCTTGAAGCTGACATTCGGCGACCACGATGCACTTGACTATCGAACGGTCCAAAACCTTCTTGTCTAAATCCATTCACATCAGGTGTTTTATGATAGCCTGCTTCAACAcc                                                            <  2:167789/148‑1 (MQ=255)                                                                                 gTTTTTACGCTTCATCGCTGGATGTAAATATGCTCTTGAAGCTGACATTCGGCGTCCACGATGCACTTGACTATCGAACGGTCCAAAACCTTCTTGTCTAAATCCATTCACATCAGGTGTTTTATGATAGCCTGCTTCAACACCTGCa                                                        <  2:331318/148‑1 (MQ=255)                                                                                 gTTTTTACGCTTCATCGCTGGATGTAAATATGCTCTTGAAGCTGACATTCGGCGACCACGATGCACTTGACTATCGAACGGTCCAAAACCTTCTTGTCTAAATCCATTCACATCAGGTGTTTTATGATAGCCTGCTTCAACACCTGCa                                                        <  1:310858/148‑1 (MQ=255)                                                                                 gTTTTTACGCTTCATCGCTGGATGTAAATATGCTCTTGAAGCTGACATTCGGCGACCACGATGCACTTGACTATCGAACGGTCCAAAACCTTCTTGTCTAAATCCATTCACATCAGGTGTTTTATGATAGCCTGCTTCAACACCTGCa                                                        <  2:129160/148‑1 (MQ=255)                                                                                 gTTTTTACGCTTCATCGCTGGATGTAAATATGCTCTTGAAGCTGACATTCGGCGACCACGATGCACTTGACTATCGAACGGTCCAAAACCTTCTTGTCTAAATCCATTCACATCAGGTGTTTTATGATAGCCTGCTTCAACACCTGCa                                                        <  1:207818/148‑1 (MQ=255)                                                                                 gTTTTTACGCTTCATCGCTGGATGTAAATATGCTCTTGAAGCTGACATTCGGCGACCACGATGCACTTGACTATCGAACGGTCCAAAACCTTCTTGTCTAAATCCATTCACAGCAGGTGTTTTATGATAGCCTGCTTCAACACCTGCa                                                        <  2:350954/148‑1 (MQ=255)                                                                                     ttACGCTTCATCGCTGGATGTAAATATGCTCTTGAAGCTGACATTCGGCGACCACGATGCACTTGACTATCGAACGGTCCAAAACCTTCTTGTCTAAATCCATTCACATCAGGTGTTTTATGATAGCCTGCTTCAACACCTGCATCaa                                                    <  2:52209/148‑1 (MQ=255)                                                                                        cGCTTCATCGCTGGATGTAAATATGCTCTTGAAGCTGACATTCGGCGACCACGATGCACTTGACTATCGAACGGTCCAAAACCTTCTTGTCTAAATCCATTCACATCAGGTGTTTTATGATAGCCTGCTTCAACACCTGCATCAAAGa                                                 <  2:9186/148‑1 (MQ=255)                                                                                        cGCTTCATCGCTGGATGTAAATATGCTCTTGAAGCTGACATTCGGCGACCACGATGCACTTGACTATCGAACGGTCCAAAACCTTCTTGTCTAAATCCATTCACATCAGGTGGTTTATGATAGCCTGCTTCAACACCTGCATCAAAGa                                                 <  2:43097/148‑1 (MQ=255)                                                                                          tttCATCGCTGGATGTAAATATGCTCTTGAAGCTGACATTCGGCGACCACGATGCACTTGACTATCGAACGGTCCAAAGCCTTCTTGTCTAAATCCATTCACATCAGGTGTTTTATGATAGCCTGCTTCAACACCTGCATCAAAGAAt                                               <  2:391932/147‑1 (MQ=255)                                                                                          cTTCATCGCTGGATGTAAATATGCTCTTGAAGCTGACATTCGGCGACCACGATGCACTTGACTATCGAACGGTCCAAAACCTTCTTGTCTAAATCCATTCACATCAGGTGTTTTATGATAGCCTGCTTCAACACCTGCATCAAAGAAt                                               <  2:52699/148‑1 (MQ=255)                                                                                          cTTCATCGCTGGATGTAAATATGCTCTTGAAGCTGACATTCGGCGACCACGATGCACTTGACTATCGAACGGTCCAAAACCTTCTTGTCTAAATCCATTCACATCAGGTGTTTTATGATAGCCTGCTTCAACACCTGCATCAAAGAAt                                               <  2:239546/148‑1 (MQ=255)                                                                                            tCATCGCTGGATGTAAATATGCTCTTGAAGCTGACATTCGGCGACCACGATGCACTTGACTATCGAACGGTCCAAAACCTTCTTGTCTAAATCCATTCACATCAGGTGTTTTATGATAGCCTGCTTCAACACCTGCATCAAAGAATGa                                             <  1:155048/148‑1 (MQ=255)                                                                                                cGCTGGATGTAAATATGCTCTTGAAGCTGACATTCGGCGACCACGATGCACTTGACTATCGAACGGTCCAAAACCTTCTTGTCTAAATCCATTCACATCAGGTGTTTTATGATAGCCTGCTTCAACACCTGCATCAAAGAATGACTgg                                         <  1:266924/148‑1 (MQ=255)                                                                                                cGCTGGATGTAAATATGCTCTTGAAGCTGACATTCGGCGACCACGATGCACTTGACTATCGAACGGTCCAAAACCTTCTTGTCTAAATCCATTCACATCAGGTGTTTTATGATAGCCTGCTTCAACACCTGCATCAAAGAATGACTgg                                         <  1:350894/148‑1 (MQ=255)                                                                                                cGCTGGATGTAAATATGCTCTTGAAGCTGACATTCGGCGACCACGATGCACTTGACTATCGAACGGTCCAAAACCTTCTTGTCTAAATCCATTCACATCAGGTGTTTTATGATAGCCTGCTTCAACACCTGCATCAAAGAATGACTgg                                         <  1:342152/148‑1 (MQ=255)                                                                                                      aTGTAAATATGCTCTTGAAGCTGACATTCGGCGACCACGATGCACTTGACTATCGAACGGTCCAAAACCTTCTTGTCTAAATCCATTCACATCAGGTGTTTTATGATAGCCTGCTTCAACACCTGCATCAAAGAATGACTGGAATaaa                                   <  2:242874/148‑1 (MQ=255)                                                                                                      aTGTAAATATGCTCTTGAAGCTGACATTCGGCGACCACGATGCACTTGACTATCGAACGGTCCAAAACCTTCTTGTCTAAATCCATTCACATCAGGTGTTTTATGATAGCCTGCTTCAACACCTGCATCAAAGAATGACTGGAATaaa                                   <  2:203392/148‑1 (MQ=255)                                                                                                      aTGTAAATATGCTCTTGAAGCTGACATTCGGCGACCACGATGCACTTGACTATCGAACGGTCCAAAACCTTCTTGTCTAAATCCATTCACATCAGGTGTTTTATGATAGCCTGCTTCAACACCTGCATCAAAGAATGACTGGAATaaa                                   <  2:258650/148‑1 (MQ=255)                                                                                                       tGTAAATATGCTCTTGAAGCTGACATTCGGCGACCACGATGCACTTGACTATCGAACGGTCCAAAACCTTCTTGTCTAAATCCATTCACATCAGGTGTTTTATGATAGCCTGCTTCAACACCTGCATCAAAGAATGACTGGAATAAAg                                  <  1:304383/148‑1 (MQ=255)                                                                                                                 ctctTGAAGCTGACATTCGGCGACCACGATGCACTTGACTATCGAACGGTCCAAAACCTTCTTGTCTAAATCCATTCACATCAGGTGTTTTATGATAGCCTGCTTCAACACCTGCATCAAAGAATGACTGGAATAAAGGATTCGTTGc                        <  2:21095/148‑1 (MQ=255)                                                                                                                     tGAAGCTGACATTCGGCGACCACGATGCACTTGACTATCGAACGGTCCAAAACCTTCTTGTCTAAATCCATTCACATCAGGTGTTTTATGATAGCCTGCTTCAACACCTGCATCAAAGAATGACTGGAATAAAGGATTCGTTGCTGGc                    >  1:357064/1‑148 (MQ=255)                                                                                                                              cATTCGGCGACCACGATGCACTTGACTATCGAACGGTCCAAAACCTTCTTGTCTAAATCCATTCACATCAGGTGTTTTATGATAGCCTGCTTCAACACCTGCATCAAAGAATGACTGGAATAAAGGATTCGTTGCTGGCCCTCGtttt           <  1:220668/148‑1 (MQ=255)                                                                                                                              cATTCGGCGACCACGATGCACTTGACTATCGAACGGTCCAAAACCTTCTTGTCTAAATCCATTCACATCAGGTGTTTTATGATAGCCTGCTTCAACACCTGCATCAAAGAATGACTGGAATAAAGGATTCGTTGCTGGCCCTCGtttt           <  1:351563/148‑1 (MQ=255)                                                                                                                                       aCCACGATGCACTTGACTATCGAACGGTCCAAAACCTTCTTGTCTAAATCCATTCACATCAGGTGTTTTATGATAGCCTGCTTCAACACCTGCATCAAAGAATGACTGGAATAAAGGATTCGTTGCTGGCCCTCGTTTTAACTTAAtt  <  1:205894/148‑1 (MQ=255)                                                                                                                                                |                                                                                                                                           TTTTTCTTATACGTAACACCAGTTGCTCTTCTACCTTCATAATGAATTTCAGTTACAAAGGCACGTGTTTCAACGGTTAAGTTTTTACGCTTCATCGCTGGATGTAAATATGCTCTTGAAGCTGACATTCGGCGACCACGATGTACTTGACTATCGAACGGTCCAAAACCTTCTTGTCTAAATCCATTCACATCAGGTGTTTTATGATAGCCTGCTTCAACACCTGCATCAAAGAATGACTGGAATAAAGGATTCGTTGCTGGCCCTCGTTTTAACTTAATT  >  CP000730/2750419‑2750700

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑