Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	325,935	G→A	P438L (CCC→CTC)	betA ←	choline dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	325,935	0	G	A	100.0%	175.8 / NA	51	P438L (CCC→CTC)	betA	choline dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (27/24);  total (27/24)

CGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTACGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGGGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACTCCTGCCAGTCCTGCTCGTGCGACATGTAGTTAAACAGAATCGCCGGATGCTGGTGCGGGTCGCGGGATTTAATCCGCACATGCCCACGG  >  NC_000913/325798‑326069                                                                                                                                          |                                                                                                                                       cgtcgtAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTACGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGAGt                                                                                                                                      >  1:565006/1‑140 (MQ=255)       aaCCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTACGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGAGt                                                                                                                                      >  1:301051/1‑134 (MQ=255)       aaCCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTACGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGAGTTGATGc                                                                                                                                >  1:320483/1‑140 (MQ=255)             ttttGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTACGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGAGTTGATGCATGAtc                                                                                                                          >  2:177930/1‑140 (MQ=255)                tGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTACGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGAGTTGATGCATGATCTcg                                                                                                                       >  2:188704/1‑140 (MQ=255)                   aGGTCCCGCACGGATGGAAGGCGGTTTCGGCGTGGTTACGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGAGTTGATGCATGATCTCGCgg                                                                                                                    <  2:330789/140‑1 (MQ=255)                   aGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTACGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGAGTTGATGCATGATCTCGCgg                                                                                                                    <  1:188704/140‑1 (MQ=255)                   aGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTACGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGAGTTGATGCATGATCTCGCgg                                                                                                                    <  1:391391/140‑1 (MQ=255)                                          gTTTCGGCGTGGTTACGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGAGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGaa                                                                                             >  2:761634/1‑140 (MQ=255)                                          gTTTCGGCGTGGTTACGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGAGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGaa                                                                                             >  2:143106/1‑140 (MQ=255)                                          gTTTCGGCGTGGTTACGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGAGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGaa                                                                                             >  1:343558/1‑140 (MQ=255)                                            ttCGGCGTGGTTACGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGAGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACt                                                                                           <  2:664405/140‑1 (MQ=255)                                                 cGTGGTTACGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGAGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACTCCTGc                                                                                      <  2:124485/140‑1 (MQ=255)                                                     gTTACGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGAGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACTCCTGCCAGt                                                                                  >  2:143266/1‑140 (MQ=255)                                                         cGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGAGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACTCCTGCCAGTCCTg                                                                              <  2:534789/140‑1 (MQ=255)                                                            aCGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGAGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACTCCTGCCAGTCCTGCTc                                                                           <  2:320483/140‑1 (MQ=255)                                                               aaCTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGAGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACTCCTGCCAGTCCTGCTCGTg                                                                        >  1:387254/1‑140 (MQ=255)                                                               aaCTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGAGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACTCCTGCCAGTCCTGCTCGTg                                                                        <  2:697967/140‑1 (MQ=255)                                                                aCTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGAGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACTCCTGCCAGTCCTGCTCGTGc                                                                       >  2:380132/1‑140 (MQ=255)                                                                 cTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGAGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACTCCTGCCAGTCCTGCTCGTGCg                                                                      <  1:788236/140‑1 (MQ=255)                                                                      cGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGAGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACTCCTGCCAGTCCTGCTCGTGCGACATg                                                                 >  1:4179/1‑140 (MQ=255)                                                                      cGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGAGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACTCCTGCCAGTCCTGCTCGTGCGACATg                                                                 <  1:380132/140‑1 (MQ=255)                                                                         gCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGAGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACTCCTGCCAGTCCTGCTCGTGCGACATGTAg                                                              <  1:761634/140‑1 (MQ=255)                                                                           tGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGAGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACTCCTGCCAGTCCTGCTCGTGCGACATGTAGtt                                                            >  1:673452/1‑140 (MQ=255)                                                                            gTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGAGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACTCCTGCCAGTCCTGCTCGTGCGACATGTAGTTa                                                           >  1:827780/1‑140 (MQ=255)                                                                              tCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGAGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACTCCTGCCAGTCCTGCTCGTGCGACATGTAGTTaaa                                                         <  2:4179/140‑1 (MQ=255)                                                                                       tGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGAGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACTCCTGCCAGTCCTGCTCGTGCGACATGTAGTTAAACAGAATCGc                                                >  2:804986/1‑140 (MQ=255)                                                                                             tCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGAGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACTCCTGCCAGTCCTGCTCGTGCGACATGTAGTTAAACAGAATCGCCGGATg                                          <  1:177930/140‑1 (MQ=255)                                                                                                    cGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGAGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACTCCTGCCAGTCCTGCTCGTGCGACATGTAGTTAAACAGAATCGCCGGATGCTGGTGc                                   >  1:299874/1‑140 (MQ=255)                                                                                                      gggCTGATTTCGCGGCCACGATACTGATCCAGCGCGAGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACTCCTGCCAGTCCTGCTCGTGCGACATGTAGTTAAACAGAATCGCCGGa                                            >  1:675459/1‑129 (MQ=255)                                                                                                      gggCTGATTTCGCGGCCACGATACTGATCCAGCGCGAGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACTCCTGCCAGTCCTGCTCGTGCGACATGTAGTTAAACAGAATCGCCGGa                                            <  2:675459/129‑1 (MQ=255)                                                                                                      gggCTGATTTCGCGGCCACGATACTGATCCAGCGCGAGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACTCCTGCCAGTCCTGCTCGTGCGACATGTAGTTAAACAGAATCGCCGGATGCTGGTGCgg                                 >  2:327736/1‑140 (MQ=255)                                                                                                         cTGATTTCGCGGCCACGATACTGATCCAGCGCGAGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACTCCTGCCAGTCCTGCTCGTGCGACATGTAGTTAAACAGAATCGCCGGATGCTGGTGCGGGTc                              <  2:343558/140‑1 (MQ=255)                                                                                                          tGATTTCGCGGCCACGATACTGATCCAGCGCGAGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACTCCTGCCAGTCCTGCTCGTGCGACATGTAGTTAAACAGAATCGCCGGATGCTGGTGCGGGTcg                             <  2:69579/140‑1 (MQ=255)                                                                                                                     ccACGATACTGATCCAGCGCGAGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACTCCTGCCAGTCCTGCTCGTGCGACATGTAGTTAAACAGAATCGCCGGATGCTGGTGCGGGTCGCGGGATTTAAt                  <  1:143266/140‑1 (MQ=255)                                                                                                                      cACGATACTGATCCAGCGCGAGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACTCCTGCCAGTCCTGCTCGTGCGACATGTAGTTAAACAGAATCGCCGGATGCTGGTGCGGGTCGCGGGATTTAATc                 >  2:506711/1‑140 (MQ=255)                                                                                                                        cGATACTGATCCAGCGCGAGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACTCCTGCCAGTCCTGCTCGTGCGACATGTAGTTAAACAGAATCGCCGGATGCTGGTGCGGGTCGCGGGATTTAATCCg               <  2:565006/140‑1 (MQ=255)                                                                                                                          aTACTGATCCAGCGCGAGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACTCCTGCCAGTCCTGCTCGTGCGACATGTAGTTAAACAGAATCGCCGGATGCTGGTGCGGGTCGCGGGATTTAATCCGca             <  1:494174/140‑1 (MQ=255)                                                                                                                           tACTGATCCAGCGCGAGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACTCCTGCCAGTCCTGCTCGTGCGACATGTAGTTAAACAGAATCGCCGGATGCTGGTGCGGGTCGCGGGATTTAATCCGcac            <  1:454068/140‑1 (MQ=255)                                                                                                                           tACTGATCCAGCGCGAGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACTCCTGCCAGTCCTGCTCGTGCGACATGTAGTTAAACAGAATCGCCGGATGCTGGTGCGGGTCGCGGGATTTAATCCGcac            <  1:99121/140‑1 (MQ=255)                                                                                                                            aCTGATCCAGCGCGAGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACTCCTGCCAGTCCTGCTCGTGCGACATGTAGTTAAACAGAATCGCCGGATGCTGGTGCGGGTCGCGGGATTTAATCCGcaca           <  1:506711/140‑1 (MQ=255)                                                                                                                              tGATCCAGCGCGAGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACTCCTGCCAGTCCTGCTCGTGCGACATGTAGTTAAACAGAATCGCCGGATGCTGGTGCGGGTCGCGGGATTTAATCCGCACATg         <  2:673452/140‑1 (MQ=255)                                                                                                                                 tCCAGCGCGAGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACTCCTGCCAGTCCTGCTCGTGCGACATGTAGTTAAACAGAATCGCCGGATGCTGGTGCGGGTCGCGGGATTTAATCCGCACATGccc      >  2:537683/1‑140 (MQ=255)                                                                                                                                     gcgcGAGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACTCCTGCCAGTCCTGCTCGTGCGACATGTAGTTAAACAGAATCGCCGGATGCTGGTGCGGGTCGCGGGATTTAATCCGCACATGCCCACgg  >  2:412241/1‑140 (MQ=255)                                                                                                                                     gcgcGAGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACTCCTGCCAGTCCTGCTCGTGCGACATGTAGTTAAACAGAATCGCCGGATGCTGGTGCGGGTCGCGGGATTTAATCCGCACATGCCCACgg  >  1:424963/1‑140 (MQ=255)                                                                                                                                     gcgcGAGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACTCCTGCCAGTCCTGCTCGTGCGACATGTAGTTAAACAGAATCGCCGGATGCTGGTGCGGGTCGCGGGATTTAATCCGCACATGCCCACgg  >  2:334806/1‑140 (MQ=255)                                                                                                                                     gcgcGAGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACTCCTGCCAGTCCTGCTCGTGCGACATGTAGTTAAACAGAATCGCCGGATGCTGGTGCGGGTCGCGGGATTTAATCCGCACATGCCCACgg  >  1:80338/1‑140 (MQ=255)                                                                                                                                     gcgcGAGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACTCCTGCCAGTCCTGCTCGTGCGACATGTAGTTAAACAGAATCGCCGGATGCTGGTGCGGGTCGCGGGATTTAATCCGCACATGCCCACgg  >  1:91247/1‑140 (MQ=255)                                                                                                                                     gcgcGAGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACTCCTGCCAGTCCTGCTCGTGCGACATGTAGTTAAACAGAATCGCCGGATGCTGGTGCGGGTCGCGGGATTTAATCCGCACATGCCCACgg  >  2:175876/1‑140 (MQ=255)                                                                                                                                      cgcgAGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACTCCTGCCAGTCCTGCTCGTGCGACATGTAGTTAAACAGAATCGCCGGATGCTGGTGCGGGTCGCGGGATTTAATCCGcaca           <  2:361494/130‑1 (MQ=255)                                                                                                                                      cgcgAGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACTCCTGCCAGTCCTGCTCGTGCGACATGTAGTTAAACAGAATCGCCGGATGCTGGTGCGGGTCGCGGGATTTAATCCGcaca           >  1:361494/1‑130 (MQ=255)                                                                                                                                          |                                                                                                                                       CGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTACGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGGGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACTCCTGCCAGTCCTGCTCGTGCGACATGTAGTTAAACAGAATCGCCGGATGCTGGTGCGGGTCGCGGGATTTAATCCGCACATGCCCACGG  >  NC_000913/325798‑326069

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: ‑