Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP009974	4,147,204	A→G	intergenic (+456/‑565)	RPPX_18425 → / → RPPX_18430	enoyl‑CoA hydratase/conjugal transfer protein TraR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009974	4,147,204	0	A	G	98.3%	89.7 / ‑4.4	59	intergenic (+456/‑565)	RPPX_18425/RPPX_18430	enoyl‑CoA hydratase/conjugal transfer protein TraR
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (30/28);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

GTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTC‑‑‑‑‑‑‑GCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAAGTGGCGCATAGTCACCTTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCC  >  CP009974/4147063‑4147348                                                                                                                                                                          |                                                                                                                                                 gTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTcgcg                                                                                                                                                >  2:5289496/1‑151 (MQ=11) gTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTcgcg                                                                                                                                                <  2:5526361/151‑1 (MQ=11) gTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTcgcg                                                                                                                                                >  1:5520189/1‑151 (MQ=11)  tAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCgcgc                                                                                                                                               <  2:4974862/151‑1 (MQ=11)   aaGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCa                                                                                                                                              >  2:4810753/1‑151 (MQ=11)   aaGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCa                                                                                                                                              <  2:1366103/151‑1 (MQ=11)   aaGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCa                                                                                                                                              <  1:6402187/151‑1 (MQ=11)   aaGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTC‑‑‑‑‑‑‑GCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTTCTTGCTCCTTGGATGTGGGCAGTCACATCGAGTAGGTATGCATGCGCGCa                                                                                                                                              >  2:1563504/1‑151 (MQ=11)     gggCGAAAAGGTGAGTCAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGCTGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCagg                                                                                                                                            >  2:1426314/1‑149 (MQ=11)     gggCGAAAAGGTGAGTCAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCagg                                                                                                                                            <  2:2128141/151‑3 (MQ=11)      ggCGAAAAGGTGAGTCAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCaggt                                                                                                                                           >  1:2159210/1‑148 (MQ=11)      ggCGAAAAGGTGAGTCAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGt                                                                                                                                           <  1:3509786/151‑1 (MQ=11)      ggCGAAAAGGTGAGTCAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGt                                                                                                                                           >  1:4452375/1‑151 (MQ=11)       gCGAAAAGGTGAGTCAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTg                                                                                                                                          >  1:5355779/1‑151 (MQ=11)        cGAAAAGGTGAGTCAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTgg                                                                                                                                         <  2:989232/151‑1 (MQ=11)        cGAAAAGGTGAGTCAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTgg                                                                                                                                         >  1:2088857/1‑151 (MQ=11)         gAAAAGGTGAGTCAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGgc                                                                                                                                        >  2:3016620/1‑151 (MQ=11)         gAAAAGGTGAGTCAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGgc                                                                                                                                        <  2:5688595/151‑1 (MQ=11)           aaaGGTGAGTCAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGgcgc                                                                                                                                      >  2:6432614/1‑151 (MQ=11)            aaGGTGAGTCAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGgcgct                                                                                                                                     >  1:5793897/1‑150 (MQ=11)            aaGGTGAGTCAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCa                                                                                                                                     <  2:6498490/151‑1 (MQ=11)             aGGTGAGTCAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCAt                                                                                                                                    <  1:6171686/151‑1 (MQ=11)                 gAGTCAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTc                                                                                                                                >  2:2679186/1‑151 (MQ=11)                   gTCAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCAc                                                                                                                              >  1:3852058/1‑151 (MQ=11)                   gTCAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCAc                                                                                                                              >  2:2054799/1‑151 (MQ=11)                   gTCAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCAc                                                                                                                              >  1:3894096/1‑151 (MQ=11)                    tCAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCAcc                                                                                                                             <  1:1906979/151‑1 (MQ=11)                     cAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCt                                                                                                                            >  2:6497106/1‑151 (MQ=11)                     cAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCt                                                                                                                            <  2:1978736/151‑1 (MQ=11)                     cAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCt                                                                                                                            >  1:4402424/1‑151 (MQ=11)                     cAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCt                                                                                                                            >  1:4595830/1‑151 (MQ=11)                       gCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCttt                                                                                                                          >  1:127402/1‑151 (MQ=11)                        cGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTc                                                                                                                         <  2:2024451/151‑1 (MQ=11)                         gTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTcc                                                                                                                        <  1:1235597/151‑1 (MQ=11)                         gTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTcc                                                                                                                        <  2:6064884/151‑1 (MQ=11)                         gTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTcc                                                                                                                        <  2:3717757/151‑1 (MQ=11)                         gTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTcc                                                                                                                        <  2:3629621/151‑1 (MQ=11)                         gTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTcc                                                                                                                        <  2:2885875/151‑1 (MQ=11)                         gTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTcc                                                                                                                        <  2:6485117/151‑1 (MQ=11)                          tCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTCCg                                                                                                                       <  1:1806109/151‑1 (MQ=11)                           cGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTCCGc                                                                                                                      >  1:552540/1‑151 (MQ=11)                           cGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTCCGc                                                                                                                      >  1:2039499/1‑151 (MQ=11)                           cGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTCCGc                                                                                                                      >  1:552404/1‑151 (MQ=11)                            gCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGTGCATAGTCACCTTTCCGcc                                                                                                                     >  1:2302800/1‑151 (MQ=11)                            gCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTCCGcc                                                                                                                     >  1:524568/1‑151 (MQ=11)                            gCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTCCGcc                                                                                                                     >  2:5843614/1‑151 (MQ=11)                             ccATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTCCGccc                                                                                                                    <  2:4803367/151‑1 (MQ=11)                              cATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTCCGCCCt                                                                                                                   >  2:6481485/1‑151 (MQ=11)                              cATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTCCGCCCt                                                                                                                   >  1:1537166/1‑151 (MQ=11)                              cATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTCCGCCCt                                                                                                                   <  1:2997509/151‑1 (MQ=11)                               aTCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTCCGCCCtt                                                                                                                  <  1:6419024/151‑1 (MQ=11)                               aTCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTC‑‑‑‑‑‑‑GCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTTCTTGCTCCTTGCTT‑‑‑‑‑‑‑CTTGGATGTGCGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTCGCGCATAGTAACCTTTCCCCCCtt                                                                                                                  >  2:3816966/1‑151 (MQ=11)                                  gCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTCCGCCCTTAcg                                                                                                               <  2:1119856/151‑1 (MQ=11)                                   cAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTCCGCCCTTAcgg                                                                                                              >  1:689441/1‑151 (MQ=11)                                   cAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTCCGCCCTTAcgg                                                                                                              <  1:4477192/151‑1 (MQ=11)                                   cAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTCCGCCCTTAcgg                                                                                                              <  1:333655/151‑1 (MQ=11)                                     aaTGAACTGTTCACTGAGCCAAAACAGACGACTTC‑‑‑‑‑‑‑GCTTCTTGCGTCTTGCTTCATACTTATTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCAATCGCGCAGGTGGCGCATAGTCCCCTTTCCGCCCTTAcggcg                                                                                                            <  2:2304970/151‑1 (MQ=11)                                        ggACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTCCGCCCTTACGGCGGGt                                                                                                         <  2:2720154/151‑1 (MQ=11)                                        ggACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTCCGCCCTTACGGCGGGt                                                                                                         <  1:6091095/151‑1 (MQ=11)                                        ggACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTCCGCCCTTACGGCGGGt                                                                                                         >  1:1116213/1‑151 (MQ=11)                                         gACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTCCGCCCTTACGGCGGGTc                                                                                                        <  1:1056206/151‑1 (MQ=1)                                         gACTGTTCACTGAGCCAAAACAGACGAGTTCGCTTCTTGCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCGCGCAGGTGGCGCATAGTCACCTTTCCGCCCTTACGGCGGGTc                                                                                                        >  1:5323102/1‑151 (MQ=11)                                                                                                                                    cttgGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAAGTGGCGCATAGTCACCTTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACTAAAAAACGCTG‑GCTCCTATCATCCGGCTCCtgcgctgcgc                                 >  1:831187/1‑151 (MQ=255)                                                                                                                                                                   ttcattcaCCCAAGTGGCGCATAGTCACCTTTCCGCCCGTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTG‑GCTCCTATCATCCGGCCCCTGCGCTTCGCTCCGGGGTCCCCTCGCCACGGGCTGGCTcac  >  1:3914216/1‑149 (MQ=11)                                                                                                                                                                          |                                                                                                                                                 GTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAATGGACTGTTCACTGAGCCAAAACAGACGAGTTC‑‑‑‑‑‑‑GCTTCTTGCTTCTTGCTTC‑‑‑‑‑‑‑TTGCTT‑‑‑‑‑‑‑CTTGCTT‑‑‑‑‑‑‑CTTGGATGTGGGCATTCACATCGAGTAGGTATTCATTCACGCAAGTGGCGCATAGTCACCTTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCC  >  CP009974/4147063‑4147348

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑