Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NZ_CP034943	2,980,160	G→A	29.2%	noncoding (641/2928 nt)	EO946_RS15430 ←	23S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP034943	2,980,160	0	G	A	29.2%	‑3.8 / 18.0	31	noncoding (641/2928 nt)	EO946_RS15430	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (5/16);  new base A (2/6);  total (7/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GAACTTCACCCTGGACATGGGTAGATCACCTGGTTTCGGGTCTACGACCACGTACTCATGCGCCCTATTCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCGTTAACGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGATCTCTTTCACTCCCCTTCCGGGGTGCTTTTCACCTTTCCCTCAC  >  NZ_CP034943/2980044‑2980291                                                                                                                     |                                                                                                                                    gatcTTCACCCTGGACATGGGTAGATCACGTGGTTTCGGGTCTACGACCACGTACTCATGCGCCCTATTCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCAtt                                                                                                              <  2:4499549/137‑1 (MQ=1)             gggcATGGGTAGATCACCTGGTTTCGGGTCTACGACCACGTACTCATGCGCCCTATTCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCa                                                                                                  <  1:3104984/137‑1 (MQ=1)             ggaccggggTAGATCACCTGGTTTCGGGGCTACGACCACGTACTCATGCGCCCTATTCAGACTCGCTTGCGCTGAGGCTCCGCATCTTCTGCTTAACCTTGCACAGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCa                                                                                                  <  2:555159/134‑1 (MQ=1)                             ccggggttCGGGGCTACGACCACGTACGCAGGCGCCCGATTCCGCCGCGCGGTCGCGGCGGCTCCGCATCGTCTGCGTAACCTTGCACCGGATCGTCACGCGCCGGTTCATTCTACAAAAGGCCCGCCCTCACCCGTTaa                                                                                  <  2:3461723/134‑1 (MQ=0)                                ggattCGGGTCTACGACCACGTACTCATTCGACCTATTCAGTCTCGCTTTCGATGCGGCTCCGCATCTTCTGCTTAACCTAGCACGGGATCGTAACTCGCCGGTTCATTCTTCAAAAGGCACGCCATCACCCGTTAACgg                                                                               >  1:2464071/4‑140 (MQ=1)                                        ctaaacaaCCACATACTAATGCGCCCTATTCATACTCTATTACTCAGCGGCTCCGCATCTTCTGCTTAACCTTGCACAGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCGTTAACGGGCTCTGac                                                                       >  1:962569/8‑140 (MQ=1)                                         tcgaCGACCACGTACTCATGCGCCCTCTTCAGACTCGCGTTCGCTGCGGCTCCGCATCGTCTGCGTAACCGTGCACGGGATCGTAACTCGCCGGTTCATACTACAAAAGGCACGCCATCACCCGTTAACGGGCTCTGact                                                                      <  1:3154655/137‑1 (MQ=1)                                                      acgcATGCGCCCTATTCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACAGGATCGTAACTCGCGGGTTCATTCTACAAAAGGCACGCCATCACCCGTTAACGGGCTCTGACTACTTGTAGGcaca                                                         <  2:4539196/137‑1 (MQ=0)                                                      acgcATGCGCCCTATTCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACAGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCGTTAACGGGCTCTGACTACTTGTAGGcaca                                                         <  2:2962679/137‑1 (MQ=1)                                                      acgcATGCGCCCTAGTCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCGTTAACGGGCTCTGACTACTTGTAGGcaca                                                         <  2:149680/137‑1 (MQ=1)                                                            gctcccTCTTCAGACTAGCTTGAGATGCTGCTCCGCCTCTGCTGCTTAACCTTGCACGGGATCGTAACACGCCGGATCATGCTAAAAAAGGCACGCCATAACACCTAAACGGGCTCTGACGAGTTGTAGCCACACGGtat                                                   >  2:4313015/4‑138 (MQ=0)                                                             cggccTATTCAGACGCGCTTTCGCTGCTGCTCCGCATCTTCTGCTTAACCTTGCAGGGGATCGTAAGTCGCCGGTTCATTCTACATAAGGCACGCCATCACCCGTTAACGGGCTCTGAATTCTTGTAGGCACACGGGTTc                                                  <  2:4272499/137‑1 (MQ=1)                                                                ccgaTTCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCGTTAACGGGCTCTGACTACTTGTAGGCACACGGTTTCAgg                                               <  1:99174/137‑1 (MQ=1)                                                                ccgaTTCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGCTTCATTCTACAAAAGGCACGCCATCACCCGTTAACGGGCTCTGACTACTTGTAGGCACACGGTTTCAgg                                               <  2:1698761/137‑1 (MQ=1)                                                                cccaTTCAGACTCGCTGTCGCAGCGGCACCGCATATTCGGCTGACCCTTGCACGGGATCGTAACTCGCAGTTTGATTCTACAAAAGTCACGCCAGCACCCGTTCACGCGCACGGACTACTTGTAGGCACAAGGTTTCAgg                                               <  2:156090/137‑1 (MQ=0)                                                                  tagtctgACACGCTTGCGCTGCGTGTCCGCATCTTCTGCTTAACCTTGCACACGATCGTAACGCGCCGGTTCATTCTACCAAAGGCACGCCATCACCCGTTAACGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGAt                                             <  2:574952/134‑1 (MQ=1)                                                                  tagtCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCGTTAACGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGAt                                             <  2:934570/137‑1 (MQ=1)                                                                  tagtCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCGTTAACGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGAt                                             <  1:4057232/137‑1 (MQ=1)                                                                  tactCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCTCGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCGTTAACGGGCTCTGACTAGTTGTAGGCACACGGTTTCAGGAt                                             <  1:392354/137‑1 (MQ=0)                                                                  tactCACACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCAGCACCCGTTAACGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGAt                                             <  2:1343307/137‑1 (MQ=0)                                                                          ctggCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTGCATTCTACAAAAGGCACGCCATCACCCGTGAACGGGATCTGACTAATTTTAGGCACACGGTGTCAGGATCTCTTTCa                                     >  2:2380131/4‑140 (MQ=1)                                                                             gccttCGCAGCGGCTCCGCATCGTCTGCGGACCCGTGCACCGGATCGTAACTCGCCGGTTCATTCTACAACAGGCACGCCATCACCCGGTAACGGGCTCTGACTACTGGTAGGCACACGGTTTCAGGAGCTCTTTCACTc                                  <  2:938585/137‑1 (MQ=1)                                                                              ctgtCGCTGCGGCTCCGCATCTTCTGCTTAACCTTTCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCGTTAACGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGATATCTTTCACTcc                                 >  2:883505/4‑140 (MQ=0)                                                                              ctgtCGCTGCGGCTCCGCATCTTATGCTTAACCTTGCACAGGATCGAAACTCGCCGGTTCATTATACAAAAGGCACTAAATAACACGTTAACGGGCTCTGACTACTTGTATGCACACGGTTTCAGGATCTCTATCAATcc                                 >  2:2330612/4‑140 (MQ=1)                                                                                       cgccTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCGTTAACGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGATCTCTTTCACTCCCCTTCCggg                        >  2:2858786/4‑140 (MQ=1)                                                                                        ggatCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCGTTAACGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGATCTCTTTCACTCCCCTTCCgggg                       <  2:677177/137‑1 (MQ=1)                                                                                             cgaaTCTTCTGCTTAACCTTGCACAGGATCGTAACTCGCCGGTTCATTTTACAAAAGGCAGGCAACCACCCGTTAACGGGCTCTTACTACTTGTAGGCACACGTTTTCAGCAGATCTTTCCCCCCCCTTCCGGGGTGCtt                  <  2:560164/137‑1 (MQ=1)                                                                                                            acctgtcACGGGATCGTCACTCGCCGGTTAATTCTACAAAAGGAACGCCATCACGAGTTAACGGGCTATGACTACTTGTAGTCACTCGGTTTCAGGTACTCTTTCACTGGCAGTCAGGGGTGCTTTTCACCTTTCGCTCa   <  2:4042528/134‑1 (MQ=0)                                                                                                            acctgccACGGGCTCGGAACTCCCCGGTTCATTCAACAAAAGGCCCGCCATCACCCGTTAACGGGCTCTGACGACTTGTAGGCACACGGTTTCAGGATCTCTGTCACTCCCCTTCCGGGGTGCTTTTCACCTGTCCCTCa   <  2:2645639/134‑1 (MQ=1)                                                                                                             ccgtGCACGGGATCCTAACTCGCCGGTTCATTCTACAAAAGGCACGCCCTCACCCGGTAACCGGCTCGGACGACGTGTAGGCCCACGGTTTCAGGATCACTTTCACTCCCCTTCCGGGGTGCTTTTCACCTTTCCCTCAc  <  1:4039954/137‑1 (MQ=1)                                                                                                             ccgtGCACAGGATCGTAACTCGCCGGTTCAGTCTACACAAGCCACGCCATCACCCGGTCACGGGCTCTGACTACTTGTAGGCACCCGGTTTCAGGATCTCTTTCACTCCCCTTCCGGGGTGCTTTTCACCTGTCCCTCAc  <  2:3023629/137‑1 (MQ=1)                                                                                                                     |                                                                                                                                    GAACTTCACCCTGGACATGGGTAGATCACCTGGTTTCGGGTCTACGACCACGTACTCATGCGCCCTATTCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCGTTAACGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGATCTCTTTCACTCCCCTTCCGGGGTGCTTTTCACCTTTCCCTCAC  >  NZ_CP034943/2980044‑2980291

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