Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NZ_CP034943	2,981,239	T→A	13.1%	noncoding (1281/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,239	0	T	A	13.1%	2.7 / 8.2	56	noncoding (1281/1550 nt)	EO946_RS15435	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (9/39);  new base A (3/4);  total (13/43)
Fisher's exact test for biased strand distribution p-value = 1.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGACGGGCGGTGTGTACAAGGCCCGGGAACGTATTCACCGCGGCATGCTGATCCGCGATTACTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTCGCTGCCCTTTGTTCTGTCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCACCTTCCTCCGGTTTGTCACCGGCAGTCACCTTAGAGTGCCCA  >  NZ_CP034943/2981105‑2981372                                                                                                                                       |                                                                                                                                      tggcGGGCGGTGTGAACAAGGCCCGGGAACGTATTCACCGCGGCATGCTGATCCGCGATTACTAGCGATTCCAGCTTCACGCAGTCGCGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCt                                                                                                                                  >  2:7184378/4‑140 (MQ=0)      gggggTGTGTACAAGGCCCGGGAACGTATTCACCGCGGCATGCTGATCCGCGATTACTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTAAACCTCGCgg                                                                                                                             <  2:1980001/137‑1 (MQ=1)        cgctgtgtACAAGGCACTGGAAATAAAACACCGCGGAATGCTGATCCGCGATTACTAGCGATTCCAGCTTCACGCAGTCGAGAAGCAGAATACGATCCTAACTAAGAACATATTAGTGGGATAGGCTAAACATCGCGGtt                                                                                                                           >  2:3653398/4‑140 (MQ=0)         ggggtgtAGAAGGCCCGGGACCGTATTCCCCGCGGCTTGCTGCTCCGCGATTACGAGCGATTCCAGCTTCACGCAGTGGTGTTGCTGACTGCGATCCGAACTGAGACCAGATTGGTGGGCTTGGCTCAACCTCGCGGtgt                                                                                                                          >  2:4003387/4‑138 (MQ=1)           tgggtAGAAGGCCCGGGAACGTATTCAGCGCGGCATGCTGATCCGCGATTACTAGCGATTCCAGCGTCACGCAGTCGAGTTGCAGACTGCGCTACGAACTGAGAACAGATGTGTGGGAGTGGCTTAACCTCGCGGttttg                                                                                                                        >  2:432013/4‑138 (MQ=1)             tgtcccaGGCCCGGTAACTTATTCACAGCGGCATGCTGATCCGCGATTACTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATACGAACTGAGAACCGATTTGTGGGATTGGCTTAACCTCGCGGTTTTGCt                                                                                                                      <  1:3297261/134‑1 (MQ=1)             tggaCAAGGCCCGGGAACGTATTCACCGCGGCCTGCGGATCCGCGATTACTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGGGGGATTGGCTTAACCTCGCGGTTTTGCt                                                                                                                      <  2:4795833/137‑1 (MQ=1)             tggaCAAGGCCCGGGAACGTATTCACCGCGGCATGCTGATCCGCGATTACTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGCCTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTCGCt                                                                                                                      <  2:5681235/137‑1 (MQ=1)             tggaCAAGGCCCGGGAACGTATTCACCGCGGCATGCTGATCCGCGATTACTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTCGCt                                                                                                                      <  1:4433002/137‑1 (MQ=1)             tggaCAAGGCCCGGGAACGTATTCACCGCGGCATGCTGATCCGCGAGTACGAGAGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTCGCt                                                                                                                      <  1:4757923/137‑1 (MQ=0)             tggaCAACGCCCGGGAACGTATTCACCGCGGCATGCTGATCCGCGATTACTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTCGCt                                                                                                                      <  1:6478204/137‑1 (MQ=1)             ggtacccggCCCGGGACCGTCTTGTCGGCGGCAGGCTGATCCGCGTTGACTAGCGATTCCAGCTTCACGCTGTCGCGTGGCAGTCTGCGCGCCGGTCTGAGAACAGCTTTGAGGGGTTTGCTTAACCTCGCGGTTTTGCt                                                                                                                      <  1:5899733/133‑1 (MQ=0)                    ggacccggAACGTATTCACCGCGGAATGATGATCCGCGATTACTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTCGCTGCCCttt                                                                                                               >  2:941552/7‑140 (MQ=0)                    ggacccggAACGTATTCAACGCGGCATAATAAAAAACGATTAATATAGATTCCAGCTACACGCAGTCGAGTTGCAGACTGCGATACGAACTCAGAACAGATATGTGGCAATGGCTTAACCTCTCGGTTTCGCTGCCCttt                                                                                                               >  2:6014742/7‑140 (MQ=0)                     gacccggAACGTATTCACCGCGTAATAAAAAACCGCGATTACTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGAGGGATTGGCTTAACCTCGCGGTTTTGCTGCCCTTTg                                                                                                              >  2:1579817/6‑140 (MQ=1)                       cccggAACGTATTCACCGCGGCATGCTGATCCGCGATTACTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTCGCTGCCCTTTGtt                                                                                                            >  2:2303666/4‑140 (MQ=1)                         ggtaaCGTATTCAGCGCGGCATGCTGATCCGCGATTACTAGCGATTCCAGCGTCACGCAGTCGAGTTGCAGACTGCGATCCGACCTGAGAACAGATTTGTGGGAGTGGCAAAACCTCGCGGTTTCGCTGCCCTTTGTTCt                                                                                                          >  2:3815396/4‑140 (MQ=1)                         ggtaaCGTATTCACCGCGGCATGCTGATCCGCGATGACTCGCGATTCCAGCTTCACGCAGTCGCGTTGCAGACTGCGATCCGACCTGAGAACCGAGTTGTGGGATTGGCTTAACCTCGCGGTTTTGCTGCCCTTTGTTCt                                                                                                          <  1:7273832/137‑1 (MQ=1)                           gaccGTATTCACCGCGGCATGCTGATCCGCGATTCCTCGCGCTTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGAGTTGTGGGATTGGCTTACCCTCGCGGTTTTGCTGCCCTTTGTTCTGt                                                                                                        <  2:2698530/137‑1 (MQ=1)                            aaagTATTCACCGCGGCATGCTTATCCGCGATTACTAGCGATTCCATCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTTAGAAAATATTTGTGGGATTGGCTTAACCTCGCGGTTTCGCTGCCCTTTGTTCTGTc                                                                                                       <  2:6811741/137‑1 (MQ=1)                               gtcttCACCGCGGCATGCTGATCCGCGATTACGAGCGAGTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTCGCTGCCCTTTGTTCTGTCCAt                                                                                                    <  2:4177758/137‑1 (MQ=0)                                       cgggggaTGCTGATCAGCGATTACTAGCGATTCCAGCTTCCCGCAGTCGAGGTGCAGACTGCGATTCGAACTGAGAACAGCTTTGTGGGATTGGCTTAACCTCGCGGTTTCGCTGCCCTTTGTTCTGTGCATTGTTGCAc                                                                                            <  1:3029056/134‑1 (MQ=1)                                       cggggCATGCTGATCCGCGATTAATTGCGATTGCAGCGTCACGGAGTCGAGTTGCAGACTTCGATCCGAACGGTCAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTCGCTGCCCTTTGTTCTGTCCATTGTAGCAc                                                                                            <  1:3982710/137‑1 (MQ=1)                                       cgaggCATGCTGATCCGCGATTACTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACAGAGAACAGATTTGTGGGATTGGCTAAACCTCGCGGTTTCGCTGCCCTTTGTTCTGTCCAGTGTAGCAc                                                                                            <  2:2079967/137‑1 (MQ=0)                                            caggCTGATCCGCGCTTACTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGCTCCGAACTGAGAACAGAGTTGTGGGATTGGCTTAACCTCGCGGTTTTGCTGCCCTTTGTTCTGTCCATGGTAGCACgtgtg                                                                                       <  1:2818076/137‑1 (MQ=1)                                            caggCTGATCCGCGATTACTAGCGATACCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTTGCTGCCCTTTGTTCTGTCCATTGTAGCACgtgtg                                                                                       <  2:3827400/137‑1 (MQ=0)                                            caagCTGATCCGCGATTACTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTTGCTGCCCTTTGTTCTGTCCATTGTAGCACgtgtg                                                                                       <  2:6850784/137‑1 (MQ=1)                                                 tgctCCGCGATTACTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTCGCTGCCCTTTGTTCTGTCCATTGTAGCACGTGTGTAGcc                                                                                  <  2:6998210/137‑1 (MQ=1)                                                 tgctCCGCGATTACTACCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTACCCTCGCGGTTTCGCTGCCCTTTGTTCTGTCCATTGTAGCACGTGTGTAGcc                                                                                  <  2:7711636/137‑1 (MQ=0)                                                 tgctCCGCGATTACGAGCGAGTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTCACCTCGCGGTTTCGCTGCCCTTTGTTCTGTCCATTGTAGCACGTGTGTAGcc                                                                                  <  2:1678746/137‑1 (MQ=0)                                                  gagccGCGATTCCTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTCACCTCGCGGTTTCGCTGCCCTTTGTTCTGTCCATTGTAGCACGTGTGTAGccc                                                                                 <  1:2301651/137‑1 (MQ=0)                                                  gagccGCGATTACTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTTGCTGCCCTTTGTTCTGTCCATTGTAGCACGTGTGTAGccc                                                                                 <  1:1946975/137‑1 (MQ=1)                                                  gagccGCGATTACTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTAAACCTCGCGGTTTCGCTGCCCTTTGTTCTGTCCATTGTAGCACGTGTGTAGccc                                                                                 <  2:5712606/137‑1 (MQ=1)                                                  gagccGCGAATACTCGCGAAGCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGAGTCGCTTAACCTCGCGGTTTCGCTGCCCTTTGTTCTGTCCATTGTAGCACGTGTGTAGccc                                                                                 <  2:4132892/137‑1 (MQ=1)                                                  gaaccGCGATTACTAGCGATTCCCGCTACACGCAGTCGAGTTGCAGACAGCGATCCGAACAGAGAACAGATTTGTGGGATTGGCTTAACCACCCGGTTTCGCTGCCCTTTGTTCTGTCCATTGTAGCACGTGTGTAGccc                                                                                 <  2:7564873/137‑1 (MQ=1)                                                  gaaccGCGATTACTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTCACCTCGCGGTTTCGCTGCCCTTTGTTCTGTCCATTGTAGCACGTGTGTAGccc                                                                                 <  2:7757878/137‑1 (MQ=1)                                                          attaaaaGCGATTCCAGCATCACGAAGTCCATTTGCAGACTCAGATCCGAACTGAGAACAGATAAGTCGGATTGGCATCACCTCCCGTTAACGCTGACCTTTGTTCAGTCCATTGTAGAACTTGAGTATCCCATGTCATa                                                                         >  1:7587349/7‑140 (MQ=0)                                                          atgaCTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGAGTTGTGGGATTGGCTTAACCTCGCGGTTTTGCTGCCCTTTGTTCTGTCCATTGTAGCACGTGTGTAGCCCAGGTCATa                                                                         <  1:7713436/137‑1 (MQ=0)                                                           ttccTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTCGCTGCCCTTTGTTCTGTCCATTGTAGCACGTGTGTAGCCCAGGTCATaa                                                                        <  2:4665544/137‑1 (MQ=1)                                                                   gagtCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTTGCTGCCCTTTGTTCTGTCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCatg                                                                <  1:6568440/137‑1 (MQ=1)                                                                   gagtCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTAAACCTCGCGGTTTCGCTGCCCTTTGTTCTGTCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCatg                                                                <  2:7009926/137‑1 (MQ=1)                                                                   gagtCCAGCTTCACGCAGTCGACTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTCGCGGCCCTTTGTTCTGTCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCatg                                                                <  1:4236331/137‑1 (MQ=0)                                                                             tcccGCTGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTTGCTGCCCTTTGTTCTGTCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGAc                                                      <  1:3357260/137‑1 (MQ=0)                                                                                        aggtGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTCGCTGCCCTTTGTTCTGTCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCAc                                           <  1:7108320/137‑1 (MQ=1)                                                                                                  tgtgATCCGAACTGAGAACCGATTTGTGGGATTGGCTTAACCTCGCGGTTTCGCTGCCCTTTGTTCTGTCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCACCTGCCTCCgg                                 >  2:6822870/4‑140 (MQ=0)                                                                                                          gcacggagaACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTTGCTGCCCTTTGTTCTGTCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCACCTTCCTCCGGTTTGTCAc                         <  2:3232004/135‑1 (MQ=0)                                                                                                          gaccTGAGAACAGATTTGTGGGATTGGCTTCACCTCGCGGTTTCGCTGCCCTTTGTTCTGTCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCACCTTCCTCCGGTTTGTCAc                         <  1:5890761/137‑1 (MQ=1)                                                                                                                aggaCAGATTTGTGGGATTGGCTTAACCTCGCGGTTTTGCTGCCCTTTGTTCTGTCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCACCTTCCTCCGGTTTGTCACCGGCAg                   <  1:1526473/137‑1 (MQ=1)                                                                                                                agcaCAGATTTGTGGGATTGGGTTAACGTCGCGGTTGTGCTGCACTTTGTTCTGTCCTTTGTAGCTCGGGTTTAGCCCAGGTCATAAGGGGCATGTTGATTTGACGTCATGCGCACGTTCCTCCGGTTTGTCACCGGctg                   <  2:5043204/137‑3 (MQ=1)                                                                                                                agcaCAGATTTGTGGGATTGGCTTAACCTCGCGGTTTCGCTGCCCTTTGTTCGGTCCATTGTAGCACGTTTGTAGCCCAGGTCATAAGGTGCATGATGATTTGACGTCATCCCCACCGTCCTCCGGTTTGTCACCGGCAg                   <  2:7668463/137‑1 (MQ=1)                                                                                                                agcaCAGATTTGTGGGATGGGCTTAACCTCGCGGTTTTGCTGCCCTTTGTTCTGTCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCACCTTCCTCCGGTTTGTCACCGGCAg                   <  2:3114984/137‑1 (MQ=0)                                                                                                                 gaccAGATTTGTGGGATTGGCTTAACCTCGCGGTTTCGCTGCCCTTTGTTCTGTCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCACCTTCCTCCGGTTTGTCACCGGCAGt                  <  2:6363961/137‑1 (MQ=1)                                                                                                                          ggtcggAGTGGCTTAACCTCGCGGTTTCGCTGCCGTTTGGTCTGTCCAGTGTAGCACGTGGGGAGCCCAGGTCATAAGGGGCATGATGATTTGACGGCCTCCCCGCCTTCCGCCGGTTTGTCACCGGCAGTCCCCTgaga         >  2:1575461/5‑140 (MQ=1)                                                                                                                             tggtttGTCTTATGCTCGCGGTTTCGCTGCACATTGTTCTGTCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCGCACCGTCCTCCGGTTTGTCACCGGCAGTCACCTTAGAGTg      <  1:4520717/136‑1 (MQ=1)                                                                                                                               gagtGGCTTAACCTCGCGGTTTCGCTGCCCTTTGTTCTGTCCAGTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCACCTTCCTCCGGTTTGTCACCGGCAGTCACCTTAGAGTGcc    <  2:855424/137‑1 (MQ=1)                                                                                                                                 ttcgCTAAACCTCGCGGGATCGCTGCCCTTTGTTCTGTCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCACCTTCCTCCGGTTTGTCACCGGCAGACACCTTAGACTGCCCa  >  2:5748777/4‑140 (MQ=1)                                                                                                                                       |                                                                                                                                      TGACGGGCGGTGTGTACAAGGCCCGGGAACGTATTCACCGCGGCATGCTGATCCGCGATTACTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTCGCTGCCCTTTGTTCTGTCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCACCTTCCTCCGGTTTGTCACCGGCAGTCACCTTAGAGTGCCCA  >  NZ_CP034943/2981105‑2981372

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