Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	790,832	Δ1 bp	coding (198/1047 nt)	galT ←	galactose‑1‑phosphate uridylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	790,832	0	G	.	100.0%	65.1 / NA	15	coding (198/1047 nt)	galT	galactose‑1‑phosphate uridylyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (9/6);  total (9/6)

CTCCTGCCAGGTTTTGACGATTTCCGTCAATGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGG‑TTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGT  >  NC_000913/790622‑791057                                                                                                                                                                                                                   |                                                                                                                                                                                                                                   cTCCTGCCAGGTTTTGACGATTTCCGTCAATGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCCCATTCCCTGCGCAGAGGa                                                                                                                                                                                            >  2:194895/1‑250 (MQ=255)              ttGACGATTTCCGTCAATGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATc                                                                                                                                                                               <  1:38399/250‑1 (MQ=255)                  cGATTTCCGTCAATGCTGCAACTCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTgcg                                                                                                                                                                          <  2:35130/251‑1 (MQ=255)                                                                gtgATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAggg                                                                                                                            >  1:93156/1‑251 (MQ=255)                                                                       ggTGAAATGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGAAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGTATAGTGCGCAGGTAACACCCGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCCGGTGCGCtt                                                                                                                      >  2:134399/1‑250 (MQ=255)                                                                        gTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAg                                                                                                                    <  1:63811/251‑1 (MQ=255)                                                                                cAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCgcggtgc                                                                                                            >  1:84134/1‑251 (MQ=255)                                                                                           gCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAAcc                                                                                                   >  2:86041/1‑249 (MQ=255)                                                                                                            tctGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACGCCTGTTTGGGTGGCGTGTCCTGCGcccc                                                                                                                                     <  2:139022/198‑1 (MQ=255)                                                                                                            tctGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGcccc                                                                                                                                     >  1:139022/1‑198 (MQ=255)                                                                                                                              gATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTTATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGG‑TTGTAGCGGCGAt                                                             <  1:37369/251‑1 (MQ=255)                                                                                                                                        ttCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCATAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGTG‑TTTTTATCGCCTGTCACCCGCACATTTCCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACCCCTGGTTGGCTGGCGTTTCCTGCGCCGCCTGTCATGGGCGGTTGGCGCGGTGCGGGGAAACCAGAATCCATTTCGCGGTGAGCGGGGTTG‑ATCGGCGATGTGgatgat                                                    <  2:111105/250‑1 (MQ=255)                                                                                                                                                                     aaGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGG‑TTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTc                      >  2:16800/1‑251 (MQ=255)                                                                                                                                                                                         aaaaCGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGCGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGG‑TTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCtgt  >  2:92784/1‑251 (MQ=255)                                                                                                                                                                                         aaaaCGTAAGTCACGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTCCCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGG‑GTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGGCGTTCCTTAATCGGGATATTCCtgt  >  2:975/1‑251 (MQ=255)                                                                                                                                                                                                                   |                                                                                                                                                                                                                                   CTCCTGCCAGGTTTTGACGATTTCCGTCAATGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGG‑TTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGT  >  NC_000913/790622‑791057

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