Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NZ_CP034943	2,980,247	A→T	14.7%	noncoding (554/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,247	0	A	T	14.7%	0.8 / 11.8	41	noncoding (554/2928 nt)	EO946_RS15430	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (3/32);  new base T (2/4);  total (5/36)
Fisher's exact test for biased strand distribution p-value = 1.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

TATTCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCG‑TTAA‑CGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGATCTCTTTCACTCCCCT‑TCCGGGGTGCTTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCC  >  NZ_CP034943/2980109‑2980378                                                                                                                                             |                                                                                                                                     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:432902/137‑1 (MQ=0)                                                                          ttatACAAAAGGTACGCAGTCACACGCCTAAGCGTGCTCCCACTGCTTGTACGTACACGGTTTCAGGTTCTTTTTCACTCCCCTCGCCGGGGTTCTTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACTAgg                                                              <  2:2181517/137‑1 (MQ=0)                                                                          ttatACAAAAGGTACGCAGTCACACGCCTAAGCGTGCTCCCACTGCTTGTACGTACACGGTTTCAGGTTCTTTTTCACTCCCCTCGCCGGGGTTCTTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACTAgg                                                              <  2:2485120/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)                                                                                                                              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)                                                                                                                                             |                                                                                                                                     TATTCAGACTCGCTTTCGCTGCGGCTCCGCATCTTCTGCTTAACCTTGCACGGGATCGTAACTCGCCGGTTCATTCTACAAAAGGCACGCCATCACCCG‑TTAA‑CGGGCTCTGACTACTTGTAGGCACACGGTTTCAGGATCTCTTTCACTCCCCT‑TCCGGGGTGCTTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCC  >  NZ_CP034943/2980109‑2980378

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