Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,288,215	A→G	L137P (CTG→CCG)	purU ←	formyltetrahydrofolate hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,288,215	0	A	G	100.0%	148.6 / NA	42	L137P (CTG→CCG)	purU	formyltetrahydrofolate hydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/17);  total (25/17)

TGATGATCTTATTCGGGAAGCGTGCCACAAATTCCGGCGTTAATACCCGCATATACTTCGCCAGCACCACGTAGTCAGGTTGATAAGCATCAATGGCATCCGCCATCTTTTGATCGTGCTCGTTGCGGGTTAACCCTTCATGGCTTACCAGCTCAAACGGAATATCAAAACGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACC  >  NC_000913/1288066‑1288355                                                                                                                                                      |                                                                                                                                             tgatgaTCTTATTCGGGAAGCGTGCCACAAATTCCGGCGTTAATACCCGCATATACTTCGCCAGCACCACGTAGTCAGGTTGATAAGCATCAATGGCATCCGCCATCTTTTGATCGTGCTCGTTGCGGGTTAACCCTTCATGGCTTACCa                                                                                                                                              <  1:439992/150‑1 (MQ=255)   atgatCTTATTCGGGAAGCGTGCCACAAATTCCGGCGTTAATACCCGCATATACTTCGCCAGCACCACGTAGTCAGGTTGATAAGCATCAATGGCATCCGCCATCTTTTGATCGTGCTCGTTGCGGGTTAACCCTTCATGGCTTACCGGCt                                                                                                                                           >  1:511809/1‑151 (MQ=255)    tgatCTTATTCGGGAAGCGTGCCACAAATTCCGGCGTTAATACCCGCATATACTTCGCCAGCACCACGTAGTCAGGTTGATAAGCATCAATGGCATCCGCCATCTTTTGATCGTGCTCGTTGCGGGTTAACCCTTCATGGCTTACCGGCTc                                                                                                                                          >  1:6481/1‑151 (MQ=255)     gatCTTATTCGGGAAGCGTGCCACAAATTCCGGCGTTAATACCCGCATATACTTCGCCAGCACCACGTAGTCAGGTTGATAAGCATCAATGGCATCCGCCATCTTTTGATCGTGCTCGTTGCGGGTTAACCCTTCATGGCTTACCGGCTc                                                                                                                                          >  1:29799/1‑150 (MQ=255)      atCTTATTCGGGAAGCGTGCCACAAATTCCGGCGTTAATACCCGCATATACTTCGCCAGCACCACGTAGTCAGGTTGATAAGCATCAATGGCATCCGCCATCTTTTGATCGTGCTCGTTGCGGGTTAACCCTTCATGGCTTACCGGCTCa                                                                                                                                         >  2:588633/1‑150 (MQ=255)          tATTCGGGAAGCGTGCCACAAATTCCGGCGTTAATACCCGCATATACTTCGCCAGCACCACGTAGTCAGGCTGATAAGCATCAATGGCATCCGCCATCTTTTGATCGTGCTCGTTGCGGGTTACCCCTTCATGGCTTACCGGCTCAAACgg                                                                                                                                    >  1:156386/1‑151 (MQ=255)              cGGGAAGCGTGCCACAAATTCCGGCGTTAATACCCGCATATACTTCGCCAGCACCACGTAGTCAGGTTGATAAGCATCAATGGCATCCGCCATCTTTTGATCGTGCTCGTTGCGGGTTAACCCTTCATGGCTTACCGGCTCAAACGGAat                                                                                                                                 >  2:194083/1‑150 (MQ=255)               gggAAGCGTGCCACAAATTCCGGCGTTAATCCCCGCCTATACTTCGCCAGCACCACGTAGTCAGGTTGATAAGCATCAATGGCATCCGCCATCTTTTGATCGTGCTCGTTGCGGGTTAACCCTTCATGGCTTACCGGCTCAAACGGAatat                                                                                                                               <  2:665022/151‑1 (MQ=255)                ggAAGCGTGCCACAAATTCCGGCGTTAATACCCGCATATACTTCGCCAGCACCACGTAGTCAGGTTGATAAGCATCAATGGCATCCGCCATCTTTTGATCGTGCTCGTTGCGGGTTAACCCTTCATGGCTTACCGGCTCAAACGGAATATc                                                                                                                              >  2:457824/1‑151 (MQ=255)                 gAAGCGTGCCACAAATTCCGGCGTTAATACCCGCATATACTTCGCCAGCACCACGTAGTCAGGTTGATAAGCATCAATGGCATCCGCCATCTTTTGATCGTGCTCGTTGCGGGTTAACCCTTCATGGCTTACCGGCTCAAACGGAATATCa                                                                                                                             <  1:599926/151‑1 (MQ=255)                     cGTGCCACAAATTCCGGCGTTAATACCCGCATATACTTCGCCAGCACCACGTAGTCAGGTTGATAAGCATCAATGGCATCCGCCATCTTTTGATCGTGCTCGTTGCGGGTTAACCCTTCATGGCTTACCGGCTCAAACGGAATATCAAAAc                                                                                                                         <  1:453067/151‑1 (MQ=255)                        gCCACAAATTCCGGCGTTAATACCCGCATATACTTCGCCAGCACCACGTAGTCAGGTTGATAAGCATCAATGGCATCCGCCATCTTTTGATCGTGCTCGTTGCGGGTTAACCCTTCATGGCTTACCGGCTCAAACGGAATATCAAAACGc                                                                                                                       <  1:313135/150‑1 (MQ=255)                        gCCACAAATTCCGGCGTTAATACCCGCATATACTTCGCCAGCACCACGTAGTCAGGTTGATAAGCATCAATGGCATCCGCCATCTTTTGATCGTGCTCGTTGCGGGTTAACCCTTCATGGCTTACCGGCTCAAACGGAATATCAAAACGCt                                                                                                                      <  2:501897/151‑1 (MQ=255)                          cccaAATTCCGGCGTTAATTCCCGCATATACTTCGCCAGCACCACGTAGTCAGGTTGATAAGCATCAATGGCATCCGTCATCTTTTGATCGTGCTCGTTGCGGGTTAACCCTTCATGGCTTACCGGCTCAAACGGAATATCAAAACGCTCa                                                                                                                    <  1:321461/149‑1 (MQ=255)                                  ccGGCGTTAATACCCGCATATACTTCGCCAGCACCACGTAGTCAGGTTGATAAGCATCAATGGCATCCGCCATCTTTTGATCGTGCTCGTTGCGGGTTAACCCTTCATGGCTTACCGGCTCAAACGGAATATCAAAACGCTCAACCagag                                                                                                             <  1:187528/150‑1 (MQ=255)                                    ggCGTTAATACCCGCATATACTTCGCCAGCACCACGTAGTCAGGTTGATAAGCATCAATGGCATCCGCCATCTTTTGATCGTGCTCGTTGCGGGTTAACCCTTCATGGCTTACCGGCTCAAACGGAATATCAAAACGCTCAACCAGAGAAc                                                                                                          >  2:247139/1‑151 (MQ=255)                                        ttAATACCCGCATATACTTCGCCAGCACCACGTAGTCAGGTTGATAAGCATCAATGGCATCCGCCATCTTTTGATCGTGCTCGTTGCGGGTTAACCCTTCATGGCTTACCGGCTCAAACGGAATATCAAAACGCTCAACCAGAGAACGTa                                                                                                       >  1:346900/1‑150 (MQ=255)                                              cccGCATATACTTCGCCAGCACCACGTAGTCAGGTTGATAAGCATCAATGGCATCCGCCATCTTTTGATCGTGCTCGTTGCGGGTTAACCCTTCATGGCTTACCGGCTCAAACGGAATATCAAAACGCTCAACCAGAGAACGTAAAGTATc                                                                                                >  2:130017/1‑151 (MQ=255)                                                    tataCTTCGCCAGCACCACGTAGTCAGGTTGATAAGCATCAATGGCATCCGCCATCTTTTGATCGTGCTCGTTGCGGGTTAACCCTTCATGGCTTACCGGCTCAAACGGAATATCAAAACGCTCAACCAGAGAACGTAAAGTATCGTGGtt                                                                                          >  1:558216/1‑151 (MQ=255)                                                      taCTTCGCCAGCACCACGTAGTCAGGTTGATAAGCATCAATGGCATCCGCCATCTTTTGATCGTGCTCGTTGCGGGTTAACCCTTCATGGCTTACCGGCTCAAACGGAATATCAAAACGCTCAACCAGAGAACGTAAAGTATCGTGGTTAc                                                                                        >  1:167260/1‑151 (MQ=255)                                                        cTTCGCCAGCACCACGTAGTCAGGTTGATAAGCATCAATGGCATCCGCCATCTTTTGATCGTGCTCGTTGCGGGTTAACCCTTCATGGCTTACCGGCTCAAACGGAATATCAAAACGCTCAACCAGAGAACGTAAAGTATCGTGGTTAcc                                                                                       >  2:253300/1‑150 (MQ=255)                                                              cAGCACCACGTAGTCAGGTTGATAAGCATCAATGGCATCCGCCATCTTTTGATCGTGCTCGTTGCGGGTTAACCCTTCATGGCTTACCGGCTCAAACGGAATATCAAAACGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACt                                                                                >  1:594350/1‑151 (MQ=255)                                                                  accacGTAGTCAGGTTGATAAGCATCAATGGCATCCGCCATCTTTTGATCGTGCTCGTTGCGGGTTAACCCTTCATGGCTTACCGGCTCAAACGGAATATCAAAACGCTCAACCAGAGAACGt                                                                                                        <  2:38848/123‑1 (MQ=255)                                                                  accacGTAGTCAGGTTGATAAGCATCAATGGCATCCGCCATCTTTTGATCGTGCTCGTTGCGGGTTAACCCTTCATGGCTTACCGGCTCAAACGGAATATCAAAACGCTCAACCAGAGAACGt                                                                                                        >  1:38849/1‑123 (MQ=255)                                                                            cAGGTTGATAAGCATCAATGGCATCCGCCATCTTTTGATCGTGCTCGTTGCGGGTTAACCCTTCATGGCTTACCGGCTCAAACGGAATATCAAAACGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGAc                                                                  >  2:666375/1‑151 (MQ=255)                                                                                  gATAAGCATCAATGGCATCCGCCATCTTTTGATCGTGCTCGTTGCGGGTTAACCCTTCATGGCTTACCGGCTCAAACGGAATATCAAAACGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCa                                                             <  1:418692/150‑1 (MQ=255)                                                                                    tAAGCATCAATGGCATCCGCCATCTTTTGATCGTGCTCGTTGCGGGTTAACCCTTCATGGCTTACCGGCTCAAACGGAATATCAAAACGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGgc                                                          >  1:297065/1‑151 (MQ=255)                                                                                            aaTGGCATCCGCCATCTTTTGATCGTGCTCGTTGCGGGTTAACCCTTCATGGCTTACCGGCTCAAACGGAATATCAAAACGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAg                                                            >  1:199417/1‑141 (MQ=255)                                                                                            aaTGGCATCCGCCATCTTTTGATCGTGCTCGTTGCGGGTTAACCCTTCATGGCTTACCGGCTCAAACGGAATATCAAAACGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAg                                                            <  2:199418/141‑1 (MQ=255)                                                                                            aaTGGCATCCGCCATCTTTTGATCGTGCTCGTTGCGGGTTAACCCTTCATGGCTTACCGGCTCAAACGGAATATCAAAACGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTa                                                   >  1:477430/1‑150 (MQ=255)                                                                                               ggCATCCGCCATCTTTTGATCGTGCTCGTTGCGGGTTAACCCTTCATGGCTTACCGGCTCAAACGGAATATCAAAACGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAAtt                                                <  1:592369/150‑1 (MQ=255)                                                                                                  aTCCGCCATCTTTTGATCGTGCTCGTTGCGGGTTAACCCTTCATGGCTTACCGGCTCAAACGGAATATCAAAACGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGc                                             >  2:19119/1‑150 (MQ=255)                                                                                                    ccGCCATCTTTTGATCGTGCTCGTTGCGGGTTAACCCTTCATGGCTTACCGGCTCAAACGGAATATCAAAACGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCtt                                           <  1:204412/150‑1 (MQ=255)                                                                                                         aTCTTTTGATCGTGCTCGTTGCGGGTTAACCCTTCATGGCTTACCGGCTCAAACGGAATATCAAAACGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTtcatc                                      >  1:443001/1‑150 (MQ=255)                                                                                                             tttGATCGTGCTCGTTGCGGGTTAACCCTTCATGGCTTACCGGCTCAAACGGAATATCAAAACGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATcaa                                    <  2:666993/148‑1 (MQ=255)                                                                                                                gATCGTGCTCGTTGCGGGTTAACCCTTCATGGCTTACCGGCTCAAACGGAATATCAAAACGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATc                              >  1:526231/1‑151 (MQ=255)                                                                                                                  tCGTGCTCGTTGCGGGTTAACCCTTCATGGCTTACCGGCTCAAACGGAATATCAAAACGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCg                             <  1:91912/150‑1 (MQ=255)                                                                                                                      gCTCGTTGCGGGTTAACCCTTCATGGCTTACCGGCTCAAACGGAATATCAAAACGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAg                        <  1:419320/151‑1 (MQ=255)                                                                                                                          gTTGCGGGTTAACCCTTCATGGCTTACCGGCTCAAACGGAATATCAAAACGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCa                     >  2:138888/1‑150 (MQ=255)                                                                                                                             gCGGGTTAACCCTTCATGGCTTACCGGCTCAAACGGAATATCAAAACGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATgc                 <  2:477433/151‑1 (MQ=255)                                                                                                                               gggTTAACCCTTCATGGCTTACCGGCTCAAACGGAATATCAAAACGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTtcatca                                     >  2:128775/1‑129 (MQ=255)                                                                                                                               gggTTAACCCTTCATGGCTTACCGGCTCAAACGGAATATCAAAACGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTtcatca                                     <  1:128771/129‑1 (MQ=255)                                                                                                                                            aTGGCTTACCGGCTCAAACGGAATATCAAAACGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGAcc  >  2:498873/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                             TGATGATCTTATTCGGGAAGCGTGCCACAAATTCCGGCGTTAATACCCGCATATACTTCGCCAGCACCACGTAGTCAGGTTGATAAGCATCAATGGCATCCGCCATCTTTTGATCGTGCTCGTTGCGGGTTAACCCTTCATGGCTTACCAGCTCAAACGGAATATCAAAACGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACC  >  NC_000913/1288066‑1288355

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