BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP009974	2,738,507	G→A	E240E (GAG→GAA)	RPPX_12165 →	benzene 1,2‑dioxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009974	2,738,507	0	G	A	99.5%	423.6 / NA	193	E240E (GAG→GAA)	RPPX_12165	benzene 1,2‑dioxygenase
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (103/89); total (103/90)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 193 total aligned reads displayed.

CGGTGCCAGCTCGTACATCTACGAAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCGGGTGACGAGATCAAGACCATGAGTGCTGGCGCCTGGGCCAAGCAGGGCGGCGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTC  >  CP009974/2738369‑2738654

tggagctaGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATg                                                                                                                                         >  1:359946/8‑151 (MQ=1)

  gagctaGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAg                                                                                                                                       <  1:2083604/146‑1 (MQ=1)

  gagctaGCTCCTACATCGACGAAGGTGACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAg                                                                                                                                       <  2:110436/146‑1 (MQ=1)

   agctaGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAgc                                                                                                                                      <  1:6709857/147‑2 (MQ=1)

    gctaGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAgcg                                                                                                                                     >  2:3047589/4‑149 (MQ=1)

    gctaGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAgcg                                                                                                                                     >  1:2869055/4‑149 (MQ=1)

                                              gagaaCGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCg                                                                                           >  1:3695646/4‑151 (MQ=0)

                                              gagaaCGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCg                                                                                           >  1:5920433/4‑151 (MQ=0)

                                              gagaaCGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCg                                                                                           >  1:6030034/4‑151 (MQ=0)

                                                   cGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCAt                                                                                      >  2:3785492/1‑151 (MQ=2)

                                                    ggCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATg                                                                                     >  2:1717807/1‑151 (MQ=2)

                                                    ggCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATg                                                                                     >  2:4213504/1‑151 (MQ=11)

                                                    ggCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATg                                                                                     <  2:7021841/151‑1 (MQ=2)

                                                    ggCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGGTTCTATTCCTTCGACCATg                                                                                     >  1:2965556/1‑151 (MQ=11)

                                                     gcgcCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATgg                                                                                    <  2:1822900/151‑1 (MQ=2)

                                                         cGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGCTTCTATTCCTTCGACCATGGacac                                                                                >  2:503022/1‑151 (MQ=2)

                                                          gACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACAcc                                                                               >  2:2642154/1‑151 (MQ=11)

                                                          gACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACAcc                                                                               >  1:4315278/1‑151 (MQ=2)

                                                          gACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACAcc                                                                               >  2:1377682/1‑151 (MQ=11)

                                                           aCGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACct                                                                              <  1:1718050/151‑1 (MQ=11)

                                                            cGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACctg                                                                             >  1:5768689/1‑151 (MQ=2)

                                                                  tcatgttaGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGctc                                                                       <  2:7240359/144‑1 (MQ=2)

                                                                    atgttaGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCCCTg                                                                     >  2:5110652/6‑151 (MQ=11)

                                                                     tgttaGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTgg                                                                    <  1:3399385/147‑1 (MQ=11)

                                                                       ttaGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGAc                                                                  >  2:1711175/3‑151 (MQ=11)

                                                                        taGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACt                                                                 >  2:4928561/2‑151 (MQ=14)

                                                                         aGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTc                                                                <  1:4686935/151‑1 (MQ=14)

                                                                          gCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCg                                                               >  1:7056799/1‑151 (MQ=14)

                                                                          gCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCg                                                               >  1:1813370/1‑151 (MQ=14)

                                                                          gCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCg                                                               <  2:6355988/151‑1 (MQ=14)

                                                                          gCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCAGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCg                                                               <  1:3357047/151‑1 (MQ=14)

                                                                           cTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTCAGCAGGGAGGTGGTCTCTATTCCTTCGACCATGGACACCTGCTGCTCTTGACGCgg                                                              >  2:6621345/1‑151 (MQ=11)

                                                                             ccGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTg                                                            >  1:5644927/1‑151 (MQ=14)

                                                                               gTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTggg                                                          <  2:2694895/151‑1 (MQ=14)

                                                                               gTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTggg                                                          >  2:5405172/1‑151 (MQ=14)

                                                                               gTGCACAGGCACTACGCCGCCACGCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGCACACCTGCTGCTCTGGACTCGGTggg                                                          <  1:983503/151‑1 (MQ=14)

                                                                                tGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGc                                                         >  2:5715830/1‑151 (MQ=14)

                                                                                tGCACTGGAACGACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGc                                                         <  2:3313660/151‑1 (MQ=14)

                                                                                 gCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCg                                                        <  2:911594/151‑1 (MQ=14)

                                                                                  cACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGa                                                       >  1:680383/1‑151 (MQ=14)

                                                                                   aCTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGaa                                                      <  2:2261322/151‑1 (MQ=14)

                                                                                   aCTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCGAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGaa                                                      >  2:735129/1‑151 (MQ=14)

                                                                                      ggAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGGCGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAAccc                                                   <  2:4752381/151‑1 (MQ=14)

                                                                                       ggaCTGCGCCGCCCCTCGGAACCAGCGCGAGCAGCGGGAGGCAGGCGACGAGATCAAGACCATGAGCGCCGGCTCTTCGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAAcccc                                                  <  2:6573150/149‑1 (MQ=255)

                                                                                         aCTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTCAGCAGGGTGGAGGTTTCGATTCCTTAGACCATGGACACCTGCTGCTCTGGGCGCGGTGGGCGAAccccaa                                                >  2:4338100/1‑149 (MQ=14)

                                                                                           tACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAgg                                              >  2:5456609/1‑151 (MQ=17)

                                                                                           tACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAgg                                              <  2:2453278/151‑1 (MQ=17)

                                                                                           tACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAgg                                              <  2:471/151‑1 (MQ=17)

                                                                                           tACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAgg                                              <  1:7007762/151‑1 (MQ=17)

                                                                                           tACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAgg                                              >  1:6419306/1‑151 (MQ=17)

                                                                                            aCGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGa                                             >  1:4643802/1‑151 (MQ=17)

                                                                                            aCGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGa                                             >  1:1521046/1‑151 (MQ=17)

                                                                                              gccgCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATc                                           <  2:5735507/148‑1 (MQ=17)

                                                                                              gccgCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATc                                           <  2:1455815/148‑1 (MQ=17)

                                                                                                  ccACTCAGAACCAGCGCAAGCAGCGCGCGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGccc                                       <  1:2462756/151‑1 (MQ=17)

                                                                                                   cACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCa                                      >  1:2267730/1‑151 (MQ=17)

                                                                                                   cACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCa                                      >  2:1299645/1‑151 (MQ=17)

                                                                                                    actcAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAg                                     >  1:2962158/4‑151 (MQ=18)

                                                                                                     ctcAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGc                                    >  1:7190507/3‑151 (MQ=17)

                                                                                                      tcAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGTACTCGGTGGGCGAACCCCGAGGATCGCGCAGcc                                   >  1:3304194/2‑150 (MQ=17)

                                                                                                      gcAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCt                                   <  2:5199194/150‑1 (MQ=18)

                                                                                                       cAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAACCCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCtt                                  <  2:72308/151‑1 (MQ=18)

                                                                                                       cAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCtt                                  >  2:1580416/1‑151 (MQ=18)

                                                                                                        aGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTa                                 >  2:6785363/1‑151 (MQ=18)

                                                                                                         gAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTAt                                >  2:5591999/1‑151 (MQ=18)

                                                                                                          aaCCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATg                               >  1:3272007/1‑151 (MQ=18)

                                                                                                           aCCAGCGCAAGCAGCGCGAGGCTGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTTTATTCCTTCGGCCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGa                              <  2:2158292/151‑1 (MQ=255)

                                                                                                           aCCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGa                              >  1:4769679/1‑151 (MQ=18)

                                                                                                           aCCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGa                              >  1:4056280/1‑151 (MQ=18)

                                                                                                             cAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAgc                            >  1:5073036/1‑151 (MQ=18)

                                                                                                               gcgcAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAgcgc                          <  2:2007705/151‑1 (MQ=255)

                                                                                                                cgcAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGcc                         <  1:6905109/151‑1 (MQ=255)

                                                                                                                cgcAAGCAGCGCGAGGCAGGCGACGAAAGCAAAACCATGAGCGCCGGCGCTTGGGCTAGGCAGGGTGGGGGTTTCTATTCCTTCGACCATGGACGCCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGATCGcc                         <  2:5962718/151‑1 (MQ=255)

                                                                                                                 gcAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCg                        >  1:6705299/1‑151 (MQ=18)

                                                                                                                  cAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGt                       <  1:4999594/151‑1 (MQ=18)

                                                                                                                   aaGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTg                      >  1:144738/1‑151 (MQ=255)

                                                                                                                    agcagcGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGa                     >  2:5880963/1‑151 (MQ=255)

                                                                                                                     gcagcGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGAc                    >  2:1632269/1‑151 (MQ=18)

                                                                                                                     gcagcGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGAc                    <  1:1397013/151‑1 (MQ=255)

                                                                                                                      cagcGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGAcc                   <  2:1719832/151‑1 (MQ=21)

                                                                                                                      cagcGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGAcc                   <  1:1929714/151‑1 (MQ=18)

                                                                                                                      cagcGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCGCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGAcc                   <  2:315575/151‑1 (MQ=18)

                                                                                                                       ggcgcGAGGCGGGCGACGAAATCAGAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCa                  <  2:3246749/150‑1 (MQ=255)

                                                                                                                       agcGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCa                  <  1:7002734/151‑1 (MQ=255)

                                                                                                                       agcGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCa                  >  2:3759549/1‑151 (MQ=255)

                                                                                                                       agcGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCa                  >  1:800493/1‑151 (MQ=255)

                                                                                                                         cgcgAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAAt                <  1:10098/151‑1 (MQ=255)

                                                                                                                          gcgAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCCGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCACCTTATGAGCGCCGTGACCactc               >  2:6423492/1‑148 (MQ=255)

                                                                                                                           cgAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTg              <  1:1181329/151‑1 (MQ=255)

                                                                                                                            gAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTgg             <  2:1008058/151‑1 (MQ=255)

                                                                                                                            gAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTgg             >  1:5434409/1‑151 (MQ=255)

                                                                                                                            gAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTgg             >  1:1624724/1‑151 (MQ=255)

                                                                                                                             aggcaggcGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGc            >  1:5099731/1‑151 (MQ=21)

                                                                                                                               gcaggcgACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTc          >  2:4983638/7‑151 (MQ=255)

                                                                                                                                caggcgACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCa         >  1:2910160/6‑151 (MQ=255)

                                                                                                                                  ggcgACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAg       >  2:169349/4‑151 (MQ=255)

                                                                                                                                   gcgACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGCGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGa      >  1:7169408/3‑151 (MQ=255)

                                                                                                                                      aCGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGAttt   >  2:6246268/1‑151 (MQ=255)

                                                                                                                                       cGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTc  >  1:4965360/1‑151 (MQ=255)

                                                                                                                                       cGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTc  >  1:3761527/1‑151 (MQ=255)

CGGTGCCAGCTCGTACATCTACGAAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCGGGTGACGAGATCAAGACCATGAGTGCTGGCGCCTGGGCCAAGCAGGGCGGCGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTC  >  CP009974/2738369‑2738654

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