Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,428,805	G→A	F478F (TTC→TTT)	purF ←	amidophosphoribosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,428,805	0	G	A	100.0%	109.0 / NA	32	F478F (TTC→TTT)	purF	amidophosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (14/18);  total (14/18)

TTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGTTGAACACCGAGCATTCAAACTGCTGGATATCCGGATTTTCAGCGCGAACGGCGTCGATCAGATCGTTCAGATCCTGGAAAATCAACCCGT  >  NC_000913/2428669‑2428941                                                                                                                                         |                                                                                                                                         ttCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGaaaa                                                                                                                                       >  1:580813/1‑140 (MQ=255)  tCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGAAAAt                                                                                                                                      <  1:369698/140‑1 (MQ=255)           gTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGACGTAATCATTACGTAAAGTATGGAGAAAATCGAGGTAGc                                                                                                                             >  2:170584/1‑140 (MQ=255)            ttGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGAAAATCGAGGTAGcc                                                                                                                            <  2:52680/140‑1 (MQ=255)            ttGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGAAAATCGAGGTAGcc                                                                                                                            <  1:387066/140‑1 (MQ=255)                            cAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGAAAATCGAGGTAGCCCTGATCAACATCTTTg                                                                                                            <  1:417221/140‑1 (MQ=255)                                        cTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGAAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACg                                                                                                <  1:364143/140‑1 (MQ=255)                                        cTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGAAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACg                                                                                                <  1:483387/140‑1 (MQ=255)                                        cTCAGCGAAGGCATCATCCTTCCTTATGCATTTCGCGATTTTCCACTTTGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGAAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACg                                                                                                <  2:188077/140‑1 (MQ=255)                                           aGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGAAAATCAAGGAAGCCTTGATCAACATCTGTGGTGACGCAGACGCCg                                                                                             >  2:323723/1‑140 (MQ=255)                                                   catcatCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGAAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGTTGAacac                                                                                     >  1:544933/1‑140 (MQ=255)                                                                 aTGCATTTCGAGATTTTTCACTTCGTTCTGACGTTTCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGAAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGTTGAACACCGAGCATTCAAAct                                                                       <  2:249439/140‑1 (MQ=255)                                                                  tGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGAAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGTTGAACACCGAGCATTCAAActg                                                                      >  1:596851/1‑140 (MQ=255)                                                                  tGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGAAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGTTGAACACCGAGCATTCAAActg                                                                      >  1:561702/1‑140 (MQ=255)                                                                                           tCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGAAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGTTGAACACCGAGCATTCAAACTGCTGGATATCCGGATTTTCAGCGCGa                                             >  1:51540/1‑140 (MQ=255)                                                                                              gACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGAAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGTTGAACACCGAGCATTCAAACTGCTGGATATCCGGATTTTCAGCGCGAACg                                          <  2:544933/140‑1 (MQ=255)                                                                                                cGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGAAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGTTGAACACCGAGCATTCAAACTGCTGGATATCCGGATTTTCAGCGCGAACGGc                                        <  1:530868/140‑1 (MQ=255)                                                                                                 gTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGAAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGTTGAACACCGAGCATTCAAACTGCTGGATATCCGGATTTTCAGCGCGAACGGCg                                       <  2:580813/140‑1 (MQ=255)                                                                                                    gCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGAAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGTTGAACACCGAGCATTCAAACTGCTGGATATCCGGATTTTCAGCGCGAACGGCGTCg                                    >  1:610325/1‑140 (MQ=255)                                                                                                          cTTTGGCGTCATCATTACGTAATGTCTCGGGAAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGTTGAACACCGAGCATTCAAACTGCTGGATATCCGGATTTTCAGCGCGAACGGCGTCGATCAGa                              <  1:323723/140‑1 (MQ=255)                                                                                                          cTTTGGCGTCATCATTACGTAACGTATCGAGAAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGTTGAACACCGAGCATTCAAACTGCTGGATATCCGGATTTTCAGCGCGAACGGCGTCGATCAGa                              <  2:561702/140‑1 (MQ=255)                                                                                                          cTTTGGCGTCATCATTACGTAACGTATCGAGAAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGTTGAACACCGAGCATTCAAACTGCTGGATATCCGGATTTTCAGCGCGAACGGCGTCGATCAGa                              <  2:596851/140‑1 (MQ=255)                                                                                                          cTTTGGCGTCATCATTACGTAACGTATCGAGAAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGTTGAACACCGAGCATTCAAACTGCTGGATATCCGGATTTTCAGCGCGAACGGCGTCGATCAGa                              <  1:363249/140‑1 (MQ=255)                                                                                                          cTTTGGCGTCATCATTACGTAACGTATCGAGAAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGTTGAACACCGAGCATTCAAACTGCTGGATATCCAGATGTTGAGCGCGAACGGCGTCGATCAGa                              >  1:183488/1‑140 (MQ=255)                                                                                                             tGGCGTCATCATTACGTAACGTATCGAGAAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGTTGAACACCGAGCATTCCAACTGCTGGATAACCGGATTTTCAGTGCGAACGGCGTCGATCAGATCg                           >  2:232456/1‑140 (MQ=255)                                                                                                               gCGTCATCATTACGTAACGTATCGAGAAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGTTGAACACCGAGCATTCAAACTGCTGGATATCCGGATTTTCAGCGCGAACGGCGTCGa                                   >  2:114049/1‑130 (MQ=255)                                                                                                               gCGTCATCATTACGTAACGTATCGAGAAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGTTGAACACCGAGCATTCAAACTGCTGGATATCCGGATTTTCAGCGCGAACGGCGTCGa                                   <  1:114049/130‑1 (MQ=255)                                                                                                                    atcatTACGTAACGTATCGAGAAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGTTGAACACCGAGCATTCAAACTGCTGGATATCCGGATTTTCAGCGCGGACAGCGACGCTCAGATCGTTCAGAt                    >  1:312925/1‑140 (MQ=255)                                                                                                                         tACGTAACGTATCGAGAAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGTTGAACACCGAGCATTCAAACTGCTGGATATCCGGATTTTCAGCGCGAACGGCGTCGATCAGATCGTTCAGATCCTgg               >  2:82554/1‑140 (MQ=255)                                                                                                                                 gtaTCGAGAAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGTTGAACACCGAGCATTCAAACTGCTGGATATCCGGATTTTCAGCGCGAACGGCGTCGATCAGATCGTTCAGATCCTGGAAAATCaa       <  1:545005/140‑1 (MQ=255)                                                                                                                                   aTCGAGAAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGTTGAACACCGAGCATTCAAACTGCTGGATATCCGGATTTTCAGCGCGAACGGCGTCGATCAGATCGTTCAGATCCTGGAAAATCAAcc     >  2:545419/1‑140 (MQ=255)                                                                                                                                   aTCGAGAAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGTTGAACACCGAGCATTCAAACTGCTGGATATCCGGATTTTCAGCGCGAACGGCGTCGATCAGATCGTTCAGATCCTGGAAAATCAAcc     >  1:596885/1‑140 (MQ=255)                                                                                                                                    tCGAGAAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGTTGAACACCGAGCATTCAAACTGCTGGATATCCGGATTTTCAGCGCGAACGGCGTCGATCAGATCGTTCAGATCCTGGAAAATCAAccc    <  1:610814/140‑1 (MQ=255)                                                                                                                                      gagaAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGTTGAACACCGAGCATTCAAACTGCTGGATATCCGGATTTTCAGCGCGAACGGCGTCGATCAGATCGTTCAGATCCTGGAAAATCAACCCGt  <  2:610325/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                         TTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACGCCGTTGAACACCGAGCATTCAAACTGCTGGATATCCGGATTTTCAGCGCGAACGGCGTCGATCAGATCGTTCAGATCCTGGAAAATCAACCCGT  >  NC_000913/2428669‑2428941

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