Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,492,470	C→A	G19V (GGC→GTC)	idnT ←	L‑idonate/5‑ketogluconate/gluconate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,492,470	0	C	A	100.0%	56.6 / NA	21	G19V (GGC→GTC)	idnT	L‑idonate/5‑ketogluconate/gluconate transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (10/11);  total (10/11)

CCAAAAGTAGCAATCAGCGTAGTGGCGATACGTTGTGCCGCACCCGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCA‑TTGCCAGCCCGCCGAG‑CGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGGGG‑CTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGG  >  NC_000913/4492248‑4492708                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 ccAAAAGTAGCAATCAGCGTAGTGGCGATACGTTGTGCCGCACCCGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCA‑TTGCCAGCCCGCCGAG‑CGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGAc                                                                                                                                                                                                                      >  1:125398/1‑251 (MQ=255)     aaGTAGCAATCAGCGTAGTGGCGATACGTTGTGCCGCACCCGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCA‑TTGCCAGCCCGCCGAG‑CGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCt                                                                                                                                                                                                                  <  2:100614/251‑1 (MQ=255)      aGTAGCAATCAGCGTAGTGGCGATACGTTGTGCCGCACCCGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCA‑TTGCCAGCCCGCCGAG‑CGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTg                                                                                                                                                                                                                 <  1:233254/251‑1 (MQ=255)         agaaacaatcggaGTGGCGATACTGTGGGCCGCACCCGTATCAGAAATCACCTGGCCTAACATGGCCCCGAAACCGATGA‑CACGTGCCAGCCCGCCGAG‑CGTGATCCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCCTCCATCCCCTCGGCAAATCCGACGACGGCAGCTACCAGATCGAGGGCAATAAAGCCGTTAACTTTAAAGACGAACAGCAGGATAAGAAGCAGCGCGACGACTGcc                                                                                                                                                                                                               <  2:226903/238‑1 (MQ=255)                        ggctaTACGTTTTGCCGCACCCTTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCA‑TTGCCAGCCCGCCGAG‑CGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTCCCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAAt                                                                                                                                                                                                <  2:171645/246‑1 (MQ=255)                          gtaTACGTTTTGCCGCACCCGTATCAGAAATCAGCTTGCCTAACATGGCCCCTAAACCGAGGATCA‑TTGCCAGCCCTCCGAG‑CGTGCTGCCGATACCATTTTTTATAGAGTGCAGGACGGCCTGTGTATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGTAATAATGATTAATGGc                                                                                                                                                                                             <  1:166297/249‑1 (MQ=255)                          ggATACGTTGTGCCGCACCCGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCA‑TTGCCAGCCCGCCGAG‑CGTACTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACTGCAGCTACCAGAACGCGGGCAATAAAGCCGTTAACTTTAATGACGATCATAAGGATAAGAAGCAGCGCGACGCCTGCCGAAATAATGATTAATGGc                                                                                                                                                                                             >  2:969/2‑251 (MQ=255)                               cGTTGGGCCGCCCCCGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCA‑TTGCCAGCCCGCCGAG‑CGGGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGt                                                                                                                                                                                        <  2:125398/251‑1 (MQ=255)                                 ttGTGCCGCACCCGTATCAGAAATCAGCTTGCCTAACATCGCCCCGAACCCGAGGATCA‑TTCCCAACCCGCCGAG‑CGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCCCATCCGACGACGGCAGCTACCAGAACGAGGGCAATCAAGCCGTCAACTTTAAAGACGTTCTTCCGGATATGACGCAGCGCGACGCCTGCCGCCCTAATGATTAATGGCATAGTCt                                                                                                                                                                                      >  2:132961/1‑251 (MQ=255)                                  tgttccGCACCCGGATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCA‑TTGCCAGCCCGCCGAG‑CGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGGCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGAGAAGAAGCAGCGCGACGCCGGCCGCAATAATGATTAATGGCATAGTCt                                                                                                                                                                                      <  2:230778/246‑1 (MQ=255)                                                                                 acattggATCA‑TGGCTACCCCGCCGAGGCGT‑CTGCCGAACCCAGTTTGTAACGAGGGCAGGACGGCCTGTGCATCCATCCCTCCGGCAAATCAGACGACGGCAGCTACCAGACCGAGGGCATTAAAGCCGTTAACGTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGGGCCAGAAAAGCAGCGCTTTCGCGCTTCACGT‑Ag                                                                                                                                     <  2:151009/246‑1 (MQ=255)                                                                                                          gAG‑CGTGCTGCAGATAACATCTTGTATAGAGTGCAGGACGGCCTGTGCATACATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATAATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGGGGGCTGTTAAACAGCCACTa                                                                                                              >  2:28743/1‑250 (MQ=255)                                                                                                                      cGATAACATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCAt                                                                                                  >  2:247125/1‑251 (MQ=255)                                                                                                                                            cAGGACGGCCTGTGCATCCATCACTTCGGCAAATCCGACGCCGGCAGCTCCCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTc                                                                                                                                                   >  2:189898/1‑181 (MQ=255)                                                                                                                                            cAGGACGGCCTGCGCATCCATCCCTTCGGCAAATCCGACGACTGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGACGATGATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCTGAAAAGCAGCGCTTTc                                                                                                                                                   <  1:189898/181‑1 (MQ=255)                                                                                                                                                                                                   aCGAGGGCAATAAAGCCGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGc                     >  2:88983/1‑251 (MQ=255)                                                                                                                                                                                                     gaggaaaaTAAAACCGTAAATTTGAAAGACGATCATAAGGGTAAAAAGCAGCGGGACGCCTGCCGCAAAAATGATTAATGGCATAGTGTTTATCTCTTATTAGTGCAAGAAAAGCAGCGCTTACGCGCTGCACGG‑AGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAAGCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGGTCCTGCgg                   <  2:70635/245‑1 (MQ=255)                                                                                                                                                                                                         gCAATAAAGCCGTTAACTTTAAAGACGAGCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGAtct               <  1:168026/251‑1 (MQ=255)                                                                                                                                                                                                                  gcGTTAACTTTACAGAAGATTATCAGGATAAGAAGCAGCCCGACGCCTGCCGCAATAATGATTAATGGCATAGCCTTTATCTCTTATTAGTACCAGAAAAGCACCGCTTTCGCGCTGCACGTCGGGGGGG‑CTGTTAAACAGCCCATAACCTGCCGCCATCAACAACCAACAGGTGGCAGTTTAAGCAATCAGAGGCTTTTAACGAAAGGAACACAGCAGCAACAATCATTTcaagcggtgttccccagcg       >  2:95364/2‑232 (MQ=255)                                                                                                                                                                                                                   cGTTAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGGGGGCTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTg     >  2:100353/1‑251 (MQ=255)                                                                                                                                                                                                                      tAACTTTAAAGACGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGGGGGCTGTTAAACAGCCACTAACATGCCGCAATAAACAAACACCAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCAACAATCAGTTCCTGCGGATCTACCCAGCGTgcgg  >  1:120297/1‑251 (MQ=255)                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 CCAAAAGTAGCAATCAGCGTAGTGGCGATACGTTGTGCCGCACCCGTATCAGAAATCAGCTTGCCTAACATGGCCCCGAAACCGAGGATCA‑TTGCCAGCCCGCCGAG‑CGTGCTGCCGATACCATTTTGTATAGAGTGCAGGACGGCCTGTGCATCCATCCCTTCGGCAAATCCGACGACGGCAGCTACCAGAACGAGGGCAATAAAGCCGTTAACTTTAAAGCCGATCATCAGGATAAGAAGCAGCGCGACGCCTGCCGCAATAATGATTAATGGCATAGTCTTTATCTCTTATTAGTGCCAGAAAAGCAGCGCTTTCGCGCTGCACGT‑AGGGGGG‑CTGTTAAACAGCCACTAACATGCCGCCATCAACAAACAACAGGTGGCCGTTTACGAAATCAGAGGCTTTTGAAGAAAGGAACACCGCAGCACCAATCAGTTCCTGCGGATCTCCCCAGCGTGCGG  >  NC_000913/4492248‑4492708

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