Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,623,014	T→C	intergenic (+142/‑30)	purN → / → ppk	phosphoribosylglycinamide formyltransferase 1/polyphosphate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,623,014	0	T	C	100.0%	125.4 / NA	35	intergenic (+142/‑30)	purN/ppk	phosphoribosylglycinamide formyltransferase 1/polyphosphate kinase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (24/11);  total (24/11)

CGTAATTAAAGCGCCAGCTCTGCCGCTGGCGTTTTTCAATTCACCTGTAAATCGCAAGCTCCAGCAGTTTTTTTCCCCCTTTTCTGGCATAGTTGGACATCTGCCAATATTGCTCGCCATAATATCCAGGCAGTGTCCCGTGAATAAAACGGAGTAAAAGTGGTAATGGGTCAGGAAAAGCTATACATCGAAAAAGAGCTCAGTTGGTTATCGTTCAATGAACGCGTGCTTCAGGAAGCGGCGGACAAATCTAACCCGCTGATTGAA  >  NC_000913/2622879‑2623145                                                                                                                                        |                                                                                                                                    cGTAATTAAAGCGCCAGCTCTGCCGCTGGCGTTTTTCAATTCACCTGTAAATCGCAAGCTCCAGCAGTTTTTTTCCCCCTTTTCTGGCATAGTTGGACATCTGCCAATATTGCTCGCCATAATATCCAGGCAGTGCCCCg                                                                                                                                 >  1:211432/1‑140 (MQ=255)          aGCGCCAGCTCTGCCGCTGGCGTTTTTCAATTCACCTGTAAATCGCAAGCTCCAGCAGTTTTTTTCCCCCTTTTCTGGCATAGTTGGACATCTGCCAATATTGCTCGCCATAATATCCAGGCAGTGCCCCGTGAATaaaa                                                                                                                        >  1:925356/1‑140 (MQ=255)             gcCAGCTCTGCCGCTGGCGTTTTTCAATTCACCTGTAAATCGCAAGCTCCAGCAGTTTTTTTCCCCCTTTTCTGGCATAGTTGGACATCTGCCAATATTGCTCGCCATAATATCCAGGCAGTGCCCCGTGAATAAAACgg                                                                                                                     <  1:917028/140‑1 (MQ=255)                            ggCGTTTTTCAATTCACCTGTAAATCGCAAGCTCCAGCAGTTTTTTTCCCCCTTTTCTGGCATAGTTGGACATCTGCCAATATTGCTCGCCATAATATCCAGGCAGTGCCCCGTGAATAAAACGGAGTAAAAGTGGTAAt                                                                                                      >  1:1009917/1‑140 (MQ=255)                            ggCGTTTTTCAATTCACCTGTAAATCGCAAGCTCCAGCAGTTTTTTTCCCCCTTTTCTGGCATAGTTGGACATCTGCCAATATTGCTCGCCATAATATCCAGGCAGTGCCCCGTGAATAAAACGGAGTAAAAGTGGTAAt                                                                                                      <  1:781827/140‑1 (MQ=255)                             gCGTTTTTCAATTCACCTGTAAATCGCAAGCTCCAGCAGTTTTTTTCCCCCTTTTCTGGCATAGTTGGACATCTGCCAATATTGCTCGCCATAATATCCAGGCAGTGCCCCGTGAATAAAACGGAGTAAAAGTGGTa                                                                                                        <  1:522442/137‑1 (MQ=255)                             gCGTTTTTCAATTCACCTGTAAATCGCAAGCTCCAGCAGTTTTTTTCCCCCTTTTCTGGCATAGTTGGACATCTGCCAATATTGCTCGCCATAATATCCAGGCAGTGCCCCGTGAATAAAACGGAGTAAAAGTGGTa                                                                                                        >  2:522442/1‑137 (MQ=255)                                    tCAATTCACCTGTAAATCGCAAGCTCCAGCAGTTTTTTTCCCCCTTTTCTGGCATAGTTGGACATCTGCCAATATTGCTCGCCATAATATCCAGGCAGTGCCCCGTGAATAAAACGGAGTAAAAGTGGTAATGGGTCAgg                                                                                              <  1:46/140‑1 (MQ=255)                                            ccTGTAAATCGCAAGCTCCAGCAGTTTTTTTCCCCCTTTTCTGGCATAGTTGGACATCTGCCAATATTGCTCGCCATAATATCCAGGCAGTGCCCCGTGAATAAAACGGAGTAAAAGTGGTAATGGGTCAGGAAAAGCta                                                                                      <  2:1009917/140‑1 (MQ=255)                                            ccTGTAAATCGCAAGCTCCAGCAGTTTTTTTCCCCCTTTTCTGGCATAGTTGGACATCTGCCAATATTGCTCGCCATAATATCCAGGCAGTGCCCCGTGAATAAAACGGAGTAAAAGTGGTAATGGGTCAGGAAAAGCta                                                                                      <  1:945932/140‑1 (MQ=255)                                                     cGCAAGCTCCAGCAGTTTTTTTCCCCCTTTTCTGGCATAGTTGGACATCTGCCAATATTGCTCGCCATAATATCCAGGCAGTGCCCCGTGAATAAAACGGAGTAAAAGTGGTAATGGGTCAGGAAAAGCTATACATCGaa                                                                             >  1:696428/1‑140 (MQ=255)                                                      gCAAGCTCCAGCAGTTTTTTTCCCCCTTTTCTGGCATAGTTGGACATCTGCCAATATTGCTCGCCATAATAACCAGGCAGTGCCCCGTGAATAAAACGGAGTAAAAGTGGTAATGGGTCAGGAAAAGCTATACATCGaaa                                                                            <  1:167108/140‑1 (MQ=255)                                                       cAAGCTCCAGCAGTTTTTTTCCCCCTTTTCTGGCATAGTTGGACATCTGCCAATATTGCTCGCCATAATATCCAGGCAGTGCCTCGCGAATAAAACGGAGAAAAAGTGGTAATGGGTCAGGTAAAGCTATACATCGcaaa                                                                           >  2:167108/1‑140 (MQ=255)                                                       cAAGCTCCAGCAGTTTTTTTCCCCCTTTTCTGGCATAGTTGGACATCTGCCAATATTGCTCGCCATAATATCCAGGCAGTGCCCCGTGAATAAAACGGAGTAAAAGTGGTAATGGGTCAGGAAAAGCTATACATCGaaaa                                                                           >  1:936732/1‑140 (MQ=255)                                                          gCTCCAGCAGTTTTTTTCCCCCTTTTCTGGCATAGTTGGACATCTGCCAATATTGCTCGCCATAATATCCAGGCAGTGCCCCGTGAATAAAACGGAGTAAAAGTGGTAATGGGTCAGGAAAAGCTATACATCGAAAAaga                                                                        >  2:474777/1‑140 (MQ=255)                                                                                       gCATAGTTGGACATCTGCCAATATTGCTCGCCATAATATCCAGGCAGTGCCCCGTGAATAAAACGGAGTAAAAGTGGTAATGGGTCAGGAAAAGCTATACATCGAAAAAGAGCTCAGTTGGTTATCGTTCAATGAAcgcg                                           <  2:925356/140‑1 (MQ=255)                                                                                        cATAGTTGGACATCTGCCAATATTGCTCGCCATAATATCCAGGCAGTGCCCCGTGAATAAAACGGAGTAAAAGTGGTAATGGGTCAGGAAAAGCTATACATCGAAAAAGAGCTCAGTTGGTTATCGTTCAATGAACGCGt                                          >  1:412873/1‑140 (MQ=255)                                                                                          tAGTTGGACATCTGCCAATATTGCTCGCCATAATATCCAGGCAGTGCCCCGTGAATAAAACGGAGTAAAAGTGGTAATGGGTCAGGAAAAGCTATACATCGAAAAAGAGCTCAGTTGGTTATCGTTCAATGAACGCGTGc                                        >  2:1051385/1‑140 (MQ=255)                                                                                                 aCATCTGCCAATATTGCTCGCCATAATATCCAGGCAGTGCCCCGTGAATAAAACGGAGTAAAAGTGGTAATGGGTCAGGAAAAGCTATACATCGAAAAAGAGCTCAGTTGGTTATCGTTCAATGAACGCGTGCTTCAGGa                                 >  1:1047962/1‑140 (MQ=255)                                                                                                     cTGCCAATATTGCTCGCCATAATATCCAGGCAGTGCCCCGTGAATAAAACGGAGTAAAAGTGGTAATGGGTCAGGAAAAGCTATACATCGAAAAAGAGCTCAGTTGGTTATCGTTCAATGAACGCGTGCTTCAGGAAgcg                             >  2:864528/1‑140 (MQ=255)                                                                                                      tGCCAATATTGCTCGCCATAATATCCAGGCAGTGCCCCGTGAATAAAACGGAGTAAAAGTGGTAATGGGTCAGGAAAAGCTATACATCGAAAAAGAGCTCAGTTGGTTATCGTTCAATGAACGCGTGCTTCAGGAAgcgg                            >  2:315064/1‑140 (MQ=255)                                                                                                      tGCCAATATTGCTCGCCATAATATCCAGGCAGTGCCCCGTGAATAAAACGGAGTAAAAGTGGTAATGGGTCAGGAAAAGCTATACATCGAAAAAGAGCTCAGTTGGTTATCGTTCAATGAACGCGTGCTTCAGGAAgcgg                            >  2:638265/1‑140 (MQ=255)                                                                                                      tGCCAATATTGCTCGCCATAATATCCAGGCAGTGCCCCGTGAATAAAACGGAGTAAAAGTGGTAATGGGTCAGGAAAAGCTATACATCGAAAAAGAGCTCAGTTGGTTATCGTTCAATGAACGCGTGCTTCAGGAAgcgg                            >  1:9947/1‑140 (MQ=255)                                                                                                       gCCAATATTGCTCGCCATAATATCCAGGCAGTGCCCCGTGAATAAAACGGAGTAAAAGTGGTAATGGGTCAGGAAAAGCTATACATCGAAAAAGAGCTCAGTTGGTTATCGTTCAATGAACGCGTGCTTCAGGAAgcggc                           >  2:1019254/1‑140 (MQ=255)                                                                                                          aaTATTGCTCGCCATAATATCCAGGCAGTGCCCCGTGAATAAAACGGAGTAAAAGTGGTAATGGGTCAGGAAAAGCTATACATCGAAAAAGAGCTCAGTTGGTTATCGTTCAATGAACGCGTGCTTCAGGAAGCGGCGGa                        <  1:1019254/140‑1 (MQ=255)                                                                                                               tGCTCGCCATAATATCCAGGCAGTGCCCCGTGAATAAAACGGAGTAAAAGTGGTAATGGGTCAGGAAAAGCTATACATCGAAAAAGAGCTCAGTTGGTTATCGTTCAATGAACGCGTGCTTCAGGAAgcggcg                          >  1:1106367/1‑133 (MQ=255)                                                                                                               tGCTCGCCATAATATCCAGGCAGTGCCCCGTGAATAAAACGGAGTAAAAGTGGTAATGGGTCAGGAAAAGCTATACATCGAAAAAGAGCTCAGTTGGTTATCGTTCAATGAACGCGTGCTTCAGGAAgcgccg                          <  2:1106367/133‑1 (MQ=255)                                                                                                               tGCTCGCCATAATATCCAGGCAGTGCCCCGTGAATAAAACGGAGTAAAAGTGGTAATGGGTCAGGAAAAGCTATACATCGAAAAAGAGCTCAGTTGGTTATCGTTCAATGAACGCGTGCTTCAGGAAGCGGCGGACAAAt                   >  2:956823/1‑140 (MQ=255)                                                                                                                gCTCGCCATAATATCCAGGCAGTGCCCCGTGAATAAAACGGAGTAAAAGTGGTAATGGGTCAGGAAAAGCTATACATCGAAAAAGAGCTCAGTTGGTTATCGTTCAATGAACGCGTGCTTCAGGAAGCGGCGGACAAATc                  >  1:653854/1‑140 (MQ=255)                                                                                                                gCTCGCCATAATATCCAGGCAGTGCCCCGTGAATAAAACGGAGTAAAAGTGGTAATGGGTCAGGAAAAGCTATACATCGAAAAAGAGCTCAGTTGGTTATCGTTCAATGAACGCGTGCTTCAGGAAGCGGCGGACAAATc                  >  2:702511/1‑140 (MQ=255)                                                                                                                  tCGCCATAATATCCAGGCAGTGCCCCGTGAATAAAACGGAGTAAAAGTGGTAATGGGTCAGGAAAAGCTATACATCGAAAAAGAGCTCAGTTGGTTATCGTTCAATGAACGCGTGCTTCAGGAAGCGGCGGACAAATCTa                >  2:671238/1‑140 (MQ=255)                                                                                                                  tCGCCATAATATCCAGGCAGTGCCCCGTGAATAAAACGGAGTAAAAGTGGTAATGGGTCAGGAAAAGCTATACATCGAAAAAGAGCTCAGTTGGTTATCGTTCAATGAACGCGTGCTTCAGGAAGCGGCGGACAAATCTa                >  1:155635/1‑140 (MQ=255)                                                                                                                      cATAATATCCAGGCAGTGCCCCGTGAATAAAACGGAGTAAAAGTGGTAATGGGTCAGGAAAAGCTATACATCGAAAAAGAGCTCAGTTGGTTATCGTTCAATGAACGCGTGCTTCAGGAAGCGGCGGACAAATCTAAccc            >  1:916092/1‑140 (MQ=255)                                                                                                                      cATAATATCCAGGCAGTGCCCCGTGAATAAAACGGAGTAAAAGTGATAATGGGTCAGGAAAAGCTATACATCGAAAAAGAGCTCAGTTGGTTATCGTTCAATGAACGCGTGCTTCAGGAAGCGGCGGAc                       >  2:451476/1‑129 (MQ=255)                                                                                                                      cATAATATCCAGGCAGTGCCCCGTGAATAAAACGGAGTAAAAGTGATAATGGGTCAGGAAAAGCTATACATCGAAAAAGAGCTCAGTTGGTTATCGTTCAATGAACGCGTGCTTCAGGAAGCGGCGGAc                       <  1:451476/129‑1 (MQ=255)                                                                                                                       ataataTCCAGGCAGTGCCCCGTGAATAAAACGGAGTAAAAGTGGTAATGGGTCAGGAAAAGCTATACATCGAAAAAGAGCTCAGTTGGTTATCGTTCAATGAACGCGTGCTTCAGGAAGCGGCGGACAAATCTAACCCg           <  1:989942/140‑1 (MQ=255)                                                                                                                                aGGCAGTGCCCCGTGAATAAAACGGAGTAAAAGTGGTAATGGGTCAGGAAAAGCTATACATCGAAAAAGAGCTCAGTTGGTTATCGTTCAATGAACGCGTGCTTCAGGAAGCGGCGGACAAATCTAACCCGCTGATTGaa  >  1:656965/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                    CGTAATTAAAGCGCCAGCTCTGCCGCTGGCGTTTTTCAATTCACCTGTAAATCGCAAGCTCCAGCAGTTTTTTTCCCCCTTTTCTGGCATAGTTGGACATCTGCCAATATTGCTCGCCATAATATCCAGGCAGTGTCCCGTGAATAAAACGGAGTAAAAGTGGTAATGGGTCAGGAAAAGCTATACATCGAAAAAGAGCTCAGTTGGTTATCGTTCAATGAACGCGTGCTTCAGGAAGCGGCGGACAAATCTAACCCGCTGATTGAA  >  NC_000913/2622879‑2623145

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: ‑