Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,264,691	Δ1 bp	coding (327/1332 nt)	tdcC ←	threonine/serine:H(+) symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,264,689	0	G	.	100.0%	187.8 / NA	44	coding (329/1332 nt)	tdcC	threonine/serine:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (28/16);  total (28/16)

AAGCCGAGCTGGTTTTCCCAGAACGTCATAAAGGTATTGGTAATAGTAACGCCATAAATCCACAGCAGTGGGCAAATCGCGAAGAAGTACAGGAACGTGATAACCACGCCGCCAGTTTTACCAAAATGCTCTTCCACCGTTTCCG  >  NC_000913/3264620‑3264764                                                                      |                                                                            aaGCCGAGCTGGTTTTCCCAGAACGTCATAAAGGTATTGGTAATAGTAACGCCATAAATCCACAGCAGT‑GGCAAATcgcg                                                                  >  1:548796/1‑80 (MQ=255) aaGCCGAGCTGGTTTTCCCAGAACGTCATAAAGGTATTGGTAATAGTAACGCCATAAATCCACAGCAGT‑GGCAAATcgcg                                                                  >  1:548798/1‑80 (MQ=255)      gAGCTGGTTTTCCCAGAACGTCATAAAGGTATTGGTAATAGTAACGCCATAAATCCACAGCAGT‑GGCAAATCGCGaagaag                                                            >  1:548781/1‑81 (MQ=255)       aGCTGGTTTTCCCAGAACGTCATAAAGGTATTGGTAATAGTAACGCCATAAATCCACAGCAGT‑GGCAAATCGCGaagaag                                                            >  2:548779/1‑80 (MQ=255)       aGCTGGTTTTCCCAGAACGTCATAAAGGTATTGGTAATAGTAACGCCATAAATCCACAGCAGT‑GGCAAATCGCGaagaag                                                            >  1:548784/1‑80 (MQ=255)           ggTTTTCCCAGAACGTCATAAAGGTATTGGTAATAGTAACGCCATAAATCCACAGCAGT‑GGCAAATCGCGAAGAAGTACAg                                                       >  1:548801/1‑81 (MQ=255)            gTTTTCCCAGAACGTCATAAAGGTATTGGTAATAGTAACGCCATAAATCCACAGCAGT‑GGCAAATCGCGAAGAAGTACAg                                                       >  1:548790/1‑80 (MQ=255)            gTTTTCCCAGAACGTCATAAAGGTATTGGTAATAGTAACGCCATAAATCCACAGCAGT‑GGCAAATCGCGAAGAAGTACAg                                                       <  2:548744/80‑1 (MQ=255)             ttttCCCAGAACGTCATAAAGGTATTGGTAATAGTAACGCCATAAATCCACAGCAGT‑GGCAAATCGCGAAGAAGTACAGGa                                                     <  1:548746/81‑1 (MQ=255)               ttCCCAGAACGTCATAAAGGTATTGGTAATAGTAACGCCATAAATCCACAGCAGT‑GGCAAATCGCGAAGAAGTACAGGaa                                                    <  2:548746/80‑1 (MQ=255)                  ccAGAACGTCATAAAGGTATTGGTAATAGTAACGCCATAAATCCACAGCAGT‑GGCAAATCGCGAAGAAGTACAGGAACGt                                                 >  2:548796/1‑80 (MQ=255)                   cAGAACGTCATAAAGGTATTGGTAATAGTAACGCCATAAATCCACAGCAGT‑GGCAAATCGCGAAGAAGTACAGGAACGTGa                                               >  2:548785/1‑81 (MQ=255)                    aGAACGTCATAAAGGTATTGGTAATAGTAACGCCATAAATCCACAGCAGT‑GGCAAATCGCGAAGAAGTACAGGAACGTGAt                                              <  2:548748/81‑1 (MQ=255)                    aGAACGTCATAAAGGTATTGGTAATAGTAACACCATAAATCCACAGCAGT‑GGCAAATCGCGAAGAAGTACAGGAACGTGAt                                              <  2:548747/81‑1 (MQ=255)                        cGTCATAAAGGTATTGGTAATAGTAACGCCATAAATCCACAGCAGT‑GGCAAATCGCGAAGAAGTACAGGAACGTGATAAcc                                          <  1:548750/81‑1 (MQ=255)                         gTCATAAAGGTATTGGTAATAGTAACGCCATAAATCCACAGCAGT‑GGCAAATCGCGAAGAAGTACAGGAACGTGATAACCa                                         >  2:548798/1‑81 (MQ=255)                          tCATAAAGGTATTGGTAATAGTAACGCCATAAATCCACAGCAGT‑GGCAAATCGCGAAGAAGTACAGGAACGTGATAACCAc                                        >  2:548802/1‑81 (MQ=255)                          tCATAAAGGTATTGGTAATAGTAACGCCATAAATCCACAGCAGT‑GGCAAATCGCGAAGAAGTACAGGAACGTGATAACCAc                                        >  2:548784/1‑81 (MQ=255)                             tAAAGGTATTGGTAATAGTAACGCCATAAATCCACAGCAGT‑GGCAAATCGCGAAGAAGTACAGGAACGTGATAACCAcgcc                                     <  2:548750/81‑1 (MQ=255)                              aaaGGTATTGGTAATAGTAACGCCATAAATCCACAGCAGT‑GGCAAATCGCGAAGAAGTACAGGAACGTGATAACCAcgccg                                    >  1:548792/1‑81 (MQ=255)                               aaGGTATTGGTAATAGTAACGCCATAAATCCACAGCAGT‑GGCAAATCGCGAAGAAGTACAGGAACGTGATAACCAcgccgc                                   <  1:548752/81‑1 (MQ=255)                                aGGTATTGGTAATAGTAACGCCATAAATCCACAGCAGT‑GGCAAATCGCGAAGAAGTACAGGAACGTGATAACCACgccgcc                                  <  2:548752/81‑1 (MQ=255)                                       ggTAATAGTAACGCCATAAATCCACAGCAGT‑GGCAAATCGCGAAGAAGTACAGGAACGTGATAACCACGCCGCCAGtttt                            <  2:548753/80‑1 (MQ=255)                                       ggTAATAGTAACGCCATAAATCCACAGCAGT‑GGCAAATCGCGAAGAAGTACAGGAACGTGATAACCACGCCGCCAGTTTTa                           >  1:548782/1‑81 (MQ=255)                                        gTAATAGTAACGCCATAAATCCACAGCAGT‑GGCAAATCGCGAAGAAGTACAGGAACGTGATAACCACGCCGCCAGTTTTAc                          >  1:548793/1‑81 (MQ=255)                                          aaTAGTAACGCCATAAATCCACAGCAGT‑GGCAAATCGCGAAGAAGTACAGGAACGTGATAACCACGCCGCCAGTTTTAcc                         >  1:548795/1‑80 (MQ=255)                                          aaTAGTAACGCCATAAATCCACAGCAGT‑GGCAAATCGCGAAGAAGTACAGGAACGTGATAACCACGCCGCCAGTTTTAcc                         >  1:548806/1‑80 (MQ=255)                                           aTAGTAACGCCATAAATCCACAGCAGT‑GGCAAATCGCGAAGAAGTACAGGAACGTGATAACCACGCCGCCAGTTTTACCaa                       <  2:548754/81‑1 (MQ=255)                                           aTAGTAACGCCATAAATCCACAGCAGT‑GGCAAATCGCGAAGAAGTACAGGAACGTGATAACCACGCCGCCAGTTTTACCa                        >  1:548813/1‑80 (MQ=255)                                              gTAACGCCATAAATCCACAGCAGT‑GGCAAATCGCGAAGAAGTACAGGAACGTGATAACCACGCCGCCAGTTTTACCAAAAt                    <  1:548756/81‑1 (MQ=255)                                              gTAACGCCATAAATCCACAGCAGT‑GGCAAATCGCGAAGAAGTACAGGAACGTGATAACCACGCCGCCAGTTTTACCAAAAt                    >  2:548788/1‑81 (MQ=255)                                               tAACGCCATAAATCCACAGCAGT‑GGCAAATCGCGAAGAAGTACAGGAACGTGATAACCACGCCGCCAGTTTTACCAAAATg                   <  2:548756/81‑1 (MQ=255)                                               tAACGCCATAAATCCACAGCAGT‑GGCAAATCGCGAAGAAGTACAGGAACGTGATAACCACGCCGCCAGTTTTACCAAAATg                   >  1:548808/1‑81 (MQ=255)                                                 aCGCCATAAATCCACAGCAGT‑GGCAAATCGCGAAGAAGTACAGGAACGTGATAACCACGCCGCCAGTTTTACCAAAATGc                  <  2:548757/80‑1 (MQ=255)                                                  cGCCATAAATCCACAGCAGT‑GGCAAATCGCGAAGAAGTACAGGAACGTGATAACCACGCCGCCAGTTTTACCAAAATGctc                <  2:548758/81‑1 (MQ=255)                                                  cGCCATAAATCCACAGCAGT‑GGCAAATCGCGAAGAAGTACAGGAACGTGATAACCACGCCGCCAGTTTTACCAAAATGctc                <  1:548760/81‑1 (MQ=255)                                                     cATAAATCCACAGCAGT‑GGCAAATCGCGAAGAAGTACAGGAACGTGATAACCACGCCGCCAGTTTTACCAAAATGCTCTTc             >  2:548808/1‑81 (MQ=255)                                                        aaaTCCACAGCAGT‑GGCAAATCGCGAAGAAGTACAGTAACGTGATAACCACTCCGCCAGTTTTACCAAAATGCTCTTCCAc          >  1:548805/1‑81 (MQ=255)                                                        aaaTCCACAGCAGT‑GGCAAATCGCGAAGAAGTACAGGAACGTGATAACCACGCCGCCAGTTTTACCAAAATGCTCTTCCAc          >  1:548802/1‑81 (MQ=255)                                                        aaaTCCACAGCAGT‑GGCAAATCGCGAAGAAGTACAGGAACGTGATAACCACGCCGCCAGTTTTACCAAAATGCTCTTCCAc          >  2:548814/1‑81 (MQ=255)                                                             cacaGCAGT‑GGCAAATCGCGAAGAAGTACAGGAACGTGATAACCACGCCGCCAGTTTTACCAAAATGCTCTTCCACCGttt     >  1:548794/1‑81 (MQ=255)                                                             cacaGCAGT‑GGCAAATCGCGAAGAAGTACAGGAACGTGATAACCACGCCGCCAGTTTTACCAAAATGCTCTTCCACCGttt     >  2:548816/1‑81 (MQ=255)                                                                agcagT‑GGCAAATCGCGAAGAAGTACAGGAACGTGATAACCACGCCGCCAGTTTTACCAAAATGCTCTTCCACCGTTTCCg  >  1:548823/1‑81 (MQ=255)                                                                 gcagT‑GGCAAATCGCGAAGAAGTACAGGAACGTGATAACCACGCCGCCAGTTTTACCAAAATGCTCTTCCACCGTTTCCg  >  2:548806/1‑80 (MQ=255)                                                                      |                                                                            AAGCCGAGCTGGTTTTCCCAGAACGTCATAAAGGTATTGGTAATAGTAACGCCATAAATCCACAGCAGTGGGCAAATCGCGAAGAAGTACAGGAACGTGATAACCACGCCGCCAGTTTTACCAAAATGCTCTTCCACCGTTTCCG  >  NC_000913/3264620‑3264764

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