Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	2,003,480	T→C	F60S (TTC→TCC)	yedK →	SOS response‑associated peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	2,003,480	0	T	C	100.0%	102.7 / NA	32	F60S (TTC→TCC)	yedK	SOS response‑associated peptidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/16);  total (16/16)

TTGGCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTTCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGTATG  >  NZ_CP009273/2003399‑2003559                                                                                  |                                                                                ttGGCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGat                                                                         >  1:156492/1‑90 (MQ=255)   ggCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGatat                                                                       <  2:313230/90‑1 (MQ=255)     cAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGc                                                                     <  1:610632/90‑1 (MQ=255)             aCGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGAt                                                             <  1:612238/90‑1 (MQ=255)               gTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGAtgg                                                           <  2:308710/90‑1 (MQ=255)                tCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGAtggt                                                          <  1:803422/90‑1 (MQ=255)                  gcgcCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGAtggtgg                                                        >  2:240643/1‑90 (MQ=255)                  gcgcCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGAtggtgg                                                        >  2:143026/1‑90 (MQ=255)                  gcgcCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGAtggtgg                                                        >  1:829733/1‑90 (MQ=255)                       gggAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATaa                                                   >  2:134447/1‑90 (MQ=255)                          aaCCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAAcc                                                <  2:641552/90‑1 (MQ=255)                               aaGTCCTACTGCTCAGTGAACGTGATGAACCCCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACcgccgc                                           <  1:187748/90‑1 (MQ=255)                                    cTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGAtt                                      >  1:12532/1‑90 (MQ=255)                                     tACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTa                                     >  1:664966/1‑90 (MQ=255)                                         gCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGc                                 <  2:156492/90‑1 (MQ=255)                                           tCAGTGACCGTGATGAACCCCTTCTTCGGGATCCGGTTTCCGGGGGATATGCCCCGGGATGGTGGGATAAACCGCCGCTGATTAACGccc                               <  1:770724/90‑1 (MQ=255)                                           tCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGccc                               >  1:552796/1‑90 (MQ=255)                                           tCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGccc                               <  2:664966/90‑1 (MQ=255)                                           tCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGccc                               <  1:270027/90‑1 (MQ=255)                                           tCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGccc                               <  2:537577/90‑1 (MQ=255)                                           tCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGccc                               >  2:545876/1‑90 (MQ=255)                                           tCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGccc                               <  2:650287/90‑1 (MQ=255)                                              gTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCcgcg                            >  2:82153/1‑90 (MQ=255)                                              gTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCcgcg                            >  2:505170/1‑90 (MQ=255)                                                   cGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGAtggtgg                                                        >  1:432412/1‑57 (MQ=255)                                                   cGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGAtggtgg                                                        <  2:432412/57‑1 (MQ=255)                                                            cacCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCa              <  1:12644/90‑1 (MQ=255)                                                              ccTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAAc            >  1:360438/1‑90 (MQ=255)                                                                 tCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAg         <  2:800019/90‑1 (MQ=255)                                                                   aTCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTc       >  1:500940/1‑90 (MQ=255)                                                                   aTCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTc       >  2:973778/1‑90 (MQ=255)                                                                        gATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGTATg  >  2:824677/1‑90 (MQ=255)                                                                                  |                                                                                TTGGCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTTCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGTATG  >  NZ_CP009273/2003399‑2003559

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