Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	174,499	A→G	G128G (GGT→GGC)	hemL ←	glutamate‑1‑semialdehyde 2,1‑aminomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	174,499	0	A	G	100.0%	97.3 / NA	30	G128G (GGT→GGC)	hemL	glutamate‑1‑semialdehyde 2,1‑aminomutase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/14);  total (16/14)

ACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCAGGCGGATGGCGCTCATGGTCGCTTCAGTGCCGGAGTTCACCATGCGCACCATATCCATGGTCGGGACCAGTTCGGTCACCAGTTGCGCCATTTTCACTTCCATTTCGGTTGGTGCACCAAAGCTTAAACCACGCTCGGCGGCTTCAATCACGGCATTGCGGATTGCCGGATGGTTATGGCCCAGCACCATCGGCCCCCAGGAACCGACATAATCGATATAGGCTTTGCCATCAAC  >  NC_000913/174266‑174744                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                      aCCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAGCCACGGGCCAGGCGGAt                                                                                                                                                                                                                                      <  2:75340/251‑1 (MQ=255) aCCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGGATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAGCCACGGGCCAGGCGGAt                                                                                                                                                                                                                                      >  1:68850/1‑251 (MQ=255)          acgaTAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAGCCACGGGCCAGGCGGATGGCGCTCAt                                                                                                                                                                                                                             >  2:123828/1‑251 (MQ=255)                ataCAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAGCCACGGGCCAGGCGGATGGCGCTCATGGTCGc                                                                                                                                                                                                                       >  1:97693/1‑251 (MQ=255)                       cAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGCGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACACTTCAACTTTAATAATTTTGTCGCGACCGGTAAAGCCACGGGCCCGGCGGATGGCGCTCATGGGCCCTTCAGTg                                                                                                                                                                                                                >  2:77059/1‑251 (MQ=255)                       cAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAGCCACGGGCCAGGCGGATGGCGCTCATGGTCGCTTCAGTg                                                                                                                                                                                                                <  1:107000/251‑1 (MQ=255)                                 cGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATGTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAGCCACGGGCCAGGCGGATGGCGCTCATGTTCGCTTCAGTGCCGGAGTTCa                                                                                                                                                                                                      <  2:149685/251‑1 (MQ=255)                                         gCTCAAATCCGGCGCGGCCCGAAGCCAGATCATTATAAGCCCAGGTTAAGGCATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAGCCACGTGCCAGGCGGATGGCGCTCATGGTCGCTTCAGTGCCGGAGTTCACCATgcgc                                                                                                                                                                                              >  2:209866/1‑251 (MQ=255)                                                 gcggcgCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAGCCACGGGCCAGGCGGATGGCGCTCATGGTCGCTTCAGTGCCGGAGTTCACCATGCGCACCatat                                                                                                                                                                                       <  2:347100/250‑1 (MQ=255)                                                   ggcgCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAGCCACGGGCCAGGCGGATGGCGCTCATGGTCGCTTCAGTGCCGGAGTTCACCATGCGCACCATATcc                                                                                                                                                                                     >  2:166164/1‑250 (MQ=255)                                                    gcgcGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAGCCACGGGCCAGGCGGATGGCGCTCATGGTCGCTTCAGTGCCGGAGTTCACCATGCGCACCATATCCAt                                                                                                                                                                                   <  1:131390/251‑1 (MQ=255)                                                            aGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAGCCACGGGCCAGGCGGATGGCGCTCATGGTCGCTTCAGTGCCGGAGTTCACCATGCGCACCATATCCATGGTCggg                                                                                                                                                                            >  1:223847/1‑250 (MQ=255)                                                            aGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAGCCACGGGCCAGGCGGATGGCGCTCATGGTCGCTTCAGTGCCGGAGTTCACCATGCGCACCATATCCATGGTCGGGa                                                                                                                                                                           <  2:223847/251‑1 (MQ=255)                                                                   aGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAGCCACGGGCCAGGCGGATGGCGCTCATGGTCGCTTCAGTGCCGGAGTTCACCATGCGCACCATATCCATGGTCGGGACCAGTTc                                                                                                                                                                    <  1:197340/251‑1 (MQ=255)                                                                                     gttgaGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAGCCACGGGCCAGGCGGATGGCGCTCATGGTCGCTTCAGTGCCGGAGTTCACCATGCGCACCATATCCATGGTCGGGACCAGTTCGGTCACCAGTTGCGCCAt                                                                                                                                                  >  1:283574/5‑251 (MQ=255)                                                                                        aaGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAGCCACGGGCCAGGCGGATGGCGCTCATGGTCGCTTCAGTGCCGGAGTTCACCATGCGCACCATATCCATGGTCGGGACCAGTTCGGTCACCAGTTGCGCCAtttt                                                                                                                                               <  1:356/251‑1 (MQ=255)                                                                                                tttGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAGCCACGGGCCAGGCGGATGGCGCTCATGGTCGCTTCAGTGCCGGAGTTCACCATGCGCACCATATCCATGGTCGGGACCAGTTCGGTCACCAGTTGCGCCATTTTCACTTcc                                                                                                                                        >  1:325489/1‑250 (MQ=255)                                                                                                                                      ggcCTAACGTGAGGGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAGCCACGGGCCAGGCTGATGGCGCTCATGGTCGCTTCAGTGCCGGAGTTCACCATGCGCACCATATCCATGGTCGGGACCAGTTCGGTCACCATTTGCGCCATTTTCACTTCCATTTCGGTTGGTGCACCAAAGCTTAAACCACGCTcggcg                                                                                                 <  2:308065/251‑1 (MQ=255)                                                                                                                                        ccTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAGCCACGGGCCAGGCGGATGGCGCTCATGGTCGCTTCAGTGCCGGAGTTCACCATGCGCACCATATCCATGGTCGGGACCAGTTCGGTCACCAGTTGCGCCATTTTCACTTCCATTTCGGTTGGTGCACCAAAGCTTAAACCACGCTCggcggc                                                                                               <  2:97693/251‑1 (MQ=255)                                                                                                                                           aaCGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAGCCACGGGCCAGGCGGATGGCGCTCATGGTCGCTTCAGTGCCGGAGTTCACCATGCGCACCATATCCATGGTCGGGACCAGTTCGGTCACCAGTTGCGCCATTTTCACTGCCATTTCGGTTGGTGCACCAAAGCTTAAACCACGCTCGGCGGCtt                                                                                             <  1:163637/250‑1 (MQ=255)                                                                                                                                                          gAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAGCCACGGGCCAGGCGGATGGCGCTCATGGTCGCTTCAGTGCCGGAGTTCACCATGCGCACCATATCCATGGTCGGGACCAGTTCGGTCACCAGTTGCGCCATTTTCACTTCCATTTCGGTTGGTGCACCAAAGCTTAAACCACGCTCGGCGGCTTCAATCACGGCATTGCg                                                                             <  1:142073/251‑1 (MQ=255)                                                                                                                                                                                 gTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAGCCACGGGCCAGGCGGATGGCGCTCATGGTCGCTTCAGTGCCGGAGTTCACCATGCGCACCATATCCATGGTCGGGACCAGTTCGGTCACCAGTTGCGCCATTTTCACTTCCATTTCGGTTGGTGCACCAAAGCTTAAACCACGCTCGGCGGCTTCAATCACGGCATTGCGGATTGCCGGATGGTTATGGccc                                                       >  1:232467/1‑250 (MQ=255)                                                                                                                                                                                                            aaaTTTAATAATTTTGTCGCGACCGGTAAAGCCACGGGCCAGGCGGATGGCGCTCATGGTCGCTTCAGTGCCGGAGTTCACCATGCGCACCATATCCATGGTCGGGACCAGTTCGGTCACCAGTTGCGCCATTTTCACTTCCATTTCGGTTGGTGCACCAAAGCTTAAACCACGCTCGGCGGCTTCAATCACGGCATTGCGGATTGCCGGATGGTTATGGCCCAGCACCATCGGCCCCCAGGAACCGACAt                           <  1:166164/251‑1 (MQ=255)                                                                                                                                                                                                                 taataaTTTTGTCGCGACCGGTAAAGCCACGGGCCAGGCGGATGGCGCTCATGGTCGCTTCAGTGCCGGAGTTCACCATGCGCACCATATCCATGGTCGGGACCAGTTAGGTCTCCAGTTGCGCCATTTTCACTTCACTTTCGGTTGGTGCACCAAAGCTTAAACCACGCTCGGCGGCCTCAACCACGGCATTGCGGATTGCCGGCTGGTTATGGCCCAGCACCATCGGCCCCCAGGAACCGCCATAATCg                      >  1:356102/1‑251 (MQ=255)                                                                                                                                                                                                                   ataatTTTGTCGCGACCGGTAAAGCCACGGGCCAGGCGGATGGCGCTCATGGTCGCTTCAGTGCCGGAGTTCACCATGCGGACCATATCCATGGTCGGGACCAGTTCGGTCACCAGTTGCGCCATTTTCACTTCCATTTGGGTTGGTGCACCAAAGCTTAAACCACGCTCGGCGGCTTCAATCACGGCATTGCGGATTGCCGGATGGTTATGGCCCAGCACCATCGGCCCCCAGGAAACGACATAATCGat                    >  2:264647/1‑251 (MQ=255)                                                                                                                                                                                                                    taatTTTGTCGCGACCGGTAAAGCCACGGGCCAGGCGGATGGCGCTCATGGTCGCTTCAGTGCCGGAGTTCACCATGCGCACCATATCCATGGTCGGGACCAGTTCGGTCACCAGTTGCGCCATTTTCACTTCCATTTCGGTTGGTGCACCAAAGCTTAAACCACGCTCGGCGGCTTCAATCACGGCATTGCGGATTGCCGGATGGTTATGGc                                                         >  2:192900/1‑213 (MQ=255)                                                                                                                                                                                                                    taatTTTGTCGCGACCGGTAAAGCCACGGGCCAGGCGGATGGCGCTCATGGTCGCTTCAGTGCCGGAGTTCACCATGCGCACCATATCCATGGTCGGGACCAGTTCGGTCACCAGTTGCGCCATTTTCACTTCCATTTCGGTTGGTGCACCAAAGCTTAAACCACGCTCGGCGGCTTCAATCACGGCATTGCGGATTGCCGGATGGTTATGGc                                                         <  1:192900/213‑1 (MQ=255)                                                                                                                                                                                                                      atTTTGTCGCGACCGGTAAAGCCACGGGCCAGGCGGATGGCGCTCATGGTCGCTTCAGTGCCGGAGTTCACCATGAGCACCATATCCATGGTCGGGACCAGTTCGGTCACCAGTTGCGCCATTTTCACTTCCATTTCGGTTGGTGCACCAAAGCTTAAACCACGCTCGGCGGCTTCAATCCCGGCATTGCGGATTGCCGGATGGTTATGGCCCAGCACCATCGGCCCCCAGGAACCGACATAATCGatatt                 >  2:36610/1‑250 (MQ=255)                                                                                                                                                                                                                             cgcgACCGGTAAAGCCACGGGCCAGGCGGATGGCGCTCCTGGTCGCTTCAGTGCCGGAGTTCACCATGCGCACCATATCCATGGTCGGGACCAGTTCGGTCACCAGTTGCGCCCTTTTCACTTCCATTTCGGTTGGTGCACCAAAGCTTAAACCACGCTCGGCGGCTTCAATCACGGCATTGCGGATTGCCGGATGGTTATGGCCCAGCACCATCGGCCCCCCGGAACCGACATAATCGATATAGGCTTTg          >  1:174414/1‑251 (MQ=255)                                                                                                                                                                                                                                      tAAAGCCACGGGCCAGGCGGATGGCGCTCATGGTCGCTTCAGTGCCGGAGTTCACCATGCGCACCATATCCATGGTCGGGACCAGTTCGGTCACCAGTTGCGCCATTTTCACTTCCATTTCGGTTGGTGCACCAAAGCTTAAACCACGCTCGGCGGCTTCAATCACGGCATTGCGGATTGCCGGATGGTTATGGCCCAGCACCATCGGCCCCCAGGCACCGACATAATCGATATAGGCTTTGCCATCAAc  >  2:321758/1‑250 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                      ACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCAGGCGGATGGCGCTCATGGTCGCTTCAGTGCCGGAGTTCACCATGCGCACCATATCCATGGTCGGGACCAGTTCGGTCACCAGTTGCGCCATTTTCACTTCCATTTCGGTTGGTGCACCAAAGCTTAAACCACGCTCGGCGGCTTCAATCACGGCATTGCGGATTGCCGGATGGTTATGGCCCAGCACCATCGGCCCCCAGGAACCGACATAATCGATATAGGCTTTGCCATCAAC  >  NC_000913/174266‑174744

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