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	98.3%	187.4 / NA	60	A122T (GCA→ACA)	fmt	10‑formyltetrahydrofolate:L‑methionyl‑tRNA(fMet) N‑formyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (27/32);  total (27/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCGCCGAACTGCAGGCTGATGTTATGGTCGTCGTCGCCTATGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTGCACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACACCGGTGATATGCTCTATAAGCTCTCCTGCCCGATTACTGCAGAAGATACCAGTGGTA  >  NC_000913/3434445‑3434712                                                                                                                                     |                                                                                                                                        gTCGCCGAACTGCAGGCTGATGTTATGGTCGTCGTCGCCTATGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCaa                                                                                                                                   >  1:643481/1‑139 (MQ=255) gTCGCCGAACTGCAGGCTGATGTTATGGTCGTCGTCGCCTATGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCaa                                                                                                                                   >  2:817512/1‑139 (MQ=255) gTCGCCGAACTGCAGGCTGATGTTATGGTCGTCGTCGCCTATGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCaa                                                                                                                                   >  1:1405546/1‑139 (MQ=255)     ccGAACTGCAGGCTGATGTTATGGTCGTCGTCGCCTATGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCa                                                                                                                               >  1:99497/1‑139 (MQ=255)       gAACTGCAGGCTGATGTTATGGTCGTCGTCGCCTATGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAAc                                                                                                                             <  2:1058541/139‑1 (MQ=255)        aaCTGCAGGCTGATGTTATGGTCGTCGTCGCCTATGGTTTAATTCTGCCGAAAGTAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACg                                                                                                                            <  1:380328/139‑1 (MQ=255)        aaCTGCAGGCTGATGTTATGGTCGTCGTCGCCTATGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACg                                                                                                                            <  2:861304/139‑1 (MQ=255)        aaCTGCAGGCTGATGTTATGGTCGTCGTCGCCTATGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACg                                                                                                                            <  1:19559/139‑1 (MQ=255)         aCTGCAGGCTGATGTTATGGTCGTCGTCGCCTATGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGc                                                                                                                           >  1:530797/1‑139 (MQ=255)          cTGCAGGCTGATGTTATGGTCGTCGTCGCCTATGGTTTAATTCTGCCGAAAGCAGGGCGGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCt                                                                                                                          <  2:189958/139‑1 (MQ=255)            gCAGGCTGATGTTATGGTCGTCGTCGCCTATGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCa                                                                                                                        <  2:1193544/139‑1 (MQ=255)                  tGATGTTATGGTCGTCGTCGCCTATGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATgg                                                                                                                  >  1:346725/1‑139 (MQ=255)                    aTGTTATGGTCGTCGTCGCCTATGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGc                                                                                                                <  1:666993/139‑1 (MQ=255)                       ttATGGTCGTCGTCGCCTATGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCTCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCggg                                                                                                             <  1:817512/139‑1 (MQ=255)                                  gtcgCCTATGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGaaa                                                                                                  <  2:643481/139‑1 (MQ=255)                                         aTGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGgtgt                                                                                           <  2:312041/139‑1 (MQ=255)                                         aTGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGgtgt                                                                                           <  2:406814/139‑1 (MQ=255)                                         aTGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGgtgt                                                                                           <  2:530797/139‑1 (MQ=255)                                                aaTTCTGCCGAAAGGAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCATCGCTCACTATGGGCGGGTGATGCAGAAACGGGTGTGACCAtt                                                                                    <  2:1281997/139‑1 (MQ=255)                                                  ttCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTAt                                                                                  <  2:1405546/139‑1 (MQ=255)                                                     tGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTAAGGGTGGGTGACGCAGAAACTGGTGTGACCATTATGCa                                                                               >  2:978580/1‑139 (MQ=255)                                                      gCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCaa                                                                              >  2:884652/1‑139 (MQ=255)                                                             gCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATg                                                                       <  2:346725/139‑1 (MQ=255)                                                             gCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATg                                                                       <  1:884652/139‑1 (MQ=255)                                                              cAGTGCTGGAGATTCCGCGTCTTGGCTGTATCAACGTCCATGGTACACAGCTGCCACGCAGGCGCGGTGCTGCACCAATCCAACGCTCAATATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGt                                                                      <  2:39714/139‑1 (MQ=255)                                                                            ccGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACACCg                                                        <  1:702518/139‑1 (MQ=255)                                                                            ccGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACACCg                                                        <  1:1106509/139‑1 (MQ=255)                                                                             cgcgTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACTCTCACTATGGGCGGGTGATGCAGACACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACACCgg                                                       >  2:125568/1‑139 (MQ=255)                                                                                 tCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACa                                                                                                                                       >  1:388982/1‑55 (MQ=255)                                                                                 tCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACACCGGTGat                                                   >  1:757164/1‑139 (MQ=255)                                                                                 tCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACACCGGTGat                                                   <  2:742318/139‑1 (MQ=255)                                                                                     ggCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGgtgt                                                                                           >  1:1516770/1‑95 (MQ=255)                                                                                     ggCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGgtgt                                                                                           <  2:1516770/95‑1 (MQ=255)                                                                                           aTCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGaaa                                                                                                  <  2:1030315/82‑1 (MQ=255)                                                                                           aTCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGaaa                                                                                                  >  1:1030315/1‑82 (MQ=255)                                                                                                cGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCggg                                                                                                             <  2:816786/66‑1 (MQ=255)                                                                                                cGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCggg                                                                                                             >  1:816786/1‑66 (MQ=255)                                                                                                    cATGGTTCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACACCGGTGATATGCTCTATAAGCTCTCCt                                >  1:1308401/1‑139 (MQ=255)                                                                                                          tCACTGCTGCCACGCTGGCGCGGTGCTACACCAATcc                                                                                                                                <  1:1003505/37‑1 (MQ=255)                                                                                                          tCACTGCTGCCACGCTGGCGCGGTGCTACACCAATcc                                                                                                                                >  2:1003505/1‑37 (MQ=255)                                                                                                          tCACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACACCGGTGATATGCTCTATAAGCTCTCCTGCCCGa                          >  1:970373/1‑139 (MQ=255)                                                                                                           cACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCa                                                                                                                               >  2:1525242/1‑37 (MQ=38)                                                                                                           cACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCa                                                                                                                               <  1:1525242/37‑1 (MQ=255)                                                                                                           cACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACACCg                                                        >  2:938822/1‑108 (MQ=255)                                                                                                           cACTGCTGCCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACACCg                                                        <  1:938822/108‑1 (MQ=255)                                                                                                              tgctgcCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACACCGGTGATATGCTCTATAAGCTCTCCTGCCCGATTAc                      >  1:378346/1‑139 (MQ=255)                                                                                                                ctgcCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACACCGGTGATATGCTCTATAAGCTCTCCTGCCCGATTACTg                    <  2:970373/139‑1 (MQ=255)                                                                                                                 tgcCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGacac                                                          <  1:557331/100‑1 (MQ=255)                                                                                                                 tgcCACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGacac                                                          >  2:557331/1‑100 (MQ=255)                                                                                                                   ccACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCggg                                                                                                             >  2:337835/1‑47 (MQ=255)                                                                                                                   ccACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCggg                                                                                                             <  1:337835/47‑1 (MQ=255)                                                                                                                   ccACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACACCGGTGATATGCTCTATAAGCTCTCCTGCCCGATTACTGCag                 >  1:1342417/1‑139 (MQ=255)                                                                                                                   ccACGCTGGCGCGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACACCGGTGATATGCTCTATAAGCTCTCCTGCCCGATTACTGCag                 <  1:987601/139‑1 (MQ=255)                                                                                                                          ggCGCGGTGCTACACCAATCCAACGCTCACTATGGGCggg                                                                                                             >  1:804294/1‑40 (MQ=255)                                                                                                                          ggCGCGGTGCTACACCAATCCAACGCTCACTATGGGCggg                                                                                                             <  2:804294/40‑1 (MQ=255)                                                                                                                           gcgcGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACACCGGTg                                                     >  1:181689/1‑95 (MQ=255)                                                                                                                           gcgcGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACACCGGTg                                                     <  2:181689/95‑1 (MQ=255)                                                                                                                           gcgcGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACACCGGTGat                                                   <  1:21890/97‑1 (MQ=255)                                                                                                                           gcgcGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACACCGGTGat                                                   >  2:21890/1‑97 (MQ=255)                                                                                                                           gcgcGGTGCTACACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACACCGGTGATATGCTCTATAAGCTCTCCTGCCCGATTACTGCAGAAGATAcc         >  2:243156/1‑139 (MQ=255)                                                                                                                                  gctacACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACACCGGTGATATGCTCTATAAGCTCTCCTGCCCGATTACTGCAGAAGATACCAGTGGTa  <  1:125568/139‑1 (MQ=255)                                                                                                                                     |                                                                                                                                        GTCGCCGAACTGCAGGCTGATGTTATGGTCGTCGTCGCCTATGGTTTAATTCTGCCGAAAGCAGTGCTGGAGATGCCGCGTCTTGGCTGTATCAACGTTCATGGTTCACTGCTGCCACGCTGGCGCGGTGCTGCACCAATCCAACGCTCACTATGGGCGGGTGATGCAGAAACTGGTGTGACCATTATGCAAATGGATGTCGGTTTAGACACCGGTGATATGCTCTATAAGCTCTCCTGCCCGATTACTGCAGAAGATACCAGTGGTA  >  NC_000913/3434445‑3434712

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