Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,283,900	A→G	S61P (TCC→CCC)	purU ←	formyltetrahydrofolate deformylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,283,900	0	A	G	100.0%	111.8 / NA	32	S61P (TCC→CCC)	purU	formyltetrahydrofolate deformylase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (14/18);  total (14/18)

TCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGAATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTACGGTGATCAACAAATTCATTGTTCTGTACGATATTTAACTCGTGCTTGTAGCAAATATTGGTAATACGTGCGAT  >  NZ_CP009273/1283771‑1284023                                                                                                                                  |                                                                                                                            tcatcaACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTaa                                                                                                                     >  2:736350/1‑138 (MQ=255)  catcaACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTaaa                                                                                                                    >  1:319241/1‑138 (MQ=255)      aacaaATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATc                                                                                                                <  2:319241/138‑1 (MQ=255)                   ggCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGt                                                                                                  <  1:415391/139‑1 (MQ=255)                   ggCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGt                                                                                                  >  2:277594/1‑139 (MQ=255)                              tcttTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATaaaa                                                                                        >  1:173101/1‑138 (MQ=255)                                  tAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATaaaaaa                                                                                      <  1:277594/136‑1 (MQ=255)                                    gTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGcc                                                                                 <  2:663320/139‑1 (MQ=255)                                    gTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGcc                                                                                 <  1:771190/139‑1 (MQ=255)                                       aCCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCgg                                                                              >  1:782676/1‑139 (MQ=255)                                       aCCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCgg                                                                              >  1:796803/1‑139 (MQ=255)                                         cAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTa                                                                            <  2:584730/139‑1 (MQ=255)                                         cAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTa                                                                            <  1:667770/139‑1 (MQ=255)                                                      cGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTACGGTGATcaaca                                                                >  2:576486/1‑138 (MQ=255)                                                         cgacgaCCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTACGGTGATCAACAAAtt                                                            >  1:42898/1‑139 (MQ=255)                                                           acgacCGGCAGGATTCAGCTCACGCACGGGGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTACGGTGATCAACAAATTCa                                                          <  2:796803/139‑1 (MQ=255)                                                           acgacCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTACGGTGATCAACAAATTCa                                                          <  2:625640/139‑1 (MQ=255)                                                           acgacCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTACGGTGATCAACAAATTCa                                                          <  2:782676/139‑1 (MQ=255)                                                           acgacCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTACGGTGATCAACAAATTCa                                                          <  2:390162/139‑1 (MQ=255)                                                              acCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTACGGTGATCAACAAATTCATTg                                                       >  1:475420/1‑139 (MQ=255)                                                                      gATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTACGGTGATCAACAAATTCATTGTTCTGTAc                                               <  2:517364/139‑1 (MQ=255)                                                                        ttCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTACGGTGATCAACAAATTCATTGTTCTGTACg                                              >  2:317146/1‑138 (MQ=255)                                                                        ttCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTACGGTGATCAACAAATTCATTGTTCTGTACg                                              <  1:317146/138‑1 (MQ=255)                                                                             cTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTACGGTGATCAACAAATTCATTGTTCTGTACGATAtt                                         >  1:380658/1‑138 (MQ=255)                                                                                 cgcacgGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTACGGTGATCAACAAATTCATTGTTCTGTACGATATTTAACt                                    >  2:623711/1‑139 (MQ=255)                                                                                                 gCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTACGGTGATCAACAAATTCATTGTTCTGTACGATATTTAACTCGTGCTTGTAGCaaa                     <  2:614419/138‑1 (MQ=255)                                                                                                 gCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTACGGTGATCAACAAATTCATTGTTCTGTACGATATTTAACTCGTGCTTGTAGCaaa                     <  2:687616/138‑1 (MQ=255)                                                                                                           aTCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTACGGTGa                                                                      <  2:280570/79‑1 (MQ=255)                                                                                                           aTCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTACGGTGa                                                                      >  1:280570/1‑79 (MQ=255)                                                                                                            tCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTACGGTGATCAACAAATTCATTGTTCTGTACGATATTTAACTCGTGCTTGTAGCAAATATTGGTAAt           <  1:214151/137‑1 (MQ=255)                                                                                                                 aTCCGCCAGCAGGGAGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGACCGGTACTGTGATCAACAAATTCATTGTTCTGTACGACATTTAACTCGTGCTCGTAGCAAATATTCGTAATACGTGc     >  2:365418/1‑138 (MQ=255)                                                                                                                    cGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTACGGTGATCAACAAATTCATTGTTCTGTACGATATTTAACTCGTGCTTGTAGCAAATATTGGTAATACGTGCGAt  <  2:475420/138‑1 (MQ=255)                                                                                                                                  |                                                                                                                            TCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGAATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTACGGTGATCAACAAATTCATTGTTCTGTACGATATTTAACTCGTGCTTGTAGCAAATATTGGTAATACGTGCGAT  >  NZ_CP009273/1283771‑1284023

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