Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,524,815	T→C	intergenic (‑54/‑28)	BW25113_RS25600 ← / → BW25113_RS07655	protein YncO/ISAs1 family transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,524,815	0	T	C	100.0%	90.5 / NA	28	intergenic (‑54/‑28)	BW25113_RS25600/BW25113_RS07655	protein YncO/ISAs1 family transposase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/13);  total (15/13)

GGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCT  >  NZ_CP009273/1524729‑1524884                                                                                       |                                                                      gggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGaacaac                                                                    >  1:1195418/1‑90 (MQ=255) gggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGaacaac                                                                    >  2:220423/1‑90 (MQ=255) gggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGaacaac                                                                    >  1:105712/1‑90 (MQ=255) gggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGaacaac                                                                    >  1:462515/1‑90 (MQ=255) gggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGaacaac                                                                    >  1:874841/1‑90 (MQ=255)      aGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTaaaa                                                               >  2:1275630/1‑90 (MQ=255)      aGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTaaaa                                                               >  2:1218858/1‑90 (MQ=255)      aGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTaaaa                                                               >  1:765441/1‑90 (MQ=255)      aGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTaaaa                                                               >  1:787410/1‑90 (MQ=255)           aTGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATa                                                          >  1:953027/1‑90 (MQ=255)           aTGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATa                                                          >  2:862480/1‑90 (MQ=255)           aTGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATa                                                          >  2:312231/1‑90 (MQ=255)           aTGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATa                                                          >  2:1339943/1‑90 (MQ=255)           aTGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATa                                                          >  2:1259754/1‑90 (MQ=255)           aTGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTCTGGAACAACTAAAAAGATa                                                          >  2:1181372/1‑90 (MQ=255)                  ttGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAgg                                                   <  2:462515/90‑1 (MQ=255)                  ttGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAgg                                                   <  2:1550442/90‑1 (MQ=255)                    gtgATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGAc                                                 <  1:935085/90‑1 (MQ=255)                    gtgATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGAc                                                 <  1:1339943/90‑1 (MQ=255)                     tgATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACg                                                <  2:1050115/90‑1 (MQ=255)                              tCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAg                                       >  1:409515/1‑90 (MQ=255)                                             ttGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCt                        <  1:26066/90‑1 (MQ=255)                                              tGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCtt                       <  1:1275630/90‑1 (MQ=255)                                              tGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCtt                       <  2:195490/90‑1 (MQ=255)                                              tGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCtt                       <  1:1401865/90‑1 (MQ=255)                                                      cGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAGAGTTTGCTTGATTatat               <  1:1575095/90‑1 (MQ=255)                                                         cTACAGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTc            >  2:917139/1‑90 (MQ=255)                                                              ggTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGa       <  1:920128/90‑1 (MQ=255)                                                                tGATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGAcc     <  2:105712/90‑1 (MQ=255)                                                                 gATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGAccc    >  1:1530018/1‑90 (MQ=255)                                                                 gATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGAccc    >  1:1322817/1‑90 (MQ=255)                                                                 gATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTAAATTTCAGTGAccc    >  2:1141107/1‑90 (MQ=255)                                                                   tCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCt  <  1:974391/90‑1 (MQ=255)                                                                                       |                                                                      GGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCT  >  NZ_CP009273/1524729‑1524884

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