Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	552,997	A→G	V35A (GTC→GCC)	purE ←	N(5)‑carboxyaminoimidazole ribonucleotide mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	552,997	0	A	G	100.0%	100.5 / NA	29	V35A (GTC→GCC)	purE	N(5)‑carboxyaminoimidazole ribonucleotide mutase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/14);  total (15/14)

ATCATGCCTGGCAGATGCGCTGCGCCGCCTGCGCCCGCAATAATCACCTGATAACCGTTCTCTTCGGCGCTTTCGGCGAAGCTGAACAGTTTATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCT  >  NC_000913/552862‑553131                                                                                                                                        |                                                                                                                                       aTCATGCCTGGCAGATGCGCTGCGCCGCCTGCGCCCGCAATAATCACCTGATAACCGTTCTCTTCGGCGCTTTCGGCGAAGCTGAACAGTTTATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGGCAtt                                                                                                                                    <  1:123245/140‑1 (MQ=255)           gCAGATGCGCTGCGCCGCCTGCGCCCGCAATAATCACCTGATAACCGTTCTCTTCGGCGCTTTCGGCGAAGCTGAACAGTTTATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGGCATTCAGGATTTCg                                                                                                                          >  1:611206/1‑140 (MQ=255)                tgcgctgcgcCGCCTGCGCCCGCAATAATCACCTGATAACCGTTCTCTTCGGCGCTTTCGGCGAAGCTGAACAGTTTATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGGCATTCAGGATTTCGAAGAt                                                                                                                     >  1:111860/1‑140 (MQ=255)                tgcgctgcgcCGCCTGCGCCCGCAATAATCACCTGATAACCGTTCTCTTCGGCGCTTTCGGCGAAGCTGAACAGTTTATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGGCATTCAGGATTTCGAAGAt                                                                                                                     >  2:505923/1‑140 (MQ=255)                tgcgctgcgcCGCCTGCGCCCGCAATAATCACCTGATAACCGTTCTCTTCGGCGCTTTCGGCGAAGCTGAACAGTTTATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGGCATTCAGGATTTCGAAGAt                                                                                                                     >  2:397362/1‑140 (MQ=255)                       cgccgcCTGCGCCCGCAATAATCACCTGATAACCGTTCTCTTCGGCGCTTTCGGCGAAGCTGAACAGTTTATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGGCATTCAGGATTTCGaa                                                                                                                        >  2:712080/1‑130 (MQ=255)                       cgccgcCTGCGCCCGCAATAATCACCTGATAACCGTTCTCTTCGGCGCTTTCGGCGAAGCTGAACAGTTTATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGGCATTCAGGATTTCGaa                                                                                                                        <  1:712080/130‑1 (MQ=255)                          cgccTGCGCCCGCAATAATCACCTGATAACCGTTCTCTTCGGCGCTTTCGGCGAAGCTGAACAGTTTATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGGCATTCAGGATTTCGAAGATTTCGgcggcg                                                                                                           <  1:397362/140‑1 (MQ=255)                               gcgcCCGCAATAATCACCTGATAACCGTTCTCTTCGGCGCTTTCGGCGAAGCTGAACAGTTTATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCTGGGCATTCAGGATTTCGAAGATTTCGGCGGCGAACTg                                                                                                      <  2:599290/140‑1 (MQ=255)                                      cAATAATCACCTGATAACCGTTCTCTTCGGCGCTTTCGGCGAAGCTGAACAGTTTATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGGCATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTa                                                                                               <  1:197703/140‑1 (MQ=255)                                                             tctTCGGCGCTTTCGGCGAAGCTGAACAGTTTATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGGCATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCAt                                                                        >  1:1080759/1‑140 (MQ=255)                                                               ttCGGCGCTTTCGGCGAAGCTGAACAGTTTATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGGCATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCa                                                                      >  1:154923/1‑140 (MQ=255)                                                                              gAAGCTGAACAGTTTATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGGCATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACAcgcg                                                       <  2:1080759/140‑1 (MQ=255)                                                                                                            aGCAGAAACCACTTCAACGTGGTGCGGGGCATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATg                                <  1:597086/133‑1 (MQ=255)                                                                                                            aGCAGAAACCACTTCAACGTGGTGCGGGGCATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATg                                >  2:597086/1‑133 (MQ=255)                                                                                                            aGCAGAAACCACTTCAACGTGGTGCGGGGCATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTa                         >  1:926708/1‑140 (MQ=255)                                                                                                            aGCAGAAACCACTTCAACGTGGTGCGGGGCATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTa                         >  1:101230/1‑140 (MQ=255)                                                                                                                    ccACTTCAACGTGGTGCGGGGCATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCt                 >  1:34227/1‑140 (MQ=255)                                                                                                                     cACTTCAACGTGGTGCGGGGCATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCtg                <  2:926708/140‑1 (MQ=255)                                                                                                                     cACTTCAACGTGGTGCGGGGCATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCtg                <  2:1072306/140‑1 (MQ=255)                                                                                                                        ttCAACGTGGTGCGGGGCATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTaaaa                      >  1:144262/1‑131 (MQ=255)                                                                                                                        ttCAACGTGGTGCGGGGCATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTaaaa                      <  2:144262/131‑1 (MQ=255)                                                                                                                            aCGTGGTGCGGGGCATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGGCCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGca         <  2:111860/140‑1 (MQ=255)                                                                                                                            aCGTGGTGCGGGGCATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGca         <  1:556538/140‑1 (MQ=255)                                                                                                                             cGTGGTGCGGGGCATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGcac        >  1:945971/1‑140 (MQ=255)                                                                                                                                  tgCGGGGCATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAActc   >  1:326224/1‑140 (MQ=255)                                                                                                                                  tgCGGGGCATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGCGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAActt   >  1:644756/1‑139 (MQ=255)                                                                                                                                   gCGGGGCATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAActct  <  2:101230/140‑1 (MQ=255)                                                                                                                                   gCGGGGCATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAActct  <  2:945971/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                       ATCATGCCTGGCAGATGCGCTGCGCCGCCTGCGCCCGCAATAATCACCTGATAACCGTTCTCTTCGGCGCTTTCGGCGAAGCTGAACAGTTTATCGGGGGTGCGGTGAGCAGAAACCACTTCAACGTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCT  >  NC_000913/552862‑553131

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