Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NZ_CP034943	2,981,253	C→T	37.7%	noncoding (1267/1550 nt)	EO946_RS15435 ←	16S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP034943	2,981,253	0	C	T	37.7%	‑2.3 / 40.6	44	noncoding (1267/1550 nt)	EO946_RS15435	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (6/21);  new base T (2/14);  total (9/35)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTACAAGGCCCGGGAACGTATTCACCGCGGCATGCTGATCCGCGATTACTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTCGCTGCCC‑TTTGTTCTGTCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCACCTTCCTCCGGTTTGTCACCGGCAGTCACCTTAGAGTGCCCAACTGAATGCTGGC  >  NZ_CP034943/2981117‑2981385                                                                                                                                         |                                                                                                                                      tgttctaGGCCCGGGATCGTATTCACCGCGGGATGCTGATCAGCGATTACTAGCGATTCCAGCTTCACGCAGTAGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTCGCt                                                                                                                                    <  1:1606662/134‑1 (MQ=1) tgtcccaGGCCCGGGAACGTATTCACCGCGGCATGCTGAGCCGCGATTACTCGCGATTCCAGCTTCACGCAGTCGCGGTGCAGACTGCGCTCCGAACTGAGAACAGATTTGGGGGATTGTCTTAACCTCGCGGTTTCGCt                                                                                                                                    <  2:3379781/134‑1 (MQ=1) tggaCAGGGCCCGGGAACGTATTCACCGCGGCATGCTGATCCGCGATTACTAGCGAGTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGAGTTGTGGGATGGGATTAACCTCGAGGTTTCGCt                                                                                                                                    >  2:4633466/4‑140 (MQ=1) tggaCAAGGCCCGGGAACGTATTCACCGCGGCCTGCTGATCCGCGATTACTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTCGCt                                                                                                                                    <  2:330691/137‑1 (MQ=1) tggaCAAGGCCCGGGAACGTATTCACCGCGGCATGCTGATCCGCGATTACTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTTGCt                                                                                                                                    <  1:1945960/137‑1 (MQ=1) cccaccaGGCCCGGGAACGAATTCACCGCGCCATGCGGCACCCCCATTACTACCGCGGCCAGCTTCACGCAGTCGACTGGCAGACTGCGATCCGAACTGAGACCCGATTTGTGGGATTGGCTCAACCTCGCGGTTTTGCt                                                                                                                                    <  2:1760886/134‑1 (MQ=0)      aacgCCCGGGAACGTAGTCACCGCGGCATGCTGATGCGCGATTACTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTTGCTGCCC‑t                                                                                                                              <  2:2907548/137‑1 (MQ=0)       agacccGGGAACGTATTCACCGCGGCATGCTGATCCGCGATTACTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTTGCTGCCC‑tt                                                                                                                             >  1:4485087/4‑140 (MQ=1)        atacccggAACCTATTCACCGCTGAATGATAAACCGCGATTACTATCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTCGCTGCCC‑ttt                                                                                                                            >  1:951396/7‑140 (MQ=1)         gcgcGGGAACGTATTCACCGCGGCATGCTGATCCGCGATTACTATCGATTCCAGCTACACGCAGTCGAGTTGCAGTCTGCGATCCGAACTGAGAACAGATTGGTGGGATTGGCTTAACCTCGCGGTTTCGCTGACC‑TTTg                                                                                                                           >  2:1656141/4‑140 (MQ=1)         gcgcGGGAACGTATTCACCGCGGCATGCTGATCCGCGATTACTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTAAACCTCGCGGTTTCGCTGCCC‑TTTg                                                                                                                           >  1:3588060/4‑140 (MQ=1)           gcgggttcGTATTGAGGGCGGTATGCTGATGCGTTTTTAGTTTCGTTTCCAGATTCTCGCAGTCGTTTTGCAGTGTGGGATCCGATCTGTTAAGAGTTTTGTGTGATTGGCTATTCCTCGCGGTTTCGCTGCCC‑TTTGtt                                                                                                                         <  2:2542549/133‑1 (MQ=0)             gggAACGTATTCACCGCGGCATGCTGATCCGCGATTACTAGCGATTCCAGCTTCACGCAGGCGAGTTGCAGACTGCGATCCGAACTGAGAACCGATGTGGCGGATTGGCTGACCCTCGCGGTTTCGCTGCCCGTTTGtct                                                                                                                        >  2:4321796/1‑138 (MQ=1)             ggaaaCGTATTCACCGCGGCATGCTGATCCGCGATTACTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTTGCTGCCC‑TTTGTTCt                                                                                                                       >  2:1626959/4‑140 (MQ=1)               gaaagaaTTCACCGCGGCATGCTGATCCGAGATTAATTGCGATTCCAGCTTCACTCAGTCGCTTTGCAGACTGCGTTCCGAACTGAGAACAGATTTGTGGGATTGGGTGACCCTCGCGGTTTTGCTGCCC‑TTTGTTCTGt                                                                                                                     <  1:2748457/134‑1 (MQ=1)                           cggggCATGCTGATCCGCGATTACTCGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGCTCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTTGCTGCCC‑TTTGTTCTGTCCATTGTAGCAc                                                                                                         <  2:3816193/137‑1 (MQ=0)                           cgaggCATGCTGATCCGCGATTACTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTCGCTGCCC‑TTTGTTCTGTCCATTGTAGCAc                                                                                                         <  1:1721203/137‑1 (MQ=1)                                caggCTGATCCGCGATTACTCGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGACCTGAGAACAGATTTGTGGGAGGGGCTTAACCTCGCGGTTTTGCTGCCC‑TTGGTTCTGTCCATTGTAGCACgtgtg                                                                                                    <  1:2205725/137‑1 (MQ=1)                                     tgctCCGCGATTACTAGCGATTCCATGTTCACGCAGTCGACTTGCAGACTGCGATCCGAACTGATAACAGTTTTGTGGGAGTGGCTTAACCTCGCGGTTTCGCTGCCC‑TTTGTTCTGTCCATTGTAGCACGTGTGTAGcc                                                                                               <  1:1238921/137‑1 (MQ=1)                                     tgctCCGCGATTACTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTCGCTGCCC‑TTTGTTCTGTCCATTGTAGCACGTGTGTAGcc                                                                                               <  2:2773620/137‑1 (MQ=1)                                     tcagccccgAGTACTCGCGAGTCCCGCTACACGCAGTCGAGGTGCAGCCTGCGCTCCGACCCCACAACAGATTTGTCGGCTGCGCGTAACCTCGCGGTTTTGCTGCCC‑TTAGTTCGGACCATTGTAGCACGTGTGTAGcc                                                                                               <  2:2806592/133‑1 (MQ=0)                                      gagccGCGATTACTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTCGCTGCCC‑TTTGTTCTGTCCATTGTAGCACGTGTGTAGccc                                                                                              <  1:299838/137‑1 (MQ=1)                                      gaaccGCGATTACTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATGGGCTTAACCTCGCGGTTTTGCTGCCC‑TTTGTTCTGTCCATTGTAGCACGTGTGTAGccc                                                                                              <  2:1021286/137‑1 (MQ=0)                                            cgtttACTAGCGATTCCAGCTTCACGCAGTAGAGTTGCAGACTGCGATCCGTACTGAGAAGAGATTTGTGGGATTGGCTAAACCTCGCGGTTTGGCTGCGC‑TGTGTTCTGTCCAGTGTCGGACGTGTGTTGCCCCGGTcg                                                                                        >  1:4248420/4‑139 (MQ=1)                                             gagtACTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTTGCTGCCC‑TTTGTTCTGTCCATTGTAGCACGTGTGTAGCCCAGGTCAt                                                                                       <  2:4451224/137‑1 (MQ=1)                                             gagtACGAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCGCGCGGTTTTGCTGCCC‑TTTGTTCTGTCCATTGTAGCACGTGTGTAGCCCAGGTCAt                                                                                       <  1:2798888/137‑1 (MQ=0)                                              ctgacgaGCCATTCCAGCTTCCCGCAGTCGAGTTGCAGCCTGCGCTCCGACCTGAGAACAGAGTTGGGGCATTGGCTTAACCTCGCGGTTTTGCTGCCC‑TGTGTTCGGTCCATTGTAGCCCGTGTGTAGCCCAGGTCATa                                                                                      <  2:4101946/134‑1 (MQ=1)                                              atgaCTAGCGATTCCAGCTTCACGCAGTCGAGTGGCAGCCTGCGCTCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTTGCTGCCC‑TTTGTTCTGTCCATGGTCGCACGTGTGTAGCCCAGGTCATa                                                                                      <  1:799644/137‑1 (MQ=1)                                              atgaCTAGCGAGTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTCGCTGCCC‑TTTGTTCTGTCCATTGTAGCACGTGTGTAGCCCAGGTCATa                                                                                      <  1:4238855/137‑1 (MQ=1)                                              atgaCTAGCGAGTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTCGCTGCCC‑TTTGTTCTGTCCATGGTAGCACGTGTGTAGCCCAGGTCATa                                                                                      <  1:1718080/137‑1 (MQ=0)                                                       gagtCCCGCTTCCCGCAGTCGAGTTGCAGACTGCGCTCCGAACTGAGAACACATGTGAGGGATTGGCTTAACCTCGCGGTTTTGCTGCCC‑TTTGTTCTGTCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCatg                                                                             <  2:4865992/137‑1 (MQ=1)                                                       gagtCCCGCTTCACGCAGTCGAGTGGCAGACTGCGATCCGACCTGAGAACCGCTTTGTGGGATTGGCTTAACCTCGCGGTTTCGCTGCCC‑TTTGTTCTGTCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCatg                                                                             <  2:1207923/137‑1 (MQ=1)                                                       gagtCCAGCTTCCCGCAGTCGAGTGGCAGACTGCGATCCGAACTCAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTCGCTGCCC‑TTTGTTCTGTCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCatg                                                                             <  2:4842902/137‑1 (MQ=0)                                                       gagtCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTCACCTCGCGGTTTCGCTGCCC‑TTTGTTCTGTCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCatg                                                                             <  2:2324411/137‑1 (MQ=1)                                                       gagtCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTCACCTCGCGGTTTCGCTGCCC‑TTTGTTCTGTCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCatg                                                                             <  1:2318074/137‑1 (MQ=1)                                                                 tcccGCAGTCGAGTTGCAGACTGCTATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTCGCTGCCC‑TTTGTTCTGTCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGAc                                                                   <  2:2420730/137‑1 (MQ=1)                                                                           gatttGCAGACTGCGATACGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTCGCTGCCC‑TTTGTTCTGTCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGTTTTGACGTCATCCCCa                                                         <  2:2677138/137‑1 (MQ=0)                                                                                              gaccTGAGAACAGATTTGTTGGATTGGCTTTACCTCGCGGTTTCGGTGCCC‑TTTGTTCTGTGCATTTTTGCTCGTGTGTAGCCCAGGTCATAAGGGTCATGTTGATTTGTCGTCATCCGATCCTTCCTGCGGTTTGTCAc                                      <  2:1077991/137‑1 (MQ=1)                                                                                                    agtaCAGATTTGTGGGATTGGCTGAACCTCGCGGTTTTGCTGCCC‑TTTGTTCTGTCAATTGTAGCCCGTGTGTAGACCAGGTCATAAGGGTCATGATGATTTGACGTCATCCCCACCTTCCTCCGGTTTGTCACCGGCAg                                <  2:4595631/137‑1 (MQ=1)                                                                                                    agcaCAGATTTGAGGGATGGGCTTAACCTCGCGGGTTCGCTGCCC‑TTTGTTCTGTCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCACCTTCCTCCGGTTTGTCACCGGCAg                                <  2:1690631/137‑1 (MQ=0)                                                                                                         agctttGGGGGATTGGCTTAACCTCGCGGTTTCGCTGCCC‑TTTGTTCTGTCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCACCTTCCTCCGGTTTGTCACCGGCAGTCAcc                           <  1:649231/137‑1 (MQ=1)                                                                                                                   gagtGGCTTAACCTCGCGGTTTCGCTGCCC‑TTTGTTCTGTCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCACCTTCCTCCGGTTTGTCACCGGCAGTCACCTTAGAGTGcc                 <  1:2448862/137‑1 (MQ=1)                                                                                                                   gagtGGCTAAACCTCGCGGTGTCGCTGCCC‑TTTGTTCTGTCCATTGTAGCACGGGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGCCGTCATCCCCACCTTCCTCCGGTTTGTCACCGGCAGTCACCTTAGAGTGcc                 <  1:1295952/137‑1 (MQ=1)                                                                                                                                  gctgTTTCGCTGCCC‑TTTGTTCTGTCAATTGTAGCACGTGTGTAGCACAGGTCATAAGGGGCATGATGTTTTGACGTCATCCCCACCTTCCTCCGGTTTGTCAACGGCAGTCACCTTAGAGTGCCCAACTGAATGCTGGc  >  2:883364/4‑140 (MQ=1)                                                                                                                                         |                                                                                                                                      TGTACAAGGCCCGGGAACGTATTCACCGCGGCATGCTGATCCGCGATTACTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTCGCTGCCC‑TTTGTTCTGTCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCACCTTCCTCCGGTTTGTCACCGGCAGTCACCTTAGAGTGCCCAACTGAATGCTGGC  >  NZ_CP034943/2981117‑2981385

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑