Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	221,335	T→C	noncoding (1078/1542 nt)	BW25113_RS01010 →	16S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	221,335	0	T	C	100.0%	114.4 / NA	39	noncoding (1078/1542 nt)	BW25113_RS01010	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/24);  total (15/24)

TTGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCC  >  NZ_CP009273/221248‑221419                                                                                        |                                                                                     ttGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGa                                                                                    <  2:499285/90‑1 (MQ=255)   gACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGaaa                                                                                  <  2:105378/90‑1 (MQ=255)      aTCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGt                                                                               <  1:418015/90‑1 (MQ=255)            aGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGtt                                                                         <  2:86777/90‑1 (MQ=255)             gagcTTGCCAGAGGTGGGTTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTa                                                                        <  1:598080/87‑1 (MQ=14)             gAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTa                                                                        <  1:954815/90‑1 (MQ=255)             gAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTa                                                                        <  1:56911/90‑1 (MQ=255)                  ttCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTcc                                                                   >  1:319094/1‑90 (MQ=255)                  ttCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTcc                                                                   >  2:705535/1‑90 (MQ=255)                       gagaTGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCaa                                                              >  2:783030/1‑90 (MQ=255)                            ggATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAgc                                                         <  1:387062/90‑1 (MQ=255)                               ttGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCa                                                      <  2:758521/90‑1 (MQ=255)                               ttGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCa                                                      <  1:900651/90‑1 (MQ=255)                                 ggTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGa                                                                                    >  2:855953/1‑58 (MQ=255)                                 ggTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGa                                                                                    <  1:855953/58‑1 (MQ=255)                                   tGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAAccc                                                  <  1:397959/90‑1 (MQ=18)                                      cTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACg                                                            <  2:926945/77‑1 (MQ=19)                                      cTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACg                                                            >  1:926945/1‑77 (MQ=25)                                             aaCTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTttg                                        <  2:75225/90‑1 (MQ=255)                                               cTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTttgtt                                      <  1:507179/90‑1 (MQ=255)                                               cTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTttgtt                                      <  1:265587/90‑1 (MQ=255)                                                 gtgAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGc                                    >  2:702221/1‑90 (MQ=18)                                                 gtgAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGc                                    >  2:161283/1‑90 (MQ=18)                                                    agaCAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAg                                 <  1:79504/90‑1 (MQ=25)                                                         ggTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTc                            >  2:535428/1‑90 (MQ=18)                                                                 aTGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCggg                    <  2:294541/90‑1 (MQ=25)                                                                 aTGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCggg                    <  1:135587/90‑1 (MQ=25)                                                                 aTGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCggg                    <  1:93965/90‑1 (MQ=18)                                                                 aTGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCggg                    <  1:472338/90‑1 (MQ=18)                                                                 aTGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCggg                    <  1:131199/90‑1 (MQ=25)                                                                 aTGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGAGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCggg                    <  1:548415/90‑1 (MQ=255)                                                                    gCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAAc                 >  1:87885/1‑90 (MQ=18)                                                                    gCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAAc                 >  1:619873/1‑90 (MQ=18)                                                                    gCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAAc                 >  2:738602/1‑90 (MQ=18)                                                                       gtcgtcAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCa              >  1:516192/1‑90 (MQ=18)                                                                         cgtcAGCTCGTGTTACGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCaaa            <  1:178365/90‑1 (MQ=255)                                                                              gCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGgaga       >  1:27571/1‑90 (MQ=18)                                                                                 cGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTg    >  2:23780/1‑90 (MQ=18)                                                                                 cGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTg    >  2:290401/1‑90 (MQ=18)                                                                                   tgttgcGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGcc  >  1:909796/1‑90 (MQ=18)                                                                                   tgttgcGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGcc  <  2:610796/90‑1 (MQ=18)                                                                                   tgttgcGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGcc  <  1:290401/90‑1 (MQ=18)                                                                                        |                                                                                     TTGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCC  >  NZ_CP009273/221248‑221419

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: ‑