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	95.2%	108.1 / ‑0.8	42	noncoding (1078/1542 nt)	BW25113_RS01010	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (24/16);  total (26/16)
Fisher's exact test for biased strand distribution p-value = 5.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.88e-01

TTGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCC  >  NZ_CP009273/221248‑221419                                                                                        |                                                                                     ttGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGa                                                                                    <  2:572335/90‑1 (MQ=255) ttGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGa                                                                                    <  1:61972/90‑1 (MQ=255) ttGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGa                                                                                    <  2:792830/90‑1 (MQ=255)        ccACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTg                                                                             >  2:530837/1‑90 (MQ=255)        ccACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTg                                                                             >  2:109722/1‑90 (MQ=255)        ccACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTg                                                                             >  2:66254/1‑90 (MQ=255)        ccACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTg                                                                             <  1:547346/90‑1 (MQ=255)             gAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTa                                                                        <  1:234596/90‑1 (MQ=255)               aCTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAg                                                                      >  1:670263/1‑90 (MQ=255)                  ttCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTcc                                                                   >  2:25461/1‑90 (MQ=255)                  ttCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTcc                                                                   >  1:295096/1‑90 (MQ=255)                       gagaTGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCaa                                                              <  1:25461/90‑1 (MQ=255)                        agaTGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCCCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAAc                                                             <  2:315210/90‑1 (MQ=255)                         gaTGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACg                                                            <  1:598498/90‑1 (MQ=255)                         gaTGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACg                                                            <  1:606410/90‑1 (MQ=255)                         gaTGGATTGGTGCCATCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGCCAGCTCGTGTTGCTAAATGTTGGGTTAAGTCCCGCAACg                                                            <  2:515967/90‑1 (MQ=14)                            ggATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAgc                                                         >  1:524296/1‑90 (MQ=255)                            ggATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAgc                                                         <  1:713833/90‑1 (MQ=255)                            ggATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAAGGTTGGGTTAAGTCCCGCCACGAgc                                                         >  2:695139/1‑90 (MQ=14)                               ttGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCa                                                      <  1:193901/90‑1 (MQ=255)                               ttGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCa                                                      <  2:64230/90‑1 (MQ=255)                                 ggTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAAc                                                    >  2:551424/1‑90 (MQ=35)                                 ggTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAAc                                                    >  2:440436/1‑90 (MQ=37)                                 ggTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAAc                                                    >  2:317254/1‑90 (MQ=37)                                            gAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCtttt                                         >  1:808803/1‑90 (MQ=255)                                              aCTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTttgt                                       >  2:3406/1‑90 (MQ=255)                                                         ggTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTc                            >  2:496935/1‑90 (MQ=18)                                                         ggTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTc                            >  2:773726/1‑90 (MQ=18)                                                         ggTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCAAGCGGTc                            >  1:316956/1‑90 (MQ=255)                                                           tgctgcATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTccg                          >  1:596667/1‑90 (MQ=18)                                                               gcATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTccggccg                      >  1:614348/1‑90 (MQ=33)                                                                  tGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGa                   <  2:327199/90‑1 (MQ=25)                                                                    gCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGTAc                 >  2:484774/1‑90 (MQ=255)                                                                     cTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACt                <  2:524296/90‑1 (MQ=25)                                                                       gtcgtcAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCa              >  2:794715/1‑90 (MQ=18)                                                                       gtcgtcAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCa              >  1:628310/1‑90 (MQ=18)                                                                       gtcgtcAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCa              >  1:305292/1‑90 (MQ=18)                                                                         cgtcAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCaaa            >  1:340609/1‑90 (MQ=18)                                                                              gCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGgaga       >  2:631519/1‑90 (MQ=18)                                                                                 cGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCggg                    >  2:319367/1‑74 (MQ=19)                                                                                 cGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCggg                    <  1:319367/74‑1 (MQ=25)                                                                                 cGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTg    >  2:431330/1‑90 (MQ=18)                                                                                   tgttgtGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGcc  >  2:256340/1‑90 (MQ=33)                                                                                   tgttgcGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGcc  <  2:618524/90‑1 (MQ=18)                                                                                   tgttgcGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGcc  >  1:602707/1‑90 (MQ=18)                                                                                        |                                                                                     TTGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCC  >  NZ_CP009273/221248‑221419

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