Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NZ_CP034943	2,980,231	G→C	14.3%	noncoding (570/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,231	0	G	C	14.3%	0.0 / 12.2	42	noncoding (570/2928 nt)	EO946_RS15430	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (6/28);  new base C (2/4);  total (8/34)
Fisher's exact test for biased strand distribution p-value = 5.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

CGTACTCATGCGCCCTATTCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCG‑TTAA‑CGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGATCTCTTTCACTCCCCT‑TCCGGGGTGCTTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGG  >  NZ_CP034943/2980094‑2980360                                                                                                                                            |                                                                                                                                   cgaaCTCATGCGCCCTATTCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCAAGCCATCACCCG‑TTAA‑CGGGCACTGACTACTTGTAGGc                                                                                                                                  >  2:1645250/4‑140 (MQ=0)    ccgcaggcgcCCTAGTCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCG‑TTAA‑CGGGCTCTGACTACTTGTAGGcaca                                                                                                                               <  2:3481786/134‑1 (MQ=0)    acgcATGCGCCCTAGTCAGACTCGCTTTCGCTGCGGCTCCGCATCGTCTGCTGAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACCAAAGGCACGCCATCACCCG‑TTCA‑CGGGCTCGGACTACTTGTAGGcaca                                                                                                                               <  1:291971/137‑1 (MQ=1)        attcgcCCTATTCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCG‑TTAA‑CGGGCTCTGACTACTTGTCGGCACACGGt                                                                                                                           >  2:5589187/4‑140 (MQ=1)        attcgcCCTATTCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGACATCACCCG‑TTAA‑CGGGCTCTGACTACTGGTAGGCACACGGt                                                                                                                           >  2:576254/4‑140 (MQ=0)                tagtCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACAGGATCGTAACTCGCCGGTTCATTCTACAACAGGCACGCCATCACCCG‑TTAA‑CGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGAt                                                                                                                   <  2:558178/137‑1 (MQ=0)                           gcgttCGCTGCTGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCG‑TTAA‑CGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGATCTCTTTCACTc                                                                                                        <  2:3355635/137‑1 (MQ=1)                           gcgttCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCG‑TTAA‑CGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGATCTCTTTCACTc                                                                                                        <  1:2477068/137‑1 (MQ=1)                           gcgttCGCTGCGGCTCCGCATCGTCTGCTTAACCGTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCG‑TTAA‑CGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGATCTCTTTCACTc                                                                                                        <  1:4841356/137‑1 (MQ=0)                           gccttCGCTGCGGCGCCGCATCTACTGCTGAACCTGGCACGGGATCGTAACACGCCGGTTCATTCTACAAAAGGCCCGCCATCACCCG‑TTAA‑CGGGCTCTGACTACTTGTAGGCACACGGTTTCACGATCTCTTACACTc                                                                                                        <  2:4331962/137‑1 (MQ=1)                           gccttCGCGGCGGCTCCGCATCTTCTGCATAACCTGGCACGGGAGCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCCCCCG‑TGAA‑CGGGCTCTGACTACATGTAGGCACACGGTTTCAGGATCGCTTTCACTc                                                                                                        <  2:866941/137‑1 (MQ=1)                            ctgtCGCTGCGGCTCCGCGTCTTCTGCTTAACCTATCACAGGATGGTAACTCGCCGGGTCATTCTACAAAAGGAACGCCATCACCCG‑TTAA‑CGGGCTCTGACTACTTGTAGGCCCACGGGTTCAGGAGCTGTTTCCCTcc                                                                                                       >  2:1077188/4‑140 (MQ=1)                            ctgtCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCG‑TTAA‑CGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGATCTCGTTCACTcc                                                                                                       >  2:697576/4‑140 (MQ=1)                                                     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:3687113/4‑138 (MQ=0)                                                                                                                                   ctgctTGTACGTACACGGTTTCAGGTTCTTTTTCACTCCCCTCGCCGGGGTTCTTTTCGCCTTTCCCTCACGGTACTGGTTCACTATCGGTCAGTCAGGAGTATTTAGCCTTGGAGGATGGTCCCCCCATATTCAGAcag  >  1:793273/4‑138 (MQ=0)                                                                                                                                            |                                                                                                                                   CGTACTCATGCGCCCTATTCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCG‑TTAA‑CGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGATCTCTTTCACTCCCCT‑TCCGGGGTGCTTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGG  >  NZ_CP034943/2980094‑2980360

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