Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NZ_CP034943	2,978,569	C→T	18.0%	noncoding (2232/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,978,569	0	C	T	18.0%	2.7 / 20.9	56	noncoding (2232/2928 nt)	EO946_RS15430	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (11/34);  new base T (3/7);  total (14/42)
Fisher's exact test for biased strand distribution p-value = 7.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01

GCTCCCTTGTGCCTTTACACTCTGCGAATGATTTCCAACCATTCTGAGGGAACCTTTGGGCGCCTCCGTTACCTTTTAGGAGGCGACCGCCCCAGTCAAACTGCCCACCTGACACTGTCTCCCCGCCCGATAAGGGCGGCGGGTTAGAAGGTCAATACAGCCAGGGTAGTATCCCACCGATGCCTCCACCGAAGCTGGCGCTCCGGTTTCCAAGGCTCCTACCTATCCTGTACAAGCTGTACCAACATTCAATATCAGGCTGCAGTAAAG  >  NZ_CP034943/2978430‑2978699                                                                                                                                            |                                                                                                                                   gcgcccTTGTGCCTGTAGACTCTGCGAATGATTTCCAACCATTCTGAGGGAACCTTTGGGCGCCTCCGTTACCTTTTAGGAGGCGACCGCCCCAGTCAAACTGCCCACCTGACACTGTCTCCCCGCCCGATAAGggcggc                                                                                                                                    <  2:2583731/137‑1 (MQ=0)     ccgtgtgCCTTTACACTCTGCGAATGATTTCCAACCATTCTGAGGGACCCTTTGGGCGCCTCCGGTACCTTTTAGGAGGCGACCGCCCCAGTCAAACTGCCCACCTGACACTGTCTCCCCGCGCGATAAGGGCGGCGGGt                                                                                                                                >  2:234315/4‑140 (MQ=0)        tttggcTTTACACTCTGCGAATGATTTACATCCATTCTGTGGGAACCTTTGGGCGGCTCCGTTACCTTTTAGGAGGCGAACGCCCCAGTCAAACTGCCCACCTGACACTTTCTCCCCGCCCGATAAGGGCGGCGGGTTAg                                                                                                                             <  2:3055086/135‑1 (MQ=1)        tgggCCTTTACACTCTGCGAATGATTTCCAACCATTCTGAGGGAACCTTTGGGCGCCTCCGTTACCTTTTAGGAGGCGACCGCCCCAGTCAAACTGCCCACCTGACACTGTCTCCCCGCCCGATAAGGGCGGTGGGTTAg                                                                                                                             <  2:3281511/137‑1 (MQ=1)        tgggCCTTTACACTCTGCGAATGATTTCCAACCATTCTGAGGGAACCTTTGGGCGCCTCCGTTACCTTTTAGGAGGCGACCGCCCCAGTCAAACTGCCCACCTGACACTGTCTCCCCGCCCGATAAGGGCGGCGGGTTAg                                                                                                                             <  1:2388599/137‑1 (MQ=1)        tgggCCTTTACACTCTGCGAATGATTTCCAACCATTCTGAGGGAACCTTTGGGCGCCTCCGTTACCTTTTAGGAGGCGACCGCCCCAGTCAAACTGCCCACCTGACACTGTCTCCCCGCCCGATAAGGGCGGCGGGTTAg                                                                                                                             <  2:2839028/137‑1 (MQ=1)        ggttccTTTCCAGTCTGCGACTGTTTTCCAAGCATTCTGAGGGAAGCTTTGGGCGCCTCCGTTACGTTTTCGGAGGCGACCGACCCAGTCAAACTGCCCACATGACACTGTCTCCCCGCCGGATTAGGGCGGCGGGTTAg                                                                                                                             <  2:6454511/136‑1 (MQ=1)             ctgtACACTCTGAGAATGATTTCCAACCATTCTGAGGGAACCTATGGGCGCCTCCGTTACCATCTCCGAGGCGACCGCCCCAGTCAAACTGCCCACCTGACACTGTCTCCCCGCCCGATAAGGGCGGCGGGTTAGAAGGt                                                                                                                        <  2:4581747/137‑1 (MQ=1)               ttgcacTCTGCGAATGATTTCCAACCATTCTGAGGGAACCTTTGGGCGCCTCCGTTACCTTTTAGGAGGCGACCGCCCCAGTCAAACTGCCCACCTGACACTGTCTCCCCGCCCGATAAGGGCGGCGGGTTAGAAGGTCa                                                                                                                      >  1:1805113/4‑140 (MQ=1)                 accctctGCGAAGGATTTCCAACAAGTCTGAGGGAACCTTGTGGCGCCTCCGTTACCTTTTAGGAGGCGACCCCCCCAGTCAAACTGCCCCCCTGACACTGTCTCCCCGCCCGATAAGGGCGGTGGGTTAGAAGGTCAAt                                                                                                                    <  2:4178988/137‑1 (MQ=1)                 aaacactGCGAATGATTTACAACCATTATGATTGAACCTTTGGGCGCCTCCGTTAACTTTTAGGAGGCGACCGACCCAGTCAAACTGCCCACCTGACACTGTCTCCCCGCCCGATAAGGGCGGCGGGTTAGAAGGTCAAt                                                                                                                    >  2:4024343/6‑140 (MQ=1)                      cttcGAATGATTTCCAAACATTCTGAGGGAACCTTTGGGGGCCTCCGTTACCTTTTAGGAGGGGACCTCGCCATTCAAACTGCCCACCTGACACTGTCTCCCCGCCCGATAAGGGCGGCGGGTTAGAAGGTCAATACAGc                                                                                                               <  1:4635130/137‑1 (MQ=1)                      ctccGAATGATTTCCAACCAATCTGAGGGAACCTTTGGGCGCCTCCGTTACCTGTGAGGAGGCGACCGCCCCAGTCAAACTGCCAACCTGACACTGTCTCCCCGCCCGAGAACCGCGGTGGGTTAGAAGGTCAATACAGc                                                                                                               >  2:6106752/4‑140 (MQ=1)                       tgagAATGATTTCCAACCATTCTGATTGAACAAATGGGCGCATACGTTACCTTTTAGGAGGCGACCGCCCCAGTCAAACTGCCCACCTGACACTGTCTCCCCGCCCGATAAGGGCGGCGGGTTAGAAGGTCAATACAGcc                                                                                                              >  2:6565692/4‑140 (MQ=1)                       tgagAATGATTTCCAACCATTCTGATTCAACCTAACCGCCCCTCCCTTACCTTTTAGTAGGCCACCGCCCCACTCAAACTGACCACCAGACACTGTCTCCCCGCCCGATAAGGGCGGCGGGTTAGACGGGCAATACAGcc                                                                                                              >  2:6958544/4‑140 (MQ=1)                       tgagAATGATTTCCAACCATTCTGAGGGAACATTAGGGCGCCTCCGTTACCTTTTAGGAGGAGACCGCCCCAGTCAAACTGCCCACCTGACACTGTCTCCCCGCCCGATAAGGGCGGCGGGTTAGAAGGTCAATACAGcc                                                                                                              >  2:2614305/4‑140 (MQ=0)                       tgagAATGATTTCAAAACATTCTGAGGGAACCTTTGGGCGCCTCCGTTACCTTTTAGGAGGCGACCGCCCCAGTCAAACTGCCCACCTGACACTGTCTCCCCGCCCGATACGGGCGGCGGGTTAGAAGGTCAATACAGca                                                                                                              >  2:2957358/4‑139 (MQ=0)                        gctaaTGATTTCCAACCATTCTGAGGGAACCTTTGGGCGCCTCCGTTACCTTTTAGGAGGCGACCGCCCCAGTCAAACTGCCCACCTGACACTGTCTCCCCGCCCGATAAGGGCGGTGGGTTAGAAGGTCAATACAGCCa                                                                                                             <  1:2495862/137‑1 (MQ=1)                        gccaaTGATTTCCAACCATTCTGAGGGAACCTTTGGGCGCCTCCGTTACCTTTTAGGAGGCGACCGCCCCAGTCAAACTGCCCACCTGACACTGTCTCCCCGCCCGATAAGGGCGGCGGGTTACAAGGTCAAGACAGCCa                                                                                                             <  1:5608274/137‑1 (MQ=0)                        gccaaTGATTTCCAACCAGTCTGAGGGAACCTTTGGGCGCCTCCGTTACCTTTTAGGAGGCGACCGCCCCAGTCAAACTGCCCACCTGACACTGTCTCCCCGCCCGATAAGGGCGGCGGGTTAGAAGGTCAATACAGCCa                                                                                                             <  2:6024282/137‑1 (MQ=1)                        gagattGATGTCCAACCCTTCTGAGGGAACCTTTGGGCGCCTCCGTTACCTTTTAGGAGGCGACCGCCCCCCTCCACCTGCCCACCTGAGACTGTCTCCCCGCCCGATAAGGGAGGCGGGTTAGAAGTTCAATACAGCCa                                                                                                             <  2:3107437/135‑1 (MQ=1)                           aacgATTTCCAACCATTCTGAGGGAACCTTTGGGCGCCTCCGTTACCTTTTAGGAGGCGACCGCCCCAGTCAAACTGCCCACCTGACACTGTCTCCCCGCCCGATAAGGGCGGCGGGTTAGAAGGTCAATACAGCCAggg                                                                                                          <  1:6635305/137‑1 (MQ=1)                             tgatggccAACCATTCTGAGGGAACCTTTGGGCGCCTCCGTTACCTTTGAGGAGGCGACCGCCCCAGTCAAACTGCCCACCTGACACTGTCTCCCCGCCCGATAAGGGCGGTGGGTTAGAAGGTCAATACAGCCAGGgta                                                                                                        <  2:7040559/134‑1 (MQ=0)                              gagttCCACCCATTCTGAGGGAACCTTTGGGCGCCTCCGTTACCTTTTAGGAGGCGACCGCCCCAGTCAAACTGCCCACCTGACACTGTCTCCCCGCCCGATAAGGGCGGCGGGTTAGAAGGTCAATACAGCCAGGgtag                                                                                                       <  1:5700536/137‑1 (MQ=1)                               atatCCAACCATTCTGAGGGAACATTTGTGACCCTCCGTTACCTTTTAGGAGGCGACCGCCCCAGTCAAACTGCCCACCTGACACTGTCTCCCCGCCCGATAAGGGCGGTGGGTTAGAAGGTCAATACAGCCAGGgtagt                                                                                                      >  2:2844827/4‑140 (MQ=1)                                   cccaCCATTCTGAGGGAACCTTTGGGCGCCTCCGTTACCTTTTAGGAGGCGACCGCCCCAGTCGAACTGCCCACCTGACACTGTCTCCCCGCCCGATAAGGGCGGCGGGTTAGAAGGTCAATACAGCCAGGGTAGTAtgc                                                                                                  <  2:5483332/137‑3 (MQ=0)                                   cccaCCAGTCTGAGGGAACCATTGGGCGCCTCCGTTACCTTTTAGGAGGCGACCGCCCCAGTCAAACTGCCCACCTGACACTGTCTCCCCGCCCGATAAGGGCGGCGGGTTAGAAGGTCAATACAGCCAGGGTAGTATcc                                                                                                  <  1:5009561/137‑1 (MQ=0)                                       cccttCTGAGGGAACCTTTGGGCGCCTCCGTTACCTTTTAGGAGGCGACCGCCCCAGTCAAACTGCCCACCTGACACTGTCTCCCCGCCCGATATGGGCGGCGGGTTAGAAGGTCAATACAGCCAGGGTAGTATCCCAcc                                                                                              <  2:1835807/137‑1 (MQ=1)                                          gtcggAGGGAACCGTGGGGCGCCGCCGTTCCCTTGTCGGAGGCGACCGCCCCAGTCACACTGCCCACCTGACACTGTCTCCCCGCCCGATAAGGGCGGCGGGTTAGAAGGTCAATACAGCCTGGGTAGTATCCCACCGAt                                                                                           <  2:4613229/136‑1 (MQ=1)                                           tcggAGGGAACCTTTGGGCGCCTCCGTTACCTTTTAGGAGGCGACCGCCCCAGTCAAACTGCCCACCTGACACTGTCTCCCCGCCCGATAAGGGCGGCGGGTTAGAAGGTCAATACAGCCAGGGTCGTATCCCACCGATg                                                                                          >  2:1579620/4‑140 (MQ=1)                                                    acatttGGGCGCCTCCGTTACCTGTTAGGAGGCGACCGCCCCAGTCAAACTGCCCACCTGACCCTGTCTCCCCGCCCGATAAGGGCGGCGGGTTAGAAGGTCAATACAGCCAGGGTAGTCAGCCACCGATGCCTCCACCg                                                                                 >  2:5340072/4‑140 (MQ=1)                                                                cacccttACATATTAGGAGTCGACCGCACCAGTCAAACTGCCCACCTGACACTCTCTCCCCGCCCGAAAAGGGCGGCGGGATAGAAGGTCAATACAGCCAGGGTAGTATCCCACAGATGCCTCCACCGAAGATGGCGCTc                                                                     >  2:186932/6‑140 (MQ=1)                                                                         ctgttAGGAGGCGACCGCCCCAGTCACACTGCCCACCTGACACTGTCTCCCCGCCCGATAAGGGCGGCGGGTGAGAAGGTCAATACAGCCCGGGTAGTATCCCACCGATGCCTCCACCGAAGCTGGCGCTCCGGTTTCCa                                                            <  1:4479147/137‑1 (MQ=0)                                                                         ctgttAGGAGGCGACCGCCCCACTCAAACTGCCCACCTGACACTGTCTCCCCGCCCCACACGGGCGGCGGGTTAGAAGGTCAATCCAGCCAGGGTAGTATCCCACCGATGCCTCCACCGAAGCTGGCGCTCCGGTTTCCa                                                            <  2:4942386/137‑1 (MQ=1)                                                                          ttgtAGGAGGCGACCGCCCCCGTCAAACGGCCCACCTGACACTGTCTCCCCGCCCGATAAGGGCGGCGGGTTAGAAGGTCAATACAGCCAGGGTAGTATCCCACCGATGCCTCCACCGAAGCTGGCGCTCCGGTTTCCaa                                                           <  2:6218085/137‑1 (MQ=0)                                                                          ttgtAGGAGGCGACCGCCCCAGTCAACCTGCCCACCTGACACTGTCTCCCCGCCCGATAAGGGCGGTGGGTTAGAAGGTCAATACAGCCAGGGTAGTATCCCACCGATGCCTCCACCGAAGCTGGCGCTCCGGTTTCCaa                                                           <  2:6601859/137‑1 (MQ=0)                                                                            tttggaggCGACCGCCCCAGTCAAACTGCCCACCTGACACTGTCTCCCCGCCCGATAAGGGCGGCGGGTTAGAAGGTCAATACAGCCAGGGTAGTATCCCACCGATGCCTCCACCGAAGCTGGCGCTCCGGTTTCCAAgg                                                         <  1:590724/137‑1 (MQ=1)                                                                              agcaggCGACCGCCCCACTCAAACTGCCCACCTGACACTGTCTCCCCGCCCGATACGGGCGGCGGGTTAGAAGGTAAATACAGCCAGGGTAGTATCCCACCGATGCCTCCACCGAAGCTGGCGCTCCGGTTTCCAAGGCt                                                       <  1:7613474/137‑1 (MQ=0)                                                                              acgacgCGACCGCCCCAGTCACACTGCCCACCTGACACGGTCTCCCCGCCCGATCACCCCGGCGGGTTAGAAGGTCAATACAGCCAGCGTAGTATCCCACCGCTGCCTCCACCGAAGCTGGCGCTCCGGTTTCCAAGGCt                                                       <  2:4243335/135‑1 (MQ=1)                                                                                        cggcccAGTCAAACTTCGCACCTGACACTGTCTCCCCGCCCGATAAGGGCGGCGGGTTAGATGGTCAATACACCCAGGGTAGTATCCCACCGATGCCTCCACCGAAGCTGGCGCTCCGGTTTCCAAGGCTCCTACCTATc                                             <  2:4238484/137‑1 (MQ=1)                                                                                        cggcccAGTCAAACTGCCCACCTGACAGTGTCTCCCCGCCCGATAAGGGCGGCGGGTTAGAAGGTCAATACAGGCAGGGTAGTATCCCACCGATGCCTCCACCGAAGCTGGCGCTCCGGTTTCCAAGGCTCCTACCTATc                                             <  2:3022430/137‑1 (MQ=0)                                                                                        cggcccAGTCAAACTGCCCACCTGACACTGTCTCCCCGCCCGATTAGGGCGGCGGGTTAGAAGGTCAATACAGCCCGGGTAGTATCCCACCGATGCCTCCACCGAAGCTGGCGCTCCGGTTTCCAAGGCTCCTACCTATc                                             <  2:1196729/137‑1 (MQ=0)                                                                                        cgacccAGTCAAACTGCCCACCTGACACTGTCTCCCCGCCCGATAAGGGCGGCGGGTTAGAAGGTCAATACAGCCAGGGTAGTATCCCACCGATGCCTCCACCGAAGCTGGCGCTCCGGTTTCCAAGGCTCCTACCTATc                                             <  2:1393657/137‑1 (MQ=1)                                                                                             cactCACACTGCCCACCTGACACTGTCTCCCCGCCCGCGCAGGGCGGAGGGTTAGAAGGTCAAGACAGCCAGGGTAGTATCCCACCGATGCCTCCACCGAAGCTGGCGCTCCGGTTTCCAAGGCTCCTACCTATCCTGTa                                        <  2:1251128/137‑1 (MQ=1)                                                                                             cactCAAACTGCCCACCTGACACTGTCTCCCCGCCCGATAAGGGCGGCGGGTTAGAAGGTCAATACAGCCAGGGTAGTATCCCACCGATGCCTCCACCGAAGCTGGCGCTCCGGTTTCCAAGGCTCCTACCTATCCTGTa                                        <  1:7639509/137‑1 (MQ=1)                                                                                               ttccaaCTGCCCACGGGACACGGTCTCCCCGCCCGATAAGGGCGGCTGGTTCGAAGGTAAATAGAGCCAGGGTAGTATCGCACCGATGCCTCCACCGAAGCTGGCGCTCCGGTTTCCAAGTGTCCTACCTATCCTGTACa                                      <  2:4426218/136‑1 (MQ=1)                                                                                                 cacacggCCCACCTGACACTGTCTCCCCGCCCGATAAGGGCGGCGGGTTAGAAGGTCAATACAGCCAGGGTAGTATCCCACCGATGCCTCCACCGAAGCTGGCGCTCCGGTTTCCAAGGCTCCTACCTATCCTGTACAAg                                    <  2:1757096/134‑1 (MQ=1)                                                                                                 cacaCTGCCCACCTGACACTGTCTCCCCGCCCGATAAGGGCGGTGGGTTAGAAGGTCAATACAGCCAGGGTAGTATCCCACCGATGCCTCCACCGAAGCTGGCGCTCCGGTTTCCAAGGCTCCTACCTATCCTGTACAAg                                    <  2:2219691/137‑1 (MQ=1)                                                                                                  aaccTGCCCACCTGACACTGTCTCCCCGCCCGATAAGGGCGGCGGGTTAGAAGGTCAATACAGCCAGGGTAGTATCCCACCGATGCCTCCACCGAAGCTGGCGCTCCGGTTTCCAAGGCTCCTACCTATCCTGTACAAGc                                   <  1:2603101/137‑1 (MQ=1)                                                                                                    acggCCCACCTGACACTGTCTCCCCGCCCGATAAGGGCGGCGGGTTAGAAGGTCAATACAGCCAGGGTAGTATCCCACCGATGCCTCCACCGAAGCTGGCGCTCCGGTTTCCAAGGCTCCTACCTATCCTGTACAAGCTg                                 <  2:1378055/137‑1 (MQ=1)                                                                                                         cccccTGACACTGTCTCCCCGCCCGATAAGGGCGGTGGGATAGAACGTCAAGACAGCCAGGGTAGTATCCCACCGATGCCTCCACCGAAGCTGGCGCTCCGGTTTCCAAGGCTCCTACCTATCCTGTACAAGCTGTACCa                            <  2:3932989/137‑1 (MQ=1)                                                                                                         cccccTGACACTGTCTCCCCGCCCGATAAGGGCGGCGGGTTAGAAGGTCAATACCGCCAGGGTAGTATCCCACCGATGCCTCCACCGAAGCTGGCGCTCCGGTTTCCAAGGCTCCTACCTATCCTGTACAAGCTGTACCa                            <  2:2417398/137‑1 (MQ=1)                                                                                                                    tttcaccccGCCCGATAAGGGCGTTTTGTTAGAAGGTCAATACAGCCAGGGTAGTATCCCACCGATGCCTCCACCGAAGCTGGCGCTCCGGTTTCCAAGGCTCCTACCTATCCTGTACAAGCTGTACCAACATTCAatat                 >  2:3373572/6‑140 (MQ=1)                                                                                                                                cgttAAGGGCGGCGGGTTCGAAGGTCAATACAGCCAGGGTAGTATCCCACCGATTCCTCCACCGAAGCTGGCGCTCCGGTTTCCAAGGCTCCTCCCTATCCTGTACAAGCTGTACCAACCTTCAATATCAGGCTGCAGTa     >  2:7626844/4‑140 (MQ=1)                                                                                                                                cgatccgggCGGCGGGTTAGAAGGTCAATACAGCCAGGGTAGTATCCCACCGATGCCTCCACCGAAGCTGGCGCTCCGGTTTCCAAGGCTCCTACCTATCCTGTACAAGCTGTACCAACATTCAATATCAGGCTGCAGTa     <  2:7317310/134‑1 (MQ=1)                                                                                                                                   tacgggCGGCGGGTTAGAAGGTCAATACAGCCAGGGTAGTATCCCACCGATGCCTCCACCGAAGCTGGCGCTCCGGTTTCCAAGGCTCCTACCTATCCTGTACAAGCTGTACCAACATTCAATATCAGGCTGCAGTAAAg  <  2:3035256/137‑1 (MQ=1)                                                                                                                                   tacgggCGGCGGGTTAGAAGGTCAATACAGCCAGGGTAGTATCCCACCGATGCCTCCACCGAAGCTGGCGCTCCGGTTTCCAAGGCTCCTACCTATCCTGTACAAGCTGTACCAACATGCAATATCAGGCTGCAGTAAAg  <  2:7500158/137‑1 (MQ=1)                                                                                                                                            |                                                                                                                                   GCTCCCTTGTGCCTTTACACTCTGCGAATGATTTCCAACCATTCTGAGGGAACCTTTGGGCGCCTCCGTTACCTTTTAGGAGGCGACCGCCCCAGTCAAACTGCCCACCTGACACTGTCTCCCCGCCCGATAAGGGCGGCGGGTTAGAAGGTCAATACAGCCAGGGTAGTATCCCACCGATGCCTCCACCGAAGCTGGCGCTCCGGTTTCCAAGGCTCCTACCTATCCTGTACAAGCTGTACCAACATTCAATATCAGGCTGCAGTAAAG  >  NZ_CP034943/2978430‑2978699

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