Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,434,577	G→A	A122T (GCA→ACA)	fmt →	10‑formyltetrahydrofolate:L‑methionyl‑tRNA(fMet) N‑formyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,434,577	0	G	A	100.0%	118.9 / NA	38	A122T (GCA→ACA)	fmt	10‑formyltetrahydrofolate:L‑methionyl‑tRNA(fMet) N‑formyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/17);  total (21/17)

AACTGCAGGCTGATGTTATGGTCGTCGTCGCCTATGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTGCACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACACCGGTGATATGCTCTATAAGCTCTCCTGCCCGATTACTGCAGAAGATACCAGTG  >  NC_000913/3434452‑3434709                                                                                                                              |                                                                                                                                     aaCTGCAGGCTGATGTTATGGTCGTCGTCGCCTATGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACg                                                                                                                         <  2:993582/139‑1 (MQ=255)     gCAGGCTGATGTTATGGTCGTCGTCGCCTATGGTTTAATTCTGCTGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCa                                                                                                                     >  1:714455/1‑139 (MQ=255)       aGGCTGATGTTATGGTCGTCGTCGCCTATGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACt                                                                                                                   <  1:842215/139‑1 (MQ=255)                  aTGGTCGTCGTCGCCTATGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTg                                                                                                        >  1:841741/1‑139 (MQ=255)                    ggTCGTCGTCGCCTATGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGAt                                                                                                      <  2:841741/139‑1 (MQ=255)                       cgtcgtcgCCTATGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACg                                                                                                                         <  2:107363/117‑1 (MQ=255)                       cgtcgtcgCCTATGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAAGCCAACg                                                                                                                         >  1:107363/1‑117 (MQ=255)                        gtcgtcgCCTATGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCt                                                                                                                       >  1:34022/1‑118 (MQ=255)                        gtcgtcgCCTATGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCt                                                                                                                       <  2:34022/118‑1 (MQ=255)                          cgtcgCCTATGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGaa                                                                                                >  1:411569/1‑139 (MQ=255)                           gtcgCCTATGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGaaa                                                                                               <  2:411569/139‑1 (MQ=255)                           gtcgCCTATGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGaaa                                                                                               <  2:421376/139‑1 (MQ=255)                             cgCCTATGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACt                                                                                             <  2:714455/139‑1 (MQ=255)                               ccTATGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTgg                                                                                           >  1:708560/1‑139 (MQ=255)                                cTATGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGgt                                                                                          >  1:918538/1‑139 (MQ=255)                                              tGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTGTCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCa                                                                            >  1:863618/1‑139 (MQ=255)                                                ccGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACg                                                                                                                         >  2:302171/1‑92 (MQ=255)                                                ccGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACg                                                                                                                         <  1:302171/92‑1 (MQ=255)                                                     aGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGAt                                                                     >  1:209073/1‑139 (MQ=255)                                                          tGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCgg                                                                >  2:299160/1‑139 (MQ=255)                                                          tGCTGGAGATGCCCCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCgg                                                                >  1:159087/1‑139 (MQ=255)                                                                    gCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACAcc                                                      >  1:175595/1‑139 (MQ=255)                                                                    gCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACAcc                                                      >  2:217843/1‑139 (MQ=255)                                                                     ccGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACACCg                                                     <  1:217843/139‑1 (MQ=255)                                                                          tCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACACCGGTGat                                                <  2:209073/139‑1 (MQ=255)                                                                                   tATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACACCGGTGATATGCTCtat                                       >  2:666649/1‑139 (MQ=255)                                                                                   tATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACACCGGTGATATGCTCtat                                       >  2:772087/1‑139 (MQ=255)                                                                                            taagagacagctgctgCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACACCGGTGATATGCTCTATAAGCTCTcc                              <  2:733466/129‑1 (MQ=255)                                                                                               tGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACACCGGTGATATGCTCTATAAGCTCTCCTGc                           <  2:708560/139‑1 (MQ=255)                                                                                                        gctgcCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGaaa                                                                                               >  1:333644/1‑62 (MQ=255)                                                                                                        gctgcCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGaaa                                                                                               <  2:333644/62‑1 (MQ=255)                                                                                                           gcCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACACCGGTGATATGCTCTATAAGCTCTCCTGCCCGATTACTGCa               <  2:175595/139‑1 (MQ=255)                                                                                                           gcCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACACCGGTGATATGCTCTATAAGCTCTCCTGCCCGATTACTGCa               >  1:513254/1‑139 (MQ=255)                                                                                                            ccACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACACCGGTGATATGCTCTATAAGCTCTCCTGCCCGATTACTGCag              >  1:550840/1‑139 (MQ=255)                                                                                                            ccACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACACCGGTGATATGCTCTATAAGCTCTCCTGCCCGATTACTGCag              >  1:133503/1‑139 (MQ=255)                                                                                                                 cTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACACCGGTGATATGCTCTATAAGCTCTCCTGCCCGATTACTGCAGAAGAt         <  2:863618/139‑1 (MQ=255)                                                                                                                 cTGGCGCGGTGCTACACCAATCCAACGATCACTATGGGCGGGTGATGCAGAAACTGGTGTGAGCATTATACAAATGGATGTCGGTGTAGACAACGGTGATATGCTCTATAAGGTCTCCTGCCCGATTAcctgaatctta         >  1:733466/1‑129 (MQ=255)                                                                                                                        ggTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACAGCGGTGATATGCTCTATAAGCTCTCCTGCCCGATTACTGCAGAAGATACCAGTg  <  2:159087/139‑1 (MQ=255)                                                                                                                           gctacACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACACCGGTGATATg                                             <  1:743054/93‑1 (MQ=255)                                                                                                                           gctacACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACACCGGTGATATg                                             >  2:743054/1‑93 (MQ=255)                                                                                                                              |                                                                                                                                     AACTGCAGGCTGATGTTATGGTCGTCGTCGCCTATGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTGCACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACACCGGTGATATGCTCTATAAGCTCTCCTGCCCGATTACTGCAGAAGATACCAGTG  >  NC_000913/3434452‑3434709

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