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%	165.2 / NA	51	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 (24/27);  total (24/27)

CCGAACTGCAGGCTGATGTTATGGTCGTCGTCGCCTATGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTGCACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACACCGGTGATATGCTCTATAAGCTCTCCTGCCCGATTACTGC  >  NC_000913/3434449‑3434695                                                                                                                                 |                                                                                                                       ccGAACTGCAGGCTGATGTTATGGTCGTCGTCGCCTATGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCa                                                                                                              <  2:833236/139‑1 (MQ=255)    aaCTGCAGGCTGATGTTATGGTCGTCGTCGCCTATGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACg                                                                                                           <  1:1078920/139‑1 (MQ=255)    aaCTGCAGGCTGATGTTATGGTCGTCGTCGCCTATGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACg                                                                                                           <  1:1061200/139‑1 (MQ=255)                 tGTTATGGTCGTCGTCGCCTATGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCa                                                                                                              <  1:1176711/123‑1 (MQ=255)                 tGTTATGGTCGTCGTCGCCTATGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCa                                                                                                              >  2:1176711/1‑123 (MQ=255)                 tGTTATGGTCGTCGTCGCCTATGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCa                                                                                                              <  1:423224/123‑1 (MQ=255)                 tGTTATGGTCGTCGTCGCCTATGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCa                                                                                                              >  2:423224/1‑123 (MQ=255)                 tGTTATGGTCGTCGTCGCCTATGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCg                                                                                              >  2:324751/1‑139 (MQ=255)                    tATGGTCGTCGTCGCCTATGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGt                                                                                           <  1:324751/139‑1 (MQ=255)                           gtcgtcgCCTATGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCa                                                                                                       >  2:439244/1‑120 (MQ=255)                           gtcgtcgCCTATGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCa                                                                                                       <  1:439244/120‑1 (MQ=255)                           gtcgtcgCCTATGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAg                                                                                    >  2:453375/1‑139 (MQ=255)                           gtcgtcgCCTATGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAg                                                                                    >  1:1355063/1‑139 (MQ=255)                                    tATGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGgtg                                                                           <  2:552709/139‑1 (MQ=255)                                     aTGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGgtgt                                                                          <  1:1289630/139‑1 (MQ=255)                                     aTGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGgtgt                                                                          <  2:897768/139‑1 (MQ=255)                                            aaTTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCAtt                                                                   >  1:646365/1‑139 (MQ=255)                                                  gCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCa                                                                     <  2:1283629/131‑1 (MQ=255)                                                  gCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCa                                                                     >  1:1283629/1‑131 (MQ=255)                                                  gCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCaa                                                             <  2:1355063/139‑1 (MQ=255)                                                   ccGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCaaa                                                            >  1:955240/1‑139 (MQ=255)                                                             tGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCa                                                                                                              >  1:813965/1‑79 (MQ=255)                                                             tGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCa                                                                                                              <  2:813965/79‑1 (MQ=255)                                                             tGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCa                                                                                                       >  1:837281/1‑86 (MQ=255)                                                             tGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCa                                                                                                       <  2:837281/86‑1 (MQ=255)                                                             tGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGt                                                                                           <  2:419044/98‑1 (MQ=255)                                                             tGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGt                                                                                           >  1:419044/1‑98 (MQ=255)                                                                       gCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACAcc                                        >  1:667092/1‑139 (MQ=255)                                                                        ccGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACACCg                                       <  1:488003/139‑1 (MQ=255)                                                                             tCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTAt                                                                 >  1:1098998/1‑108 (MQ=255)                                                                             tCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTAt                                                                 <  2:1098998/108‑1 (MQ=255)                                                                                tGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACACCGGTGATATg                               <  2:667092/139‑1 (MQ=255)                                                                                        tCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACACCGGTGATATGCTCTATaa                       >  2:842921/1‑139 (MQ=255)                                                                                            cGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCa                                                                                                              >  1:351863/1‑48 (MQ=255)                                                                                            cGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCa                                                                                                              <  2:351863/48‑1 (MQ=255)                                                                                              ttCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGc                                                                                      <  2:948655/70‑1 (MQ=255)                                                                                              ttCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGc                                                                                      >  1:948655/1‑70 (MQ=255)                                                                                                cATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCt                                                                                                         <  2:88939/49‑1 (MQ=255)                                                                                                cATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCt                                                                                                         >  1:88939/1‑49 (MQ=255)                                                                                                cATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGcctgtc                                                          <  2:1666517/96‑6 (MQ=255)                                                                                                cATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGcctgtc                                                          >  1:1666517/1‑91 (MQ=255)                                                                                                   ggTTCACTGCTGCCACGCTGACGCGGTGCTACACCAATCCAACg                                                                                                           <  1:56705/44‑1 (MQ=38)                                                                                                   ggTTCACTGCTGCCACGCTGACGCGGTGCTACACCAATCCAACg                                                                                                           >  2:56705/1‑44 (MQ=255)                                                                                                        aCTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGaaa                                                                                 <  1:479648/65‑1 (MQ=255)                                                                                                        aCTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGaaa                                                                                 >  2:479648/1‑65 (MQ=255)                                                                                                        aCTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGgtg                                                                           >  1:698850/1‑71 (MQ=255)                                                                                                        aCTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGgtg                                                                           <  2:698850/71‑1 (MQ=255)                                                                                                          tgctgcCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGa                                                                        <  1:1540419/72‑1 (MQ=255)                                                                                                          tgctgcCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGa                                                                        >  2:1540419/1‑72 (MQ=255)                                                                                                            ctgcCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACACCGGTGATATGCTCTATAAGCTCTCCTGCCCGATTACTg   <  2:955240/139‑1 (MQ=255)                                                                                                             tgcCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACACCGGTGATATGCTCTATAAGCTCTCCTGCCCGATTACTGc  >  1:1667839/1‑139 (MQ=255)                                                                                                                             tgctacACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCAt                                                                    >  1:1326398/1‑57 (MQ=255)                                                                                                                             tgctacACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCAt                                                                    <  2:1326398/57‑1 (MQ=255)                                                                                                                             tgctacACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACACCGGTGATATGCTCTATAAGCTCTCCTGCCCGATTACTg   >  1:1658750/1‑122 (MQ=255)                                                                                                                             tgctacACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACACCGGTGATATGCTCTATAAGCTCTCCTGCCCGATTACTg   <  2:1658750/122‑1 (MQ=255)                                                                                                                                 |                                                                                                                       CCGAACTGCAGGCTGATGTTATGGTCGTCGTCGCCTATGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTGCACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACACCGGTGATATGCTCTATAAGCTCTCCTGCCCGATTACTGC  >  NC_000913/3434449‑3434695

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