Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_vch_tpiA	1,197,502	C→T	G11R (GGA→AGA)	ymfD ←	e14 prophage; putative SAM‑dependent methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_vch_tpiA	1,197,502	0	C	T	100.0%	95.2 / NA	32	G11R (GGA→AGA)	ymfD	e14 prophage; putative SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/13);  total (19/13)

TTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACT  >  NC_000913_3_vch_tpiA/1197365‑1197644                                                                                                                                          |                                                                                                                                               ttCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATACTTCTAATATACTCACAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCCtattat                                                                                                                                     >  1:90771/1‑149 (MQ=255)     tcaGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGtt                                                                                                                                 >  1:287001/1‑149 (MQ=255)         aaTATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCaaa                                                                                                                             <  1:18927/149‑1 (MQ=255)             tatgacttttCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTTTTTGTTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCaa                                                                                                                         <  1:257434/143‑1 (MQ=255)                   ttttCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTTTCAATGCTTCTAATATACTCGCAAAGAAATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTTTGAGCGAATATTAATTCTGTGCATATTCTATTTCAAAGAAAGttacc                                                                                                                   >  2:74704/1‑145 (MQ=255)                    ttttggCAGCCAATATTAAGGGCTGGTCCATTTTGCTCGATGCTTCTCATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTGGTTGCATTCTCTGAGCGAATATTAACTCTGCGCATATTATAGTTCAAAGCAAGTACCAt                                                                                                                  <  2:185088/144‑1 (MQ=255)                           aGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTa                                                                                                           <  1:31306/149‑1 (MQ=255)                                  aTCAAGGGCGTGGCCATTTTTCTTAATGCTTTTAATCTACTAGCAAAGATATATTTAGGGCATAATATTAGGATTTGCTGCATTATCTGAGCGAATATTAAATGTGTGCATATTATAGTTTATAGCAAGTATAATTCAATTATTTTTAt                                                                                                    >  2:239587/1‑149 (MQ=255)                                     aaGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTAtttt                                                                                                 <  2:287001/149‑1 (MQ=255)                                      aGGACGTAGCCATTTTTCTCAATGCTTCTAATATACTCACAAAAATATCTAAAGGGCATCGTATGAGGGTTTGCTGCATTCTCTGAGCGAATATTAACTTTGTGCATATTATAGTTCAAACAAAATACAATTAAATTAATATTAttttt                                                                                                >  2:149698/1‑149 (MQ=255)                                                        tcAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAg                                                                                              <  1:175664/133‑1 (MQ=255)                                                        tcAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAg                                                                                              >  2:175664/1‑133 (MQ=255)                                                        tcAATGCTTCTAATATACTCGCAAAGATATCTAGAGAGCATCGTATAAGGTTTTACTGCATTCTCTGAGCGCATATTAACTCTATGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCttttat                                                                              >  1:100293/1‑149 (MQ=255)                                                          aaTGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGATTAATAACTGCATTATATGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCttttattt                                                                            >  1:263874/1‑149 (MQ=255)                                                                tCTAATATACTCGCAAAGATATCTAGAGGGCCTCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGt                                                                      >  1:529061/1‑149 (MQ=255)                                                                    atatACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCgg                                                                  >  2:299231/1‑149 (MQ=255)                                                                       taCTCGCAAAGGTAACGAGAGGGCATCGTCTGAGGTTTTGCTGCACTCTGTGAGCGAATGTTGACTCTGTGCATATTATAGTGCAGAGCAAGTACAATTACATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATTGTTCggggt                                                               >  2:522344/1‑147 (MQ=255)                                                                       taCTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATAGTCCGGGAt                                                               >  2:371942/1‑149 (MQ=255)                                                                          tCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGAtttt                                                            >  2:242412/1‑149 (MQ=255)                                                                            gCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTCGTGAAAAATTCTTTTATTTTATCGTCTGGGATTTTTc                                                          >  2:573459/1‑149 (MQ=255)                                                                            gCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTc                                                          >  1:197960/1‑149 (MQ=255)                                                                                aGATGTTTATAGGGCATCGTATGAGGTTTTGCCGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTAACATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCgac                                                      <  1:161949/149‑1 (MQ=255)                                                                                         gagGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTc                                             <  1:503801/149‑1 (MQ=255)                                                                                               aTCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCt                                       >  1:584487/1‑149 (MQ=255)                                                                                                      gAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTATGTGCATATAATAGTTCAAAGCAAGTACCATTCAATTGATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCAGACTTTACTTTCAGTTCTGGGGTTa                                >  2:149871/1‑149 (MQ=255)                                                                                                        ggTTGTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAAc                              <  1:440927/149‑1 (MQ=255)                                                                                                         gTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACt                             <  1:74704/149‑1 (MQ=255)                                                                                                             tgttgcATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAACATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTca                         <  2:252015/146‑1 (MQ=255)                                                                                                               ctgcATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTcaaa                       <  1:299231/149‑1 (MQ=255)                                                                                                                 gcATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATAGCTGATGTTTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTcaaagc                     >  1:415086/1‑149 (MQ=255)                                                                                                                      ctctGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGa                <  2:283926/149‑1 (MQ=255)                                                                                                                                    tAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCTGACTTGACTTTCAGTTCTGGTCTTAACTGTCAAAGCAAAGATATTACcactcact  >  2:447571/1‑149 (MQ=255)                                                                                                                                       ctctGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTaa                               >  2:19698/1‑117 (MQ=255)                                                                                                                                       ctctGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTaa                               <  1:19698/117‑1 (MQ=255)                                                                                                                                          |                                                                                                                                               TTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACT  >  NC_000913_3_vch_tpiA/1197365‑1197644

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