Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,417,788	(C)6→7	coding (69/120 nt)	BW25113_RS26130 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,417,782	1	.	C	100.0%	229.8 / NA	64	G25G (GGA→GGG)	BW25113_RS26130	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (23/41);  total (23/41)

TGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCT‑CCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGCCCAACTTGAA  >  NZ_CP009273/1417700‑1417862                                                                                    |                                                                                 tGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCT‑ccccccc                                                                           >  2:245645/1‑89 (MQ=18)  ggTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCT‑cccccccg                                                                          >  1:2132123/1‑88 (MQ=18)  ggTCGTCGCCGGTTGCACATCAAAGCAGTCTGTCCGTCAGTGCGTCAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCT‑cccccccg                                                                          <  1:732756/90‑3 (MQ=11)    tcgtcgtcgGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCCCCACCGCCTCCGGCGTGGATAATGCAGCCT‑cccccccgac                                                                        <  2:2972537/90‑5 (MQ=255)    tcgtcgtcgGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCT‑cccccccgac                                                                        <  1:1207973/90‑5 (MQ=18)    tcgtcgtcgGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCT‑cccccccgac                                                                        <  1:3010235/90‑5 (MQ=18)            ggTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGac                                                                 <  1:407828/90‑1 (MQ=255)            ggTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGac                                                                 <  1:328925/90‑1 (MQ=255)                gCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCg                                                             >  1:1697312/1‑90 (MQ=255)                      cAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAAc                                                       <  2:918630/90‑1 (MQ=255)                      cAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAAc                                                       <  2:2115428/90‑1 (MQ=255)                      cAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAAc                                                       <  1:949532/90‑1 (MQ=255)                           ccgTCTGTCAGTAGACGGGTGAAGCCACCCCCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGat                                                  <  1:3010219/88‑1 (MQ=11)                           cAGTCTGTCAGTCAGTGCGTGAAGCCCCCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGat                                                  <  1:2103321/90‑1 (MQ=11)                           cAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGat                                                  <  2:1163713/90‑1 (MQ=255)                             gTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGCGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGatta                                                >  1:1342066/1‑90 (MQ=11)                               cTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGattatt                                              >  1:2859350/1‑90 (MQ=255)                                   cagtcagtGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCAc                                          <  2:1377598/90‑1 (MQ=255)                                     gtcagtGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCAccc                                        <  1:2212122/90‑1 (MQ=255)                                     gtcagtGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCAccc                                        <  1:2615619/90‑1 (MQ=255)                                       cagtGCGTGAAGCCCCCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTc                                      <  1:1811877/90‑1 (MQ=11)                                       cagtGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTc                                      >  2:2079590/1‑90 (MQ=255)                                       cagtGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTc                                      >  1:770030/1‑90 (MQ=255)                                       cagtGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTc                                      >  1:2476801/1‑90 (MQ=255)                                       cagtGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTc                                      >  2:895552/1‑90 (MQ=255)                                         gtGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCg                                    >  1:214088/1‑90 (MQ=255)                                           gCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCggg                                  >  2:2397355/1‑90 (MQ=255)                                           gCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCggg                                  >  2:2128232/1‑90 (MQ=255)                                            cGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGa                                 <  2:2943318/90‑1 (MQ=255)                                              tGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAAc                               <  2:73788/90‑1 (MQ=255)                                              tGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAAc                               <  1:1672685/90‑1 (MQ=255)                                               gAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGat                                                  >  2:1976375/1‑70 (MQ=255)                                               gAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGat                                                  <  1:1976375/70‑1 (MQ=255)                                               gAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACg                              >  1:2945814/1‑90 (MQ=255)                                                aaGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGa                             <  2:1891733/90‑1 (MQ=255)                                                 aGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGAc                            <  1:2397355/90‑1 (MQ=255)                                                  gCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACt                           <  2:817320/90‑1 (MQ=255)                                                  gCCACCACCGCCTCCGGCGCGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACt                           <  2:1342066/90‑1 (MQ=17)                                                   ccaccaCCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTg                          <  2:642641/90‑1 (MQ=255)                                                   ccaccaCCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTg                          <  2:1721731/90‑1 (MQ=255)                                                   ccaccaCCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTg                          >  1:1152583/1‑90 (MQ=255)                                                   ccaccaCCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTg                          >  1:1721731/1‑90 (MQ=255)                                                     accaccGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGt                        >  1:2657001/1‑90 (MQ=255)                                                        cccGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAat                     <  2:227955/89‑1 (MQ=255)                                                        accGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAat                     <  2:2476801/90‑1 (MQ=255)                                                        accGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAat                     <  1:895552/90‑1 (MQ=255)                                                        accGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAat                     <  1:2568219/90‑1 (MQ=255)                                                        accGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAat                     <  2:533860/90‑1 (MQ=255)                                                           gCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAatgat                  >  1:2661680/1‑90 (MQ=255)                                                            ccTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAatgatg                 >  2:2682387/1‑90 (MQ=255)                                                             cTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGc                <  2:770030/90‑1 (MQ=255)                                                             cTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGc                <  2:2657001/90‑1 (MQ=255)                                                             cTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGc                <  1:587778/90‑1 (MQ=255)                                                               ccGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGAc                            <  1:1784932/76‑1 (MQ=255)                                                               ccGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGAc                            >  2:1784932/1‑76 (MQ=255)                                                                cGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAgg             >  1:1614454/1‑90 (MQ=255)                                                                   cGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGccc          <  1:930166/90‑1 (MQ=255)                                                                   cGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGccc          <  2:1031156/90‑1 (MQ=255)                                                                   cGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGccc          <  2:2776515/90‑1 (MQ=255)                                                                   cGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGccc          <  2:2839573/90‑1 (MQ=255)                                                                   cGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGccc          <  1:1839951/90‑1 (MQ=255)                                                                   cGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGccc          <  1:2682387/90‑1 (MQ=255)                                                                      ggATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAat                     >  1:2424499/1‑76 (MQ=255)                                                                      ggATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAat                     <  2:2424499/76‑1 (MQ=255)                                                                      ggATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGCCCAAc       >  1:1364885/1‑90 (MQ=255)                                                                           aTGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAat                     >  1:2468039/1‑71 (MQ=255)                                                                           aTGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAat                     <  2:2468039/71‑1 (MQ=255)                                                                           aTGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGCCCAACTTGaa  <  2:1565161/90‑1 (MQ=255)                                                                           aTGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGCCCAACTTGaa  <  2:1803749/90‑1 (MQ=255)                                                                           aTGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGCCCAACTTGaa  >  2:1687733/1‑90 (MQ=255)                                                                                    |                                                                                 TGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCT‑CCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGCCCAACTTGAA  >  NZ_CP009273/1417700‑1417862

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