Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	214,743	A→C	Y174D (TAC→GAC)	proS ←	proline‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	214,743	0	A	C	100.0%	132.6 / NA	35	Y174D (TAC→GAC)	proS	proline‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (20/15);  total (20/15)

GGTGTCGGCTTGTACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTAGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGAAAC  >  NZ_CP009273/214666‑214826                                                                              |                                                                                    ggTGTCGGCTTGTACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAg                                                                         >  1:685086/1‑90 (MQ=255)     tCGGCTTGTACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGAt                                                                     <  2:69752/90‑1 (MQ=255)     tCGGCTTGTACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGAt                                                                     <  2:562572/90‑1 (MQ=255)      cGGCTTGTACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGAtt                                                                    >  1:616221/1‑90 (MQ=255)        gCTTGTACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTcc                                                                  <  2:23994/90‑1 (MQ=255)            gTACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAggg                                                                       <  2:994097/81‑1 (MQ=255)            gTACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAggg                                                                       >  1:994097/1‑81 (MQ=255)            gTACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTgag                                                              >  1:125462/1‑90 (MQ=255)              aCGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGaa                                                            <  1:617135/90‑1 (MQ=255)                ggCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGt                                                          <  2:1411310/90‑1 (MQ=255)                 gcgcGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTa                                                         >  1:121875/1‑90 (MQ=255)                 gcgcGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTa                                                         >  2:1573323/1‑90 (MQ=255)                 gcgcGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTa                                                         >  1:404056/1‑90 (MQ=255)                 gcgcGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTa                                                         >  1:17364/1‑90 (MQ=255)                  cgcgGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTAt                                                        <  2:1287037/90‑1 (MQ=255)                      gAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGaa                                                    <  2:49158/90‑1 (MQ=255)                      gAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGaa                                                    <  2:1512574/90‑1 (MQ=255)                       aaaTCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGaaa                                                   <  1:1154597/90‑1 (MQ=255)                             aGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTaa                                             <  2:872847/90‑1 (MQ=255)                                  cATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATc                                        >  1:477908/1‑90 (MQ=255)                                  cATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGC‑GGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCt                                       >  2:1186976/1‑90 (MQ=255)                                    tGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGACTCCTGAGAAGTATGGAAAGAGTAAGCATCtt                                      <  2:1208677/90‑1 (MQ=255)                                                ttttGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAAtt                          >  2:1552124/1‑90 (MQ=255)                                                          ggCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACg                >  2:468081/1‑90 (MQ=255)                                                            ccGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCa              >  2:1370483/1‑90 (MQ=255)                                                            ccGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCa              >  2:456650/1‑90 (MQ=255)                                                              gCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATg            >  2:1224036/1‑90 (MQ=255)                                                               cATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGa           >  1:1293957/1‑90 (MQ=255)                                                                    aTTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCg      >  1:1306105/1‑90 (MQ=255)                                                                     ttGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGa     >  1:939245/1‑90 (MQ=255)                                                                       gCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGaaa   >  2:494395/1‑90 (MQ=255)                                                                       gCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGaaa   >  2:993994/1‑90 (MQ=255)                                                                        cATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGAAAc  <  1:1370483/90‑1 (MQ=255)                                                                        cATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGAAAc  <  2:16776/90‑1 (MQ=255)                                                                        cATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGAAAc  <  1:1418888/90‑1 (MQ=255)                                                                              |                                                                                    GGTGTCGGCTTGTACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTAGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGAAAC  >  NZ_CP009273/214666‑214826

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