Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_vch_tpiA	1,196,232	C→T	T370T (ACC→ACT)	icd →	isocitrate dehydrogenase; e14 prophage attachment site; tellurite reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_vch_tpiA	1,196,232	0	C	T	100.0%	68.1 / NA	27	T370T (ACC→ACT)	icd	isocitrate dehydrogenase; e14 prophage attachment site; tellurite reductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (11/16);  total (11/16)

GCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGA  >  NC_000913_3_vch_tpiA/1196089‑1196311                                                                                                                                                |                                                                                gCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGc                                                                            >  1:262977/1‑149 (MQ=255) gCCCCTGCTGCAAACATCGGGGGGGGGGGGGGCCGGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAGGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGc                                                                            <  2:833459/149‑1 (MQ=255)   cccTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCgg                                                                          <  1:302014/149‑1 (MQ=255)                      gACGAATGCGCCCTGTCTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAACCGCTGTCCAAGTTGTTAAAgg                                                       >  1:726504/1‑149 (MQ=255)                         gAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTAt                                                    <  1:394898/149‑1 (MQ=38)                              cgcCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAg                                               <  1:404455/149‑1 (MQ=255)                              cgcCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAg                                               <  1:62119/149‑1 (MQ=37)                               gcCCTGTTTGAAGTCACCCCCGGTACTGCGCCGAAATATGCCGGTCAGTACAAAGTAAATCCTGGCTTTATTGTTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAgg                                              <  2:603290/149‑1 (MQ=25)                                cccTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGgc                                             <  2:139862/149‑1 (MQ=255)                                   tGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTATTAAAGGTATGGAAGgagcc                                          >  2:225373/1‑148 (MQ=37)                                         aaaccacccacGGTTCTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCaat                                    <  2:455425/146‑2 (MQ=25)                                         aaGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTCAAAGGTATGGAAGGCGCAATCaat                                    <  1:722297/149‑2 (MQ=25)                                          aGCCAACAACGGTACTGCCGCGAAATATGCCGGTCAGGACTGAGTTAACCCTGGCTGTTTTGTTCTCTCCGCTGGGAGGATGCTGCGCCATAGGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCaatg                                   <  2:579990/149‑3 (MQ=17)                                            ccacccacGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATTCTAGATCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATGTGGCTTGGACTGAAGCGGCTGACCAGACTGTTAAAGGAATGGAAGACGCAATccatgcc                                 >  2:581586/1‑143 (MQ=255)                                               cccacGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCaag                              >  2:577166/1‑148 (MQ=21)                                                   cGGTACTGCCTCGAAATAGACCGGGCACGACACAGGAAGGCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCg                          <  1:528287/149‑1 (MQ=14)                                                           cgcCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTAt                  <  1:29099/149‑1 (MQ=17)                                                             ccGAAATATGCCGGTCAGGACAAAGTCACTACCGGGTCTATTGTTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGa                <  1:777964/149‑1 (MQ=11)                                                             ccGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGa                >  1:683496/1‑149 (MQ=17)                                                               aaaaTATGCCGTTCAGGACAAAGTAAATCCTGGCTCCATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGCATTTCGCAATCAATGCCAAGACCGTAACTTATGACt              <  2:389117/148‑1 (MQ=14)                                                                 aaTATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCACTCAATGCCAAGACCGTCACTTATGACTTc            >  1:134999/1‑149 (MQ=17)                                                                    atGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGAGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAACGCGCAATCAATGCCCAACCCGTAACTTATGAATTctaa         >  1:694276/1‑146 (MQ=17)                                                                      gCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTCACTTATGACTTCGAACg       >  2:477393/1‑149 (MQ=14)                                                                       ccGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGt      >  1:216945/1‑149 (MQ=14)                                                                           tCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTga  <  2:215023/149‑1 (MQ=14)                                                                                aCAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGc                                                                            >  2:761073/1‑70 (MQ=37)                                                                                aCAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGc                                                                            <  1:761073/70‑1 (MQ=37)                                                                                                                                                |                                                                                GCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGA  >  NC_000913_3_vch_tpiA/1196089‑1196311

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