Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	326,693	A→G	G185G (GGT→GGC)	betA ←	choline dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	326,693	0	A	G	100.0%	154.0 / NA	45	G185G (GGT→GGC)	betA	choline dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (20/25);  total (20/25)

TTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCAC  >  NC_000913/326558‑326815                                                                                                                                        |                                                                                                                           tttGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGGCCaa                                                                                                                        <  2:90923/140‑1 (MQ=255)     ccGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGGCCAAAAcc                                                                                                                    <  2:255038/140‑1 (MQ=255)       gTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGGCCAAAACCtt                                                                                                                  >  2:516587/1‑140 (MQ=255)                gatgTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGGCCAAAACCTTCCTGCTGAt                                                                                                         <  1:714845/140‑1 (MQ=255)                 atgTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGGCCAAAACCTTCCTGCTGATa                                                                                                        <  1:516587/140‑1 (MQ=255)                 atgTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGGCCAAAACCTTCCTGCTGATa                                                                                                        <  1:946853/140‑1 (MQ=255)                 atgTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGGCCAAAACCTTCCTGCTGATa                                                                                                        <  1:850331/140‑1 (MQ=255)                 atgTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGGCCAAAACCTTCCTGCTGATa                                                                                                        <  1:752460/140‑1 (MQ=255)                     gATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGGCCAAAACCTTCCTGCTGATAACCg                                                                                                    <  2:427353/140‑1 (MQ=255)                     gATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGGCCAAAACCTTCCTGCTGATAACCg                                                                                                    <  1:710232/140‑1 (MQ=255)                           tCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGGCCAAAACCTTCCTGCTGATAACCGTTgaga                                                                                              <  2:280549/140‑1 (MQ=255)                                 cGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGGCCAAAACCTTCCTGCTGATAACCGTTGAGATCGTcc                                                                                        >  2:727563/1‑140 (MQ=255)                                 cGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGAGCCAAAACCTTCCTGCTGATAACCGTTGAGATCGTcc                                                                                        >  1:693654/1‑140 (MQ=255)                                     aGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGGCCAAAACCTTCCTGCTGATAACCGTTgaga                                                                                              >  2:604142/1‑130 (MQ=255)                                     aGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGGCCAAAACCTTCCTGCTGATAACCGTTgaga                                                                                              <  1:604142/130‑1 (MQ=255)                                     aGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGGCCAAAACCTTCCTGCTGATAACCGTTGAGATc                                                                                            >  2:89204/1‑132 (MQ=255)                                     aGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGGCCAAAACCTTCCTGCTGATAACCGTTGAGATc                                                                                            <  1:89204/132‑1 (MQ=255)                                     aGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGGCCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTgc                                                                                    >  2:635157/1‑140 (MQ=255)                                     aGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGGCCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTgc                                                                                    >  1:191819/1‑140 (MQ=255)                                     aGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGGCCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTgc                                                                                    >  1:615909/1‑140 (MQ=255)                                     aGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGGCCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTg                                                                                     >  1:78174/1‑139 (MQ=255)                                     aGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGGCCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTg                                                                                     <  2:78174/139‑1 (MQ=255)                                         cGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGGCCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGcgcg                                                                                 <  1:652281/139‑1 (MQ=255)                                         cGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGGCCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGcgcg                                                                                 >  2:652281/1‑139 (MQ=255)                                               gTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGGCCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGcc                                                                          >  1:820672/1‑140 (MQ=255)                                               gTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGGCCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGcc                                                                          >  2:64909/1‑140 (MQ=255)                                               gTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGGCCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGcc                                                                          >  2:315528/1‑140 (MQ=255)                                                tCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGGCCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGccc                                                                         <  1:605162/140‑1 (MQ=255)                                                tCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCAGGCCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGccc                                                                         >  2:618598/1‑140 (MQ=255)                                                         gACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGGCCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACg                                                                <  1:574382/140‑1 (MQ=255)                                                              gATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGGCCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGc                                                           >  1:187995/1‑140 (MQ=255)                                                              gATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGGCCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGc                                                           >  1:203905/1‑140 (MQ=255)                                                              gATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGGCCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGc                                                           >  2:611039/1‑140 (MQ=255)                                                                     ccTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGGCCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAatc                                                    >  1:493210/1‑140 (MQ=255)                                                                      cTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGGCCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAatca                                                   <  1:611039/140‑1 (MQ=255)                                                                      cTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGGCCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAatca                                                   <  1:727563/140‑1 (MQ=255)                                                                          tCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGGCCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGc                                               <  2:615909/140‑1 (MQ=255)                                                                               aTAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGGCCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCaa                                          <  1:315528/140‑1 (MQ=255)                                                                                tAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGGCCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCaaa                                         <  1:618598/140‑1 (MQ=255)                                                                                   ccACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGGCCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAg                                      >  2:22300/1‑140 (MQ=255)                                                                                   ccACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGGCCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAg                                      <  1:64909/140‑1 (MQ=255)                                                                                        cgcgGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGGCCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGAt                                 <  1:21083/140‑1 (MQ=255)                                                                                              gCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGGCCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGc                           <  2:820672/140‑1 (MQ=255)                                                                                                                cggcgTGACGGTGCGATCCATCGGGCCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGa         >  1:713638/1‑140 (MQ=255)                                                                                                                       aCGGTGCGATCCATCGGGCCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCAc  <  1:635157/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                           TTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCAC  >  NC_000913/326558‑326815

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