Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NZ_CP034943	2,982,148	G→T	18.9%	noncoding (372/1550 nt)	EO946_RS15435 ←	16S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP034943	2,982,148	0	G	T	18.9%	5.7 / 19.5	68	noncoding (372/1550 nt)	EO946_RS15435	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (4/50);  new base T (3/10);  total (8/60)
Fisher's exact test for biased strand distribution p-value = 1.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGCTTTCTGGTTAGGTACCGTCAAGGTACCGCCCTATTCGAACGGTACTTGTTCTTCCCTAACAACAGAGCTTTACGATCCGAAAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCCGATCACCCTCTCAGGTCGGCTACGCATCGTTGCCTTGGTGAGCCATTACCTCACCAACTAGCTAATGCGCCGCGGGTC  >  NZ_CP034943/2982010‑2982283                                                                                                                                           |                                                                                                                                        tgtcTTTCTGGTTAGGTACCGTCAAGGTACCGCCCTATTCGAACGGTACTTGTTCTTCCCTACCAACAGAGCTTTACGATCCGAAAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGa                                                                                                                                        <  2:2334433/137‑1 (MQ=1) tgtcTTTCTGGTTAGGTACCGTCAAGGTACCCCCCTATTCGAACGGTACTTGTTCTTACCTAACAACAGAGCTTTACGATCCGAAAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGa                                                                                                                                        >  1:4302192/4‑140 (MQ=0) gggcgtgcTGGATAGGTACCGTCAAGGTACCGCCCTAGACCAACGGAACTTGTACTTCCCTAACAACAGAGCTTTACGATCGGACAACCGTCATCACTCCCGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGa                                                                                                                                        <  1:5016158/133‑1 (MQ=1)   gcgttCTGGTTAGGTACCGTAACGGTACCGCCCTTTTCGCACTGTTCTTGTTGTTCCCTAACAAGAGAGCTTTACGATCCGAAAACCTTCATCACTCAGGCGGCGGTGCTCCTTCAGACTTTCGTCCATTGCGGAAGAtt                                                                                                                                      <  2:5265075/137‑1 (MQ=1)        cttgTTAGGTACCGTCAAGGTACCGCCCTATTCGAACGGTACTTGTTCTTCCCTAACAACAGAGCTTTACGATCCGAAAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTa                                                                                                                                 >  1:6839940/4‑140 (MQ=1)             gagctACCGTCAAGGTACCGCCCTATTCGAACGGTACTTGTTCTTCCCTAACAACAGAGCTTTACGATCCGAAAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTActgct                                                                                                                            <  2:7719650/136‑1 (MQ=1)              cgggaCCGTCAAGGTACCGCCCTATTCGAACGGTACTTGTTCTTCCCTACCAACAGAGCTTTACGATCCGAAAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTActgctg                                                                                                                           <  2:4313006/136‑1 (MQ=0)              agttACCGTCAAGGTTCCGCACTATTAGAACGGTACTTGTTCTTCGCTTTCAACAGATCTTTACGATCCGAAAACCTTCATCTCTCACGCGGCGTTGCTCCGTCTGTTTTTCTTTATTTTCGTATGATTCCCTActtgtg                                                                                                                           <  2:7738121/137‑5 (MQ=0)              agttACCGTCAAGGTACCGCACAATTCGAACGGTACTTGTTCTTCCCTAACAACAGAGCTTTACGATCCGAAAAACTTAATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTActgctg                                                                                                                           <  2:2714304/137‑1 (MQ=1)              aggttacGTCAAGGTACCGCCCTATTCGAACGGTACTTGTTCTTCCCTAACAACAGAGCTTTACGATCCGAAAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTActgctg                                                                                                                           <  2:6433105/134‑1 (MQ=1)                gttccGTCAAGGTACCGCCATATTCGAAAGGTAATTTTTCTTACCTAACAACATAGCTTTATGATCCGAAAAGCTTCTTCACTAACGCGGCGTTGCTCCGTTAGACTTTCGTCCATTGCGGAATATTCCCTACTGCTGcc                                                                                                                         <  2:7468129/137‑1 (MQ=1)                     gttaaGGTCCCGCCCTATTCGAACGGTACTTGTTCTTCCCTAACAACAGAGCTTTACGATCCCAAAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAATATTCCCTACTGCTGCCTCCCg                                                                                                                    >  2:5408291/4‑140 (MQ=0)                       catggTACCGCCCTATTCGCACGGTACGTGTTCTTCCCTAACACCAGAGCTGTCCGATCCGAAAACCTTCATCACTCACCCGGCGTTGCTCCGTCAGACGTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTa                                                                                                                  <  2:3364310/137‑1 (MQ=1)                       catggTACCGCCCTATTCGACCGGTACTTGTTCTTCCCTAACACCAGAGCTTTACGATCCGAAAACCTTCATCACTCACGCTGCGTTGCTCCGTCAGACTTTCGTCCATTGCGTAAGATTCCCTACTGCTGCCTCCCGTa                                                                                                                  <  2:5710315/137‑1 (MQ=1)                       cacggTACCGCCCTATTCGAACGGTACTTGTTCTTCCCTAACAACAGAGCTTTACGATCCGAAAACCTTGATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTa                                                                                                                  <  1:7043588/137‑1 (MQ=1)                       caagtgaCCGCCCTATTCGAACGGTACGTGTTCTTCCCTAACAACAGAGCTTTAGGATCCGAAAACCTTGATCACTGACGCGGCGTTGCTCGGTCAGACTTGCCTCCATTGCGGAAGTTTACCTACTGCTGCCTCCCCTa                                                                                                                  <  1:4831669/134‑1 (MQ=1)                       aaacgTACCGCCCTATTCGAACGTTACTTGTTCTTACCTAACAACAGAGCTTTACGATCCGAAAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTa                                                                                                                  <  1:4417676/136‑1 (MQ=0)                        aagggtccGACCTATTAGAACGGTAATTGTTCTTCCCTAACAACAGAGCTTTACGATCCGAAAACCTTCATGACTCTCGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTAg                                                                                                                 <  1:1612112/134‑1 (MQ=1)                           gtacatcccTATTCGAAAGGTAATTTTTATAAACTAACAACAGAGCTAAACCATCCTAAAACCTTCATAACTCACGCGGCCTAGCTCCGTCAGACTTTCTTCCATTTAGGACAATTCCCTACTGCTGCCTCCCGTATGAg                                                                                                              >  2:6382907/7‑140 (MQ=0)                                 ccatATTCGAACGGTACTTGTTCTTCCCTAACAACAGAGCTTTACGATCCGAAAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTggg                                                                                                        <  1:5945063/137‑1 (MQ=1)                                    tagtCGAACGGTACTTGTTCTTCCCTAACAACAGAGCTTTACGATCCGAAAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCg                                                                                                     <  1:4623431/137‑1 (MQ=1)                                       tcgcagggTACTTGTTCTTCCCTAAGAAGAGAGCTTTACGATCCGAAAACCTTCATCACTCACGCGGCGTTTCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCgtgt                                                                                                  <  1:295176/134‑1 (MQ=1)                                         gaccGGTACTTGTTCTTCCCTCACAACCGAGCTGTACGATCCGAAAACCGTCATCACTCACGCGGCGGTGCTCCGTCAGCCTTTCGTCCCTTGCGGAAGATTCCCGCCTGCTGCCTCCCGGAGGAGTCTGGGCCGTGtct                                                                                                <  2:5750610/137‑1 (MQ=1)                                         gaccGGTACTTGTTCTTCCCTAACACCAGAGCTTTACGATCCGAAAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGtct                                                                                                <  1:7652076/137‑1 (MQ=1)                                         gaccGGTACTTGTTCTTCCCTAACAACAGAGCTTTACGATCCGAAAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGtct                                                                                                <  2:2850456/137‑1 (MQ=1)                                         gaccGGTACTTGTTCTTCCCTAACAACAGAGCTTTACGATCCGAAAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGtct                                                                                                <  2:4699424/137‑1 (MQ=1)                                         gaccGGTACTTGTTCTTCCCTAACAACAGAGCTTTACGATCCGAAAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGtct                                                                                                <  1:4669546/137‑1 (MQ=1)                                         gaccGGTACTTGTTCTTCCCTAACAACAGAGCTTTACGATCCGAAAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGtct                                                                                                <  2:1790950/137‑1 (MQ=1)                                         gaccGGTACTTGTTCTTCCCTAACAACAGAGCTTTACGATCCGAAAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCAGTGCGGAATAGTCCCTACTGCGGCCTCCCGTAGGAGTCTGGGCCGTGtct                                                                                                <  2:2997442/137‑1 (MQ=1)                                         gaccGGTACTTGTTCTTCCCTAACAACAGAGCTTTACGATCCGAAAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTGTCGTCCAGTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGtct                                                                                                <  2:4179321/137‑1 (MQ=0)                                         gaccGGTACTTGTTCTTCCCGACCAACAGAGCTTTACGATCCGAAAACCTTGATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAATATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGtct                                                                                                <  1:6587921/137‑1 (MQ=1)                                         gaccGGTACTTGTTCTTCCCGACCAACAGAGCTTTACGATCCGAAAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGtca                                                                                                <  2:3824388/137‑2 (MQ=0)                                         gaccGGTACTTGGTCTTCCCTACCAACAGAGCTTGACGATCCGAAAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAATATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGtct                                                                                                <  2:6442589/137‑1 (MQ=1)                                         gaccGGTACTTGGTCTTCCCTAACAACAGAGCTTTACGATCCGAAAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGtct                                                                                                <  2:6400080/137‑1 (MQ=1)                                         gaccGGTACGTGGTCGTCCCTAACAACAGAGCTTGACGATCCGAAACCCTTCATCACTCACGCGGCGGTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGtct                                                                                                <  2:7567585/137‑1 (MQ=1)                                         gaccGGGACTTGGTCGTACCGAACTACAGAGCTTTACGATCAGAAAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCAGTGCGGAAGATTCCCGCCTGCTGCCTCCCGTAGGAGTCTGGGCCGTGtct                                                                                                <  1:5597383/137‑1 (MQ=1)                                             gggaCTTGTTCTTCCCTAACAACAGAGCTTTACGATCCGAAAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAATATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGt                                                                                            <  1:1233795/137‑1 (MQ=1)                                             gggaCTTGGTCTTCCCTAACAACAGAGCTTTACGATCCGAAAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGt                                                                                            <  2:7739636/137‑1 (MQ=1)                                               tagttGTTCTTCCCTAACAACAGAGCATTACGATCCGAACACCTTCATCACTCCCGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTCCTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTcc                                                                                          <  2:2528029/137‑1 (MQ=1)                                               gacgtGTTCTGCCCTAACACCCGAGCTTTACGATCCGAAAACCGTCATCCCTCCCGCGGCGGTGGACCGTCAGACTTCCCTCCATCGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTcc                                                                                          <  2:4954970/136‑1 (MQ=1)                                                  tttttctgcCCTAACAACACAGCTTTACGATCCGAAAACCGTCATCACTCACGCGGCGGTGCTGGCTCAGAGTTTCGTGGATGGCTGAATATTCCCTACTGCTGCCTCCCGTTTGAGTCTGGGCCGTGTCTGAGTCCCAg                                                                                       <  2:4485901/137‑1 (MQ=1)                                                     ttattcCCTAACAACAGAGCTTTACGATCCGAAAACCATCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGAGTAAGATTACCTACTGCTGCCTCCCCTACGAGTCTCGGCCGTGTCTCAGTCCCAgtgt                                                                                    >  2:876704/4‑140 (MQ=1)                                                      tcgtcCCTACCTACAGAGCTTTACGATCCGAAAACCTTCATAACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGGTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGtgtg                                                                                   <  2:3736999/137‑1 (MQ=1)                                                      tcgtcCCTAACAACAGAGCTTTACGATCCGAAAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGtgtg                                                                                   <  2:5466386/137‑1 (MQ=1)                                                         tcccgcacaacaGAGCTTTACGATCCGAAAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGGAGGTGTCTGGGCCGTGTCTCAGTCCCAGTTTGGcc                                                                                >  2:2281386/7‑140 (MQ=1)                                                           ccgaacaccAGAGCGTTACGATCCGAAAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGACCTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCCGTCCCAGTGTGGCCGa                                                                              <  1:2661424/137‑1 (MQ=1)                                                           ccgaacaacATAGCTTTACTATCCGAAAACCATCATCACTAACGCTGCGTTGATCCGTCAGAATTTCTTCAATTGCGGAATATTCCCTACTGCTGCCTCCCGTAGGAGTCAGGGCCGGGTCTCCGTCCCAGTGAGGACGa                                                                              >  2:7271869/4‑140 (MQ=1)                                                            ctcacaccaGAGCTCGACGCTCCGACACCCGTCATCACTCACGCGGCGGTGCTCCGTCAGTCTTTCGTCCATTGCGGAAGATTCCCTCCTGCGGCCTCCCGTAGGAGTCTGGGCCGTGTCTCACTCCCACTGTGGCCGAt                                                                             <  2:4297545/137‑1 (MQ=1)                                                             taccaacaGAGCTTTCCGATCCGAAAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGCCTGTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCCGATc                                                                            <  1:2839728/137‑1 (MQ=0)                                                             taccaacaGAGCTTTACGATCCGAAACCCTTCATCACTCACGCGGCGTTGCTGCGTCAGACTTTCGTCCATTGAGGAAGATTCCCTACTGCTGCCTCGCGTTGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCCGATc                                                                            <  2:1352682/137‑1 (MQ=1)                                                             taccaacaGAGCTTTACGATCCGAAAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAATATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCCGATc                                                                            <  1:5900882/137‑1 (MQ=1)                                                             taccaacaGAGCTTTACGATCCGAAAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCCGATc                                                                            <  2:4563897/137‑1 (MQ=1)                                                             taccaacaGAGCTTTACGATCCGAAAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCCGATc                                                                            <  1:6849108/137‑1 (MQ=1)                                                              cccaccaGAGCTCTACGATCCGAAAACCTTCATCACCCCCGCCGCGTTGCTCCGTCAGACGTGCCTCCATTGCGGAAGAGTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCCGCTCa                                                                           <  2:106641/135‑1 (MQ=1)                                                               accacaGAGCGTGACGATCCGAATACCTGCATCACTCCCGCGGCGTTGCTCCGTCAGACATTCGTCCATTCCGGAAGATTCCCTACTGCTGCCTCCCTTAGGAGTGTGGGCCGTGTCTCAGTCCCAGTGTGGCCGATCAc                                                                          <  2:4263044/137‑1 (MQ=1)                                                                      aggtttCCGATCCGAAAACCTTCATCACTCACGCGGCGTTGCGCCGTCAGACGTTCGTCCATTGCGGAATATTCCCTACTGCTGCCTCCCGTGGGAGTCTGGGCCGTGTCTCAGTCCCGGTGTGGCCGCTCACCGtctct                                                                   >  2:6690800/4‑139 (MQ=1)                                                                       gccttACGATCCGAAACCCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCCGATCACCCTCTCAg                                                                  <  2:4252318/137‑1 (MQ=1)                                                                       gcattACGATCCGAAAACCTTCCTCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCCGATCACCCTCTCAg                                                                  <  2:2858916/137‑1 (MQ=1)                                                                         ttgaCGATCCGACACCCTTCATCACTCACGCGGCGTTGCGCCGTCAGACTGGCGTCCAGTGCGGAAGATTCCCTACGGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCCGATCACCCTCTCCGGt                                                                <  2:4073141/137‑1 (MQ=1)                                                                                      aagcTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAATATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCCGATCACCCTCTCAGGTCGGCTACGCATCg                                                   <  1:5620863/137‑1 (MQ=1)                                                                                                acgcacGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGGGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCCGATCACCCTCTCAGGTCGGCTACGAATCGTCGCCTTGGt                                         <  2:6037727/137‑1 (MQ=0)                                                                                                acgcacGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCCGATCACCCTCTCAGGTCGGCTACGCATCGTCGCCTTGGt                                         <  2:5927573/137‑1 (MQ=1)                                                                                                               ctgcGTCAGACTTTCGTCCATTGCGGAATATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCCGATCACCCTCTCAGGTCGGCTACGCATCGTTGCCTTGGTGAGCCATTACCTCAc                          <  2:5345833/137‑1 (MQ=1)                                                                                                                  cgacAGACTCTCGTCCAGTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGACTCGGGGCCGTGTCTCAGTCCCAGTGTCGCCGATCACCCTCTCAGGTCGGCTACGCAACGGTGCCTTGGTGAGCCATTACCTCCCCaa                       <  2:5720961/137‑1 (MQ=1)                                                                                                                      aggcTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTAGGTGAGCCGTTACCCCACCTACTAg                   <  1:1766021/137‑1 (MQ=0)                                                                                                                            tcttCCATTGAGTAAGATTCCCTTCTTCTGACTCCCGTAGGATTCTGGGGATTGTCTCAGTCGGTGTTTGGCCGATCACTCTCTCATGTCGTCTTCGCATATTTGGATTGGTGAGCCATTACCTCTCCAACTAGCTAATg             <  2:2512735/137‑1 (MQ=0)                                                                                                                                 cagtGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGGCTGGGCCGTGTCTCAGTCCCAGTGTGGCCGATCACCCTCTCAGGTCGGCTACGCATCGTTGCCTTGGTGAGCCATTACCTCACCAACTAGCTAATGcgccg        <  2:2228395/137‑1 (MQ=1)                                                                                                                                       ggcaGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCCGATCACCCTCTCAGGTCGGCTACGCATCGTTGCCTTGGTGAGCCATTACCTCACCAACTAGCTAATGCGCCGCGGGTc  <  2:7365959/137‑1 (MQ=1)                                                                                                                                           |                                                                                                                                        TGGCTTTCTGGTTAGGTACCGTCAAGGTACCGCCCTATTCGAACGGTACTTGTTCTTCCCTAACAACAGAGCTTTACGATCCGAAAACCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCCGATCACCCTCTCAGGTCGGCTACGCATCGTTGCCTTGGTGAGCCATTACCTCACCAACTAGCTAATGCGCCGCGGGTC  >  NZ_CP034943/2982010‑2982283

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