Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	4,073,163	A→T	D26V (GAT→GTT)	BW25113_RS20185 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	4,073,163	0	A	T	100.0%	111.4 / NA	30	D26V (GAT→GTT)	BW25113_RS20185	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (17/13);  total (17/13)

GAGATATTTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGATCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCAAA  >  NZ_CP009273/4073078‑4073243                                                                                      |                                                                                 gagaTATTTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTtctc                                                                              >  1:317526/1‑90 (MQ=255)    atatTTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAg                                                                           >  2:142583/1‑90 (MQ=255)        ttATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGc                                                                       <  1:220907/90‑1 (MQ=255)             aaCTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGt                                                                  <  1:201806/90‑1 (MQ=255)                tctctGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTGTCGAATGAGGTGATTAAGCAACTTTATGTTCTCGAGGTGCAGCgcagt                                                               >  1:305114/1‑86 (MQ=255)                 ctctGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATc                                                              >  1:121121/1‑90 (MQ=255)                        aGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCg                                                       >  2:452542/1‑90 (MQ=255)                              aGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGa                                                 <  2:317526/90‑1 (MQ=255)                                atatGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGAcc                                               <  1:430262/90‑1 (MQ=255)                                   tGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGTTGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAt                                            >  2:440173/1‑90 (MQ=255)                                   tGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAt                                            >  2:562197/1‑90 (MQ=255)                                   tGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAt                                            <  2:468725/90‑1 (MQ=255)                                   tGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAt                                            >  2:331125/1‑90 (MQ=255)                                    gggCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAtt                                           >  2:569926/1‑90 (MQ=255)                                    gggCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAtt                                           >  2:526296/1‑90 (MQ=255)                                    gggCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAtt                                           >  1:114977/1‑90 (MQ=255)                                    gggCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAtt                                           >  1:200003/1‑90 (MQ=255)                                     ggCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTa                                          <  1:544872/90‑1 (MQ=255)                                         gAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGa                                                 >  2:94325/1‑79 (MQ=255)                                         gAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGa                                                 <  1:94325/79‑1 (MQ=255)                                           aTTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGaa                                    <  1:331125/90‑1 (MQ=255)                                                aCTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGt                               <  1:578351/90‑1 (MQ=255)                                                 cTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTa                              <  2:470483/90‑1 (MQ=255)                                                   cGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGa                            >  2:301868/1‑90 (MQ=255)                                                                         aaGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCg      >  2:564036/1‑90 (MQ=255)                                                                           gCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCa    >  1:20480/1‑90 (MQ=255)                                                                             aaCTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCaaa  <  1:509321/90‑1 (MQ=255)                                                                             aaCTTAATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCaaa  <  1:376434/90‑1 (MQ=255)                                                                               cTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTa                              >  2:45874/1‑60 (MQ=255)                                                                               cTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTa                              <  1:45874/60‑1 (MQ=255)                                                                                      |                                                                                 GAGATATTTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGATCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCAAA  >  NZ_CP009273/4073078‑4073243

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