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%	223.7 / NA	59	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 (30/29);  total (30/29)

GGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGACAACA  >  NZ_CP009273/1524729‑1524898                                                                                       |                                                                                    gggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGaacaac                                                                                  >  2:1227609/1‑90 (MQ=255) gggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGaacaac                                                                                  >  2:1950891/1‑90 (MQ=255)      aGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTaaaa                                                                             >  2:2095793/1‑90 (MQ=255)      aGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTaaaa                                                                             >  2:2238271/1‑90 (MQ=255)           aTGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATa                                                                        >  2:2219280/1‑90 (MQ=255)           aTGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATa                                                                        >  2:1600360/1‑90 (MQ=255)           aTGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATa                                                                        >  1:2093055/1‑90 (MQ=255)           aTGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATa                                                                        >  1:1329981/1‑90 (MQ=255)             gAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATaaa                                                                      <  1:1187238/90‑1 (MQ=255)                aaTTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATaaaaaa                                                                   <  1:1846460/90‑1 (MQ=255)                 aTTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAg                                                                  >  1:2805883/1‑90 (MQ=255)                    gtgATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGAc                                                               <  1:1654945/90‑1 (MQ=255)                           aaaTCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGt                                                        >  1:2432599/1‑90 (MQ=255)                             aTCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGa                                                      <  2:1078710/90‑1 (MQ=255)                             aTCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGa                                                      <  1:1494857/90‑1 (MQ=255)                              tCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAg                                                     <  2:2093055/90‑1 (MQ=255)                              tCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAg                                                     >  1:2161409/1‑90 (MQ=255)                                aCGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTa                                                   >  1:2885957/1‑90 (MQ=255)                                   aTTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTc                                                >  1:1526618/1‑90 (MQ=255)                                            tttGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGc                                       <  2:761380/90‑1 (MQ=255)                                             ttGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCt                                      >  2:1760693/1‑90 (MQ=255)                                             ttGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCt                                      >  1:2088111/1‑90 (MQ=255)                                             ttGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCt                                      >  1:972807/1‑90 (MQ=255)                                              tGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCtt                                     <  1:1600360/90‑1 (MQ=255)                                              tGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCtt                                     <  1:2509084/90‑1 (MQ=255)                                              tGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCtt                                     <  1:2219280/90‑1 (MQ=255)                                                 tttAACGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGAt                                  <  2:1743631/90‑1 (MQ=255)                                                 tttAACGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGAt                                  <  2:2523048/90‑1 (MQ=255)                                                  ttAACGCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGAtt                                 >  2:2188327/1‑90 (MQ=255)                                                       gCCTACAGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATAtt                            <  2:778736/90‑1 (MQ=255)                                                         cTACAGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTc                          >  2:1757662/1‑90 (MQ=255)                                                         cTACAGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTc                          >  1:1939484/1‑90 (MQ=255)                                                             aGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTg                      <  2:2885957/90‑1 (MQ=255)                                                             aGGTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTg                      <  2:2161409/90‑1 (MQ=255)                                                              ggTGATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGa                     <  1:2978202/90‑1 (MQ=255)                                                                tGATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGAcc                   <  1:1197491/90‑1 (MQ=255)                                                                tGATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGAcc                   <  2:2833862/90‑1 (MQ=255)                                                                tGATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGAcc                   <  1:2585157/90‑1 (MQ=255)                                                                tGATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGAcc                   <  1:418327/90‑1 (MQ=255)                                                                 gATCAAAATTTAACTTTTGGAACAACTAAAAATATAAAAAAGGACTCCAGGTAAGTATTCAAAGTTTGCTTGATTATAGTTCAGTGAccc                  >  1:1205396/1‑90 (MQ=255)                                                                 gATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTg                      >  1:1857590/1‑86 (MQ=255)                                                                 gATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTg                      <  2:1857590/86‑1 (MQ=255)                                                                 gATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGAccc                  >  2:1719146/1‑90 (MQ=255)                                                                 gATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGAccc                  <  1:1751305/90‑1 (MQ=255)                                                                 gATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGAccc                  >  1:570568/1‑90 (MQ=255)                                                                 gATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGAccc                  >  2:1559742/1‑90 (MQ=255)                                                                 gATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGAccc                  >  1:944735/1‑90 (MQ=255)                                                                 gATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGAccc                  >  2:1751305/1‑90 (MQ=255)                                                                 gATCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGAccc                  >  1:582549/1‑90 (MQ=255)                                                                  aTCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGAcccc                 <  1:1797538/90‑1 (MQ=255)                                                                   tCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGAt                                  <  1:225171/72‑1 (MQ=255)                                                                   tCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGAt                                  >  2:225171/1‑72 (MQ=255)                                                                   tCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCt                <  1:2184292/90‑1 (MQ=255)                                                                   tCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCt                <  1:2430030/90‑1 (MQ=255)                                                                   tCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCt                <  2:2805883/90‑1 (MQ=255)                                                                   tCAAAATTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCt                <  1:1227609/90‑1 (MQ=255)                                                                      aaaTTTAACCTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGat             >  2:45711/1‑90 (MQ=255)                                                                               cTTTGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGacaa    <  1:2188327/90‑1 (MQ=255)                                                                                tttGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGacaac   >  1:1271214/1‑90 (MQ=255)                                                                                tttGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGacaac   >  2:1875589/1‑90 (MQ=255)                                                                                 ttGGAACAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGacaaca  >  2:660269/1‑90 (MQ=255)                                                                                       |                                                                                    GGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGACAACA  >  NZ_CP009273/1524729‑1524898

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