Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NZ_CP034943	2,980,263:1	+C	11.3%	noncoding (538/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,263	1	.	C	11.3%	4.0 / 27.9	53	noncoding (538/2928 nt)	EO946_RS15430	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base . (3/44);  new base C (2/4);  total (5/48)
Fisher's exact test for biased strand distribution p-value = 9.25e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.46e-01

GCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCG‑TTAA‑CGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGATCTCTTTCACTCCCCT‑TCCGGGGTGCTTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATC  >  NZ_CP034943/2980146‑2980390                                                                                                                         |                                                                                                                                gcgtAACCGTGCACGGGATCGTAACTCGCCGGTTCAGTCTACAAAAGGCACGCCATCACCCG‑TTAA‑CGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGATCTCTTTCACTCCCCT‑TCCGGGGTGCTTTTCACCTTTc                                                                                                           <  2:5435919/137‑1 (MQ=0)  ctgaaCCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCG‑TTAA‑CGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGATCTCTTTCACTCCCCT‑TCCGGGGTGCTTTTCACCTTTcc                                                                                                          <  1:4517959/137‑1 (MQ=1)    gaaccgggCACGGGCTCGCACCCCGCCGCTTGAGTCTACAAAACGCCCGCCATCCCCCG‑TTAA‑CGGGCTCGGACTCCTGGTATGCACACGGGTTCAGGATCTCTTGCACTCCCCT‑TCCGGGGTGCTTTTCACCGGTCCCt                                                                                                        <  2:2941715/133‑1 (MQ=0)       ccgtGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCG‑TTAA‑CGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGATCTCTTTCACTCCCCT‑TCCGGGGTGCTTTTCACCTTTCCCTCAc                                                                                                     <  2:2113600/137‑1 (MQ=1)       ccgtGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCG‑TTAA‑CGGGCTCTGACTACTTGTAGGCACACGGGTTCAGGATCTCTTTCACTCCCCT‑TCCGGGGTGCTTTTCACCTTTCCCTCAc                                                                                                     <  1:1576572/137‑1 (MQ=1)       ccctGCACGGGATCGTAACTCTCCGGTTCATTCTAGAAAAGGCACGCCATCACCCG‑TTAA‑CGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGATCTCTTTCACTCCCCT‑TCCGGGGTGCTTTTCACCTTTCCCTCAc                                                                                                     <  2:2280523/137‑1 (MQ=0)          tgaaCGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCG‑TTAA‑CGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGATCTCTTTCACTCCCCT‑TCCGGGGTGCTTTTCACCTTTCCCTCACGGt                                                                                                  >  2:3347368/4‑140 (MQ=1)            ccccggctCGTACCTCGCCGCTCCCTCCCACCCACGGCACGCCAGCACCCC‑CTAA‑CGGGCGCGGACTACTGGTAGGCACACGGGGTCCGGATCTCTTTCACTCCCCG‑TCCGGGGTGCTTGTCACCTTTCCCTCACGGTAc                                                                                                <  2:535284/133‑1 (MQ=0)            catgggATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCG‑TTAA‑CGGGCTCTGACTACTGGTAGGCACACGGTTTCAGGATCTCTTTCACTCCCCT‑TCCGGGGTGCTTTTCACCTGTCCCTCACGGTAc                                                                                                <  2:2461933/137‑1 (MQ=0)                        acgcGCCGGTTCATTCTACAAAAGGCACGCCATCACCCG‑TTAA‑CGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGATCTCTTTCACTCCCCT‑TCCGGGGTGCTTTTCACCTTTCCCTCACGGTACTGGTTCACTATc                                                                                    <  1:812291/137‑1 (MQ=1)                           cggcGGTTCATTCTACAAAAGGCACGCCATCACCCG‑TTAA‑CGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGATCTCTTTCACTCCCCT‑TCCGGGGTGCTTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGt                                                                                 <  1:2937219/137‑1 (MQ=1)                                   cagtCTACAAAAGGCACGCCATCACGCG‑TTAA‑CGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGATCTCTTTCACTCCCCT‑TCCGGGGTGCTTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCTCTAggg                                                                         <  2:3360999/137‑1 (MQ=0)                                   cagtCTACAAAAGGCACGCCATCACCCG‑TTAA‑CGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGATCTCTTTCACTCCCCT‑TCCGGGGTGCTTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACTAggg                                                                         <  2:5527893/137‑1 (MQ=1)                                     ttatACAAAAGGTACGCAGTCACACGCCTAAGCGTGCTCCCACTGCTTGTACGTACACGGTTTCAGGTTCTTTTTCACTCCCCTCGCCGGGGTTCTTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACTAgg                                                                          <  2:1151961/137‑1 (MQ=0)                                     ttatACAAAAGGTACGCAGTCACACGCCTAAGCGTGCTCCCACTGCTTGTACGTACACGGTTTCAGGTTCTTTTTCACTCCCCTCGCCGGGGTTCTTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACTAgg                                                                          <  2:2485120/137‑1 (MQ=0)                                     ttatACAAAAGGTACGCAGTCACACGCCTAAGCGTGCTCCCACTGCTTGTACGTACACGGTTTCAGGTTCTTTTTCACTCCCCTCGCCGGGGTTCTTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACTAgg                                                                          <  2:432902/137‑1 (MQ=0)                                     ttatACAAAAGGTACGCAGTCACACGCCTAAGCGTGCTCCCACTGCTTGTACGTACACGGTTTCAGGTTCTTTTTCACTCCCCTCGCCGGGGTTCTTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACTAgg                                                                          <  2:2181517/137‑1 (MQ=0)                                         agtaatggCAAGACATCACCCC‑TGAA‑CGGGCTCGGTCTACTTGTAGGCACACGCTTTCAGGATCTCTTTCACTCGCCT‑TCCGTGTTGCTTTTCACCTTTCCCTCACGGTACTGGTTCACTATCCGTCACTAGGGAGTAtt                                                                   <  2:4526634/134‑1 (MQ=1)                                          ccacagccACGCCATCACCCG‑TTAA‑CGGGCTCTGACTACGTGTAGGCACACGGTTTCAGGATCTCTTTCACTCCCCT‑TCCGGGGTGCTTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCCCTAGGGAGTAttt                                                                  <  1:338107/133‑1 (MQ=0)                                          cacaaGGCACGCCATCACCCG‑TTAC‑CGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGATCTCTTTCACTCCCCG‑TCCGGGGTGCTTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACTAGGGAGTAttt                                                                  <  1:5092325/137‑1 (MQ=0)                                          cacaaGGCACGCCATCACCCG‑TTAA‑CGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGATCTCTTTCACTCCCCT‑TCCGGGGTGCTTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACTAGGGAGTAttt                                                                  <  1:5232586/137‑1 (MQ=1)                                          cacaaGGCACGCCATCACCCG‑GTCA‑CGGGCTCTGACTACCGGTAGGCACACGGGGTCAGGATCGCTTTCCCTCCCCT‑GCCGGGGTGCTTGTCACCTTTCCCTCCCGGTACTGGGTCACTATCGGTCACTAGGGAGTAttt                                                                  <  2:5638803/137‑1 (MQ=1)                                           aacaggcaggcCATCACCCG‑TTAA‑CGGGCTCTGACTAATTGTATGCACACGTTTTCATGATCTCTTTCACTGCACT‑TCCGGGTTGCTTTTCTCCTTTCACTCAATGTACTGGTTCACTATCGGTCATTTTGGAGTAtttt                                                                 <  2:2077190/137‑2 (MQ=1)                                              agtcACGCCATCACCCG‑TTAA‑CGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGATCTCTTTCACTCCCCT‑TCCGGGGTGCTTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACTAGGGAGTATTTAGcc                                                              <  1:262488/137‑1 (MQ=1)                                                gcccGCCATCACCCG‑TGAA‑CGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGATCTCTTTCACTCCCCT‑TCCGGGGTGCTTTTCACCTGTCCCTCACGGTACTGGTTCACTATCGGTCACTAGGGAGTATTTAGCCtt                                                            <  2:4451644/137‑1 (MQ=0)                                                            cccgggaa‑CGGGCCCTGACTACACGTAGGCACACGGCCTCAGCAACTCCTCCCCTCCCCG‑ACCGCGCTGCAATTCACCTGACCCCCCCCGTACCGGTTCACTATCGGTCACTAGGGAGTATTAAGCCTAGCGAGATGGtc                                                 <  2:4266759/134‑1 (MQ=0)                                                                            tgtctacgTGTAGGCACACGTTTTCATTATCTCTTTCACTCCCCT‑TCCGGGGTGCTTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGAc                                  <  1:4889940/137‑1 (MQ=1)                                                                               ctgctTGTACGTACACGGTTTCAGGTTCTTTTTCACTCCCCTCGCCGGGGTTCTTTTCGCCTTTCCCTCACGGTACTGGTTCACTATCGGTCAGTCAGGAGTATTTAGCCTTGGAGGATGGTCCCCCCATATTCAGAcag                                >  1:793273/4‑138 (MQ=0)                                                                               ctgctTGTACGTACACGGTTTCAGGTTCTTTTTCACTCCCCTCGCCGGGGTTCTTTTCGCCTTTCCCTCACGGTACTGGTTCACTATCGGTCAGTCAGGAGTATTTAGCCTTGGAGGATGGTCCCCCCATATTCAGAcag                                >  1:3687113/4‑138 (MQ=0)                                                                                         gaactcGGTTTCAGGATCTCTTTCACTGCCCT‑TCCGTGGTGCTTTTGAGCTTTCCCTCCCGGTACTTGTTAACTATCGGTCCGTAGGGAGTATTTTGCCTTGGGTGATGTTCCTCACGTATTCAGACGGAATTTCACgtg                     <  1:2671576/135‑1 (MQ=1)                                                                                              cgctttCAGGATCTCTTTCACTCCCCT‑TCCGTGGTTCTTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTccg                <  1:5277467/137‑1 (MQ=0)                                                                                                gtgtCAGGATCTCTTTCACTCCCCT‑TCCGGGGTGCTTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGGGTTccgcc              <  2:1050339/137‑1 (MQ=1)                                                                                                gtgtCAGGATCTCTTTCACTCCCCT‑TCCGGGGTGCTTTTCACCTTGCCCTCACGGTACTGGTTCACTATCGGTCACTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTccgcc              <  1:4814224/137‑1 (MQ=1)                                                                                                         tcactatCACTCCCCT‑TAAGGGCTGCTTTTCACATATCCCTCACGGAAATGGTACACTATAGGACACTCGGGAGTATTTAGACTTGGGAGATGGTACTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAgc     >  2:2054883/7‑139 (MQ=1)                                                                                                           tcgttCACTCCCCT‑TCCGGGGTGCTTTTCACCTTTCCCTCCCGGTACTGGTTCACGATCGGTCACTAGGGAGTATTTAGGCTTGGGAGATGGTACTCCCGGATTCCGACGGAATTTCAAGTGTTCCGCCGTACTCAGGAt   <  1:3150491/137‑1 (MQ=1)                                                                                                           tcgttCACTCCCCT‑TCCGGGGTGCTTTTCACCTTTCCCTCACGGTTCTGGTTCACTATCGGTCACTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGTATTTCCCGTGTTCCGCCGTACTCAGGAt   <  1:2435142/137‑1 (MQ=0)                                                                                                           tcgttCACTCCCCT‑TCCGGGGTGCTTTTCACCTTTCCCTCACGGTCCTGGTTCACTATCGGTCACTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGAt   <  2:2638311/137‑1 (MQ=1)                                                                                                           tcgttCACTCCCCT‑TCCGGGGTGCTTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGAt   <  2:2834512/137‑1 (MQ=1)                                                                                                           tcgttCACTCCCCT‑TCCGGGGTGCTTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGAt   <  2:3405327/137‑1 (MQ=1)                                                                                                           tcgttCACTCCCCT‑TCCGGGGTGCTTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGAt   <  2:4723850/137‑1 (MQ=1)                                                                                                           tcgttCACTCCCCT‑TCCGGGGTGCTTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGAt   <  1:5081934/137‑1 (MQ=1)                                                                                                           tcgttCACTCCCCT‑TCCGGGGTGCTTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGAt   <  1:4998985/137‑1 (MQ=1)                                                                                                           tcgttCACTCCCCT‑GCCGGGGTGCTTTTCACCTGGCCCACCCGGGACTGGTTCACTATCGGTCACTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGAt   <  1:5548619/137‑1 (MQ=1)                                                                                                           tcgttCACTCCCCG‑TCCGGGGTGCGTTTCCCCTGTCCCGCCCGGTACTGGGTCACTCTCGGTCACGAGGGAGTATTTCGCCGTGGGAGATGGGCCGCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGag   <  2:861735/137‑2 (MQ=1)                                                                                                           tcgttCACGCCCCT‑TCCGGGGTGCTTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACTAGGGAGTATTGAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGAt   <  1:3611859/137‑1 (MQ=0)                                                                                                           tccttCACTCCCCT‑TCCGGGGTGCTTGTCACCTTTCCCTCCCGGTACTGGTTCACTATCGGTCACTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGAGTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGAt   <  1:1603150/137‑1 (MQ=0)                                                                                                           tcattCACTCCCCT‑TCCGGGGTGCTTTTCACCTTTCCCTCCCGGTACTGGTTCACTATCGGTCACTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGACTTACACGTGTGACGCCGGACTCAGGag   >  2:3386150/4‑139 (MQ=1)                                                                                                           tcattCACTCCCCT‑TCCGGGGTGCTTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGAt   <  2:3250385/137‑1 (MQ=1)                                                                                                           gctgtcccTCCCCG‑TCCGGGGCGCTTTCCACCTGGCCCTCACGGTACTCGGTCACTAGCGGGCACTAGGGAGTATTTAGCCGTGGGAGATGGTCCTCCCGGATTCCGACGGAATTACACGTGTTCCGCCGTACTCAGGAt   <  2:5509059/133‑1 (MQ=1)                                                                                                           gctgtCACTCCCCT‑TCCGGGGTGCTTTTCACCTTTCCCTCACGGTAGTGGGTCACTATCGGTCACTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGAt   <  1:4789433/136‑1 (MQ=0)                                                                                                           gctgtCACGCCCCT‑TCCGGGGTGCTTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACTAGGGAGTATTGAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGAt   <  1:1168960/136‑1 (MQ=0)                                                                                                           ccttgcccTCCCCT‑GCCGGGGTGCGTTTCACCTTGCCCTCACGGTACTGGTTCACTATCGGTCACTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGAt   <  1:3318331/133‑1 (MQ=1)                                                                                                            ctgtCACGCCCCT‑TCCGGGGTGCTTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATc  <  1:5789845/137‑1 (MQ=1)                                                                                                                         |                                                                                                                                GCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCG‑TTAA‑CGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGATCTCTTTCACTCCCCT‑TCCGGGGTGCTTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATC  >  NZ_CP034943/2980146‑2980390

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