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	97.7%	121.0 / NA	43	noncoding (1078/1542 nt)	BW25113_RS01010	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (28/14);  minor base G (1/0);  total (29/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGAC  >  NZ_CP009273/221248‑221415                                                                                        |                                                                                 ttGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGa                                                                                <  2:780736/90‑1 (MQ=255) ttGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGa                                                                                <  1:1157545/90‑1 (MQ=255) ttGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGa                                                                                <  1:20078/90‑1 (MQ=255)      aTCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGt                                                                           <  2:437468/90‑1 (MQ=255)       tCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGtt                                                                          >  2:152249/1‑90 (MQ=255)       tCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGtt                                                                          >  2:58909/1‑90 (MQ=255)       tCCACAGAACTTTCCAGAGATCGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGtt                                                                          >  2:481070/1‑90 (MQ=255)        ccACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTg                                                                         >  2:242429/1‑90 (MQ=255)               aCTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAg                                                                  >  2:1204701/1‑90 (MQ=255)               aCTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAg                                                                  <  2:196108/90‑1 (MQ=255)               aCTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGGCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAg                                                                  >  1:715096/1‑90 (MQ=255)                   tCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTccc                                                              <  2:60230/90‑1 (MQ=255)                        agaTGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAAc                                                         <  2:981614/90‑1 (MQ=255)                         gaTGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACg                                                        >  1:266938/1‑90 (MQ=255)                         gaTGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACg                                                        <  1:116295/90‑1 (MQ=255)                            ggATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACg                                                        <  1:349828/87‑1 (MQ=255)                            ggATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACg                                                        >  2:349828/1‑87 (MQ=255)                            ggATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAgc                                                     >  2:507368/1‑90 (MQ=255)                            ggATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAgc                                                     >  2:957663/1‑90 (MQ=255)                            ggATTGGTGCCTTCGGGAACTGGGAGACAGGGGCTGCATGGCTGTCGTCAGCTCGTGTTGGGAAATGTTGGGTTAAGTCCCGCAACGAgc                                                     >  1:1012234/1‑90 (MQ=14)                            ggATTGGTGACTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGaaa                                                                              >  2:1197498/1‑65 (MQ=17)                            ggATTGGTGACTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGaaa                                                                              <  1:1197498/65‑1 (MQ=14)                                 ggTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAAc                                                >  1:62852/1‑90 (MQ=37)                                               cTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTttgtt                                  >  2:1007582/1‑90 (MQ=255)                                                 gtgAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGc                                >  2:1050022/1‑90 (MQ=18)                                                 gtgAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGc                                >  1:1114493/1‑90 (MQ=18)                                                 gtgAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGc                                >  2:637807/1‑90 (MQ=18)                                                 gtgAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGc                                >  1:707539/1‑90 (MQ=18)                                                 gtgAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGc                                >  1:361811/1‑90 (MQ=18)                                                 gtgAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGc                                >  2:1091363/1‑90 (MQ=18)                                                 gtgAGACAGGTGCTGCATGGCTGTCCTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGc                                >  1:1007435/1‑90 (MQ=255)                                                   gagaCAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCa                              >  2:435822/1‑90 (MQ=18)                                                         ggTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTc                        >  1:231836/1‑90 (MQ=18)                                                         ggTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTc                        >  2:182246/1‑90 (MQ=18)                                                         ggTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTc                        >  2:163223/1‑90 (MQ=18)                                                              tgcATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAAAGAGCGAAACCCTTATCTTTTGTTGGCAGCGGTccggcc                   >  2:58158/1‑90 (MQ=11)                                                               gcATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTccggccg                  >  2:321698/1‑90 (MQ=18)                                                               gcATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTccggccg                  <  1:343940/90‑1 (MQ=25)                                                                 aTGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCggg                <  2:438030/90‑1 (MQ=25)                                                                       gtcgtcAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCa          >  1:842758/1‑90 (MQ=18)                                                                         cgtcAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCaaa        <  1:367960/90‑1 (MQ=18)                                                                              gCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGgaga   >  2:1099894/1‑90 (MQ=18)                                                                               cTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGAc  <  1:1050022/90‑1 (MQ=18)                                                                                        |                                                                                 TTGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGAC  >  NZ_CP009273/221248‑221415

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