Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,622,952	(T)7→6	intergenic (+80/‑92)	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,622,946	0	T	.	100.0%	120.4 / NA	29	intergenic (+74/‑98)	purN/ppk	phosphoribosylglycinamide formyltransferase 1/polyphosphate kinase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (14/15);  total (14/15)

TGGCTGGATGGTCAACGTCTGCCGCCGCAGGGCTACGCTGCCGACGAGTAATACCCCCGTAATTAAAGCGCCAGCTCTGCCGCTGGCGTTTTTCAATTCACCTGTAAATCGCAAGCTCCAGCAGTTTTTTTCCCCCTTTTCTGGCATAGTTGGACATCTGCCAATATTGCTCGCCATAATATCCAGGCAGTGTCCCGTGAATAAAACGGAGTAAAAGTGGTAATGGGTCAGGAAAAGCTATACATCGAAA  >  NC_000913/2622822‑2623071                                                                                                                             |                                                                                                                              tGGCTGGATGGTCAACGTCTGCCGCCGCAGGGCTACGCTGCCGACGAGTAATACCCCCGTAATTAAAGCGCCAGCTCTGCCGCTGGCGTTTTTCAATTCACCTGTAAATCGCAAGCTCCAGCAG‑TTTTTTCCCCCTTTTc                                                                                                               >  2:856697/1‑140 (MQ=255) tGGCTGGATGGTCAACGTCTGCCGCCGCAGGGCTACGCTGCCGACGAGTAATACCCCCGTAATTAAAGCGCCAGCTCTGCCGCTGGCGTTTTTCAATTCACCTGTAAATCGCAAGCTCCAGCAG‑TTTTTTCCCCCTTTTc                                                                                                               >  2:855030/1‑140 (MQ=255)           gTCAACGTCTGCCGCCGCAGGGCTACGCTGCCGACGAGTAATACCCCCGTAATTAAAGCGCCAGCTCTGCCGCTGGCGTTTTTCAATTCACCTGTAAATCGCAAGCTCCAGCAG‑TTTTTTCCCCCTTTTCTGGCATAGtt                                                                                                     >  2:618445/1‑140 (MQ=255)              aaCGTCTGCCGCCGCAGGGCTACGCTGCCGACGAGTAATACCCCCGTAATTAAAGCGCCAGCTCTGCCGCTGGCGTTTTTCAATTCACCTGTAAATCGCAAGCTCCAGCAG‑TTTTTTCCCCCTTTTCTGGCATAGTTGGa                                                                                                  >  2:963864/1‑140 (MQ=255)                cGTCTGCCGCCGCAGGGCTACGCTGCCGACGAGTAATACCCCCGTAATTAAAGCGCCAGCTCTGCCGCTGGCGTTTTTCAATTCACCTGTAAATCGCAAGCTCCAGCAG‑TTTTTTCCCCCTTTTCTGGCATAGTTGGACa                                                                                                <  1:963864/140‑1 (MQ=255)                cGTCTGCCGCCGCAGGGCTACGCTGCCGACGAGTAATACCCCCGTAATTAAAGCGCCAGCTCTGCCGCTGGCGTTTTTCAATTCACCTGTAAATCGCAAGCTCCAGCAG‑TTTTTTCCCCCTTTTCTGGCATAGTTGGACa                                                                                                <  1:980750/140‑1 (MQ=255)                cGTCTGCCGCCGCAGGGCTACGCTGCCGACGAGTAATACCCCCGTAATTAAAGCGCCAGCTCTGCCGCTGGCGTTTTTCAATTCACCTGTAAATCGCAAGCTCCAGCAG‑TTTTTTCCCCCTTTTCTGGCATAGTTGGACa                                                                                                <  1:636073/140‑1 (MQ=255)                          cgcAGGGCTACGCTGCCGACGAGTAATACCCCCGTAATTAAAGCGCCAGCTCTGCCGCTGGCGTTTTTCAATTCACCTGTAAATCGCAAGCTCCAGCAG‑TTTTTTCCCCCTTTTCTGGCATAGTTGGACATCTGCCAata                                                                                      >  2:30240/1‑140 (MQ=255)                          cgcAGGGCTACGCTGCCGACGAGTAATACCCCCGTAATTAAAGCGCCAGCTCTGCCGCTGGCGTTTTTCAATTCACCTGTAAATCGCAAGCTCCAGCAG‑TTTTTTCCCCCTTTTCTGGCATAGTTGGACATCTGCCAata                                                                                      <  1:496526/140‑1 (MQ=255)                            cAGGGCTACGCTGCCGACGAGTAATACCCCCGTAATTAAAGCGCCAGCTCTGCCGCTGGCGTTTTTCAATTCACCTGTAAATCGCAAGCTCCAGCAG‑TTTTTTCCCCCTTTTCTGGCATAGTTGGACATCTGCCGATAtt                                                                                    >  2:106224/1‑140 (MQ=255)                               ggCTACGCTGCCGACGAGTAATACCCCCGTAATTAAAGCGCCAGCTCTGCCGCTGGCGTTTTTCAATTCACCTGTAAATCGCAAGCTCCAGCAG‑TTTTTTCCCCCTTTTCTGGCATAGTTGGACATCTGCCAATATTGCt                                                                                 <  1:56200/140‑1 (MQ=255)                                gCTACGCTGCCGACGAGTAATACCCCCGTAATTAAAGCGCCAGCTCTGCCGCTGGCGTTTTTCAATTCACCTGTAAATCGCAAGCTCCAGCAG‑TTTTTTCCCCCTTTTCTGGCATAGTTGGACATCTGCCAATATTGCTc                                                                                <  1:855030/140‑1 (MQ=255)                                    cGCTGCCGACGAGTAATACCCCCGTAATTAAAGCGCCAGCTCTGCCGCTGGCGTTTTTCAATTCACCTGTAAATCGCAAGCTCCAGCAG‑TCTTTTCCCCCTTTTCTGGCATAGTTGGACATCTGCCAATATTGCTCGCCa                                                                            >  2:592287/1‑140 (MQ=255)                                     gCTGCCGACGAGTAATACCCCCGTAATTAAAGCGCCAGCTCTGCCGCTGGCGTTTTTCAATTCACCTGTAAATCGCAAGCTCCAGCAG‑TTTTTTCCCCCTTTTCTGGCATAGTTGGACATCTGCCGATATTGCTCGCCat                                                                           <  1:106224/140‑1 (MQ=255)                                       tGCCGACGAGTAATACCCCCGTAATTAAAGCGCCAGCTCTGCCGCTGGCGTTTTTCAATTCACCTGTAAATCGCAAGCTCCAGCAG‑TTTTTTCCCCCTTTTCTGGCATAGTTGGACATCTGCCAATATTGCTCGCCataa                                                                         <  1:25008/140‑1 (MQ=255)                                         ccGACGAGTAATACCCCCGTAATTAAAGCGCCAGCTCTGCCGCTGGCGTTTTTCAATTCACCTGTAAATCGCAAGCTCCAGCAG‑TTTTTTCCCCCTTTTCTGGCATAGTTGGACATCTGCCAATATTGCTCGCCataata                                                                       >  2:938090/1‑140 (MQ=255)                                          cgacgaGTAATACCCCCGTAATTAAAGCGCCAGCTCTGCCGCTGGCGTTTTTCAATTCACCTGTAAATCGCAAGCTCCAGCAG‑TTTTTTCCCCCTTTTCTGGCATAGTTGGACATCTGCCAATATTGCTCGCCATAatat                                                                      <  1:218063/140‑1 (MQ=255)                                                   aTACCCCCGTAATTAAAGCGCCAGCTCTGCCGCTGGCGTTTTTCAATTCACCTGTAAATCGCAAGCTCCAGCAG‑TTTTTTCCCCCTTTTCTGGCATAGTTGGACATCTGCCAATATTGCTCGCCATAATATCCAGGCAgt                                                             >  2:713398/1‑140 (MQ=255)                                                           gTAATTAAAGCGCCAGCTCTGCCGCTGGCGTTTTTCAATTCACCTGTAAATCGCAAGCTCCAGCAG‑TTTTTTCCCCCTTTTCTGGCATAGTTGGACATCTGCCAATATTGCTCGCCATAATATCCAGGCAGTGTCCCGTg                                                     >  2:228649/1‑140 (MQ=255)                                                                          gCTCTGCCGCTGGCGTTTTTCAATTCACCTGTAAATCGCAAGCTCCAGCAG‑TTTTTTCCCCCTTTTCTGGCATAGTTGGACATCTGCCAATATTGCTCGCCATAATATCCAGGCAGTGTCCCGTGAATAAAACGGAGTaa                                      <  1:452637/140‑1 (MQ=255)                                                                                   cTGGCGTTTTTCAATTCACCTGTAAATCGCAAGCTCCAGCAG‑TCTTTTCCCCCTTTTCTGGCATAGTTGGACATCTGCCAATATTGCTCGCCATAATATCCAGGCAGTGTCCCGTGAATAAAACGGAGTAAAAGTGGTaa                             <  1:592287/140‑1 (MQ=255)                                                                                       cGTTTTTCAATTCACCTGTAAATCGCAAGCTCCAGCAG‑TTTTTTCCCCCTTTTCTGGCATAGTTGGACATCTGCCAATATTGCTCGCCATAATATCCAGGCAGTGTCCCGTGAATAAAACGGAGTAAAAGTGGTAATggg                         <  1:713398/140‑1 (MQ=255)                                                                                           tttCAATTCACCTGTAAATCGCAAGCTCCAGCAG‑TTTTTTCCCCCTTTTCTGGCATAGTTGGACATCTGCCAATATTGCTCGCCATAATATCCAGGCAGTGTCCCGTGAATAAAACGGAGTAAAAGTGGTAATGGGTCAg                     >  1:687107/1‑140 (MQ=255)                                                                                              cAATTCACCTGTAAATCGCAAGCTCCAGCAG‑TTTTTTCCCCCTTTTCTGGCATAGTTGGACATCTGCCAATATTGCTCGCCATAATATCCAGGCAGTGTCCCGTGAATAAAACGGAGTAAAAGTGGTAATGGGTCAGGaa                  <  1:30240/140‑1 (MQ=255)                                                                                                      cTGTAAATCGCAAGCTCCAGCAG‑TTTTTTCCCCCTTTTCTGGCATAGTTGGACATCTGCCAATATTGCTCGCCATAATATCCAGGCAGTGTCCCGTGAATAAAACGGAGTAAAAGTGGTAAt                            >  2:129658/1‑122 (MQ=255)                                                                                                      cTGTAAATCGCAAGCTCCAGCAG‑TTTTTTCCCCCTTTTCTGGCATAGTTGGACATCTGCCAATATTGCTCGCCATAATATCCAGGCAGTGTCCCGTGAATAAAACGGAGTAAAAGTGGTAAt                            <  1:129658/122‑1 (MQ=255)                                                                                                              cGCAAGCTCCAGCAG‑TTTTTTCCCCCTTTTCTGGCATAGTTGGACATCTGCCAATATTGCTCGCCATAATATCCAGGCAGTGTCCCGTGAATAAAACGGAGTAAAAGTGGTAATGGGTCAGGAAAAGCTATACATCGaaa  >  2:246220/1‑140 (MQ=255)                                                                                                                cAAGCTCCAGCAG‑TTTTTTCCCCCTTTTCTGGCATAGTTGGACATCTGCCAATATTGCTCGCCATAATATCCAGGCAGTGTCCCGTGAATAAAACGGAGTAAAAGTGGTAATGGGTCAGTAAAAGCTATACATCGaaa  <  2:417479/138‑1 (MQ=255)                                                                                                                cAAGCTCCAGCAG‑TTTTTTCCCCCTTTTCTGGCATAGTTGGACATCTGCCAATATTGCTCGCCATAATATCCAGGCAGTGTCCCGTGAATAAAACGGAGTAAAAGTGGTAATGGGTCAGTAAAAGCTATACATCGaaa  >  1:417479/1‑138 (MQ=255)                                                                                                                             |                                                                                                                              TGGCTGGATGGTCAACGTCTGCCGCCGCAGGGCTACGCTGCCGACGAGTAATACCCCCGTAATTAAAGCGCCAGCTCTGCCGCTGGCGTTTTTCAATTCACCTGTAAATCGCAAGCTCCAGCAGTTTTTTTCCCCCTTTTCTGGCATAGTTGGACATCTGCCAATATTGCTCGCCATAATATCCAGGCAGTGTCCCGTGAATAAAACGGAGTAAAAGTGGTAATGGGTCAGGAAAAGCTATACATCGAAA  >  NC_000913/2622822‑2623071

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