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%	125.4 / NA	34	S61P (TCC→CCC)	purU	formyltetrahydrofolate deformylase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (11/23);  total (11/23)

GAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGAATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTACGGTGATCAACAAATTCATTGTTCTGTACGAT  >  NZ_CP009273/1283813‑1283981                                                                                        |                                                                                  gAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATc                                                                                <  2:185486/91‑1 (MQ=255)         aTCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTaaaaa                                                                        <  1:380050/91‑1 (MQ=255)           ccGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATc                                                                      >  1:413367/1‑91 (MQ=255)            cGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATcc                                                                     <  1:163530/91‑1 (MQ=255)            cGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATcc                                                                     <  1:843164/91‑1 (MQ=255)            cGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATcc                                                                     <  1:783251/91‑1 (MQ=255)            cGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATcc                                                                     <  1:486817/91‑1 (MQ=255)               cgacgaCCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCtt                                                                  >  1:190823/1‑91 (MQ=255)                   gacCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAg                                                              >  2:343735/1‑91 (MQ=255)                                     cacgcacgGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATaaacaa                                            >  2:682768/1‑91 (MQ=255)                                     cacgcacgGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATaaaaaa                                            >  1:349581/1‑91 (MQ=255)                                           cgGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCGGTTCCGTGCGCATAAAAAAGCGccc                                      <  2:413367/91‑1 (MQ=255)                                           cgGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGccc                                      <  1:672040/91‑1 (MQ=255)                                           cgGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGccc                                      <  1:378147/91‑1 (MQ=255)                                           cgGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGccc                                      <  2:195102/91‑1 (MQ=255)                                                ccTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTAc                                 <  1:404589/91‑1 (MQ=255)                                                ccTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTAc                                 <  2:438525/91‑1 (MQ=255)                                                ccTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTAc                                 <  1:738179/91‑1 (MQ=255)                                                ccTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGAGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTAc                                 <  1:377287/91‑1 (MQ=255)                                                         aaTGCGCTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTACGGTGATcaa                        >  1:763771/1‑91 (MQ=255)                                                            gcgcTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTACGGTGATcaacaa                     >  2:423207/1‑91 (MQ=255)                                                             cgcTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTACGGTGATCAACaaa                    <  2:150309/91‑1 (MQ=255)                                                             cgcTATCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTACGGTGATCAACaaa                    <  1:231182/91‑1 (MQ=255)                                                                  tCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTACGGTGATCAACAAATTCAt               <  2:190823/91‑1 (MQ=255)                                                                  tCGAGATCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTACGGTGATCAACAAATTCAt               <  1:557466/91‑1 (MQ=255)                                                                       aTCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTACGGTGATCAACAAATTCATTGTTc          <  2:135402/91‑1 (MQ=255)                                                                       aTCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTACGGTGATCAACAAATTCATTGTTc          <  1:702677/91‑1 (MQ=255)                                                                       aTCCGCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTACGGTGATCAACAAATTCATTGTTc          <  2:763771/91‑1 (MQ=255)                                                                           gCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTACGGTGATCAACAAATTCATTGTTCTGTa      >  2:460624/1‑91 (MQ=255)                                                                           gCCAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTACGGTGATCAACAAATTCATTGTTCTGTa      >  2:652749/1‑91 (MQ=255)                                                                            ccAGCAGGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTACGGTGATCAACAAATTCATTGTTCTGTAc     <  2:328744/91‑1 (MQ=255)                                                                               gcagGGTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTACGGTGATCAACAAATTCATTGTTCTGTACGat  >  2:53810/1‑91 (MQ=255)                                                                                  gggTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTAc                                 >  1:111301/1‑57 (MQ=255)                                                                                  gggTGGGATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTAc                                 <  2:111301/57‑1 (MQ=255)                                                                                        |                                                                                  GAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGAATCATTAAAAATCCCTTCCAGTTCCGTGCGCATAAAAAAGCGCCCGGTACGGTGATCAACAAATTCATTGTTCTGTACGAT  >  NZ_CP009273/1283813‑1283981

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