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	100.0%	371.8 / NA	165	E240E (GAG→GAA)	RPPX_12165	benzene 1,2‑dioxygenase
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (70/95); total (70/95)

Only 100 of 165 total aligned reads displayed.

GGTGCCAGCTCGTACATCTACGAAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCGGGTGACGAGATCAAGACCATGAGTGCTGGCGCCTGGGCCAAGCAGGGCGGCGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTC  >  CP009974/2738370‑2738654

ggagctaGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGa                                                                                                                                        >  1:4090766/7‑151 (MQ=1)

  agctaGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAgc                                                                                                                                      >  1:1887559/5‑150 (MQ=1)

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

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

   gctaGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAgcg                                                                                                                                     <  1:172489/148‑3 (MQ=1)

   gctaGCTCCTACAACTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGAAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAgcg                                                                                                                                     <  1:5227194/148‑3 (MQ=1)

         tcctACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGgct                                                                                                                               <  1:4932243/148‑2 (MQ=0)

                    ggagggtaaCGGGAAGCTGACCGCCGAGAACGGGGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCCAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAg                                                                                                                    <  2:3519538/144‑1 (MQ=0)

                                             gagaaCGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCg                                                                                           <  1:4866971/148‑1 (MQ=0)

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

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

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

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

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

                                                     cgcCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGa                                                                                   <  2:427459/151‑1 (MQ=11)

                                                     cgcCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGa                                                                                   >  2:1301533/1‑151 (MQ=2)

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

                                                          aCGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACct                                                                              <  1:3640066/151‑1 (MQ=2)

                                                           cGGCTATCATGTTGGCTCCGTGCACTGGAACTACGCCGCCCCTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGCCGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTCAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACctg                                                                             <  2:1196736/151‑1 (MQ=2)

                                                           cGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACctg                                                                             <  2:1074243/151‑1 (MQ=2)

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

                                                           cGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACctg                                                                             <  1:1018359/151‑1 (MQ=2)

                                                           cGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACctg                                                                             <  1:3343224/151‑1 (MQ=2)

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

                                                                 tcatgttaGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGctc                                                                       >  1:5151866/8‑151 (MQ=11)

                                                                  catgttaGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGctct                                                                      >  2:27921/7‑151 (MQ=11)

                                                                      ttaGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGAc                                                                  <  2:2550562/149‑1 (MQ=14)

                                                                       taGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACt                                                                 <  2:951691/150‑1 (MQ=14)

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

                                                                           tCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGt                                                             <  1:4656216/151‑1 (MQ=14)

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

                                                                                gCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCg                                                        >  1:3322156/1‑151 (MQ=14)

                                                                                  aCTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTGAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGaa                                                      <  2:1656047/151‑1 (MQ=14)

                                                                                   cTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAAc                                                     >  1:1516255/1‑151 (MQ=17)

                                                                                      gAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAAcccc                                                  >  2:4255058/1‑151 (MQ=14)

                                                                                        aCTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGa                                                <  1:1609330/151‑1 (MQ=17)

                                                                                         cTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCGTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAg                                               <  1:3508599/151‑1 (MQ=17)

                                                                                           aCGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCTAACCCCGAGGa                                             >  1:4905252/1‑151 (MQ=14)

                                                                                           aCGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGa                                             <  2:1501916/151‑1 (MQ=17)

                                                                                            cgccgcCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGAt                                            <  2:77199/151‑1 (MQ=17)

                                                                                            cgccgcCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGAt                                            >  1:4991762/1‑151 (MQ=17)

                                                                                             gccgCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGTGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATc                                           >  1:2275347/4‑151 (MQ=17)

                                                                                             gccgCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATc                                           >  1:3026126/4‑151 (MQ=17)

                                                                                             gccgCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATc                                           >  2:1075860/4‑151 (MQ=17)

                                                                                              ccgCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAACACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCg                                          <  2:2704931/149‑1 (MQ=17)

                                                                                               cgCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGc                                         >  1:4060049/2‑151 (MQ=17)

                                                                                               cgCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGc                                         <  2:1442260/150‑1 (MQ=17)

                                                                                                gCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGcc                                        <  1:3212642/151‑1 (MQ=18)

                                                                                                gCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGcc                                        >  1:5229121/1‑151 (MQ=17)

                                                                                                 ctactcAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGGCTCGGTGGGCGGACCCCGAGGATCGccc                                       >  2:2990509/6‑151 (MQ=17)

                                                                                                  cACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCa                                      <  1:1008482/151‑1 (MQ=18)

                                                                                                  cACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGACGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCa                                      <  2:1887583/151‑1 (MQ=18)

                                                                                                      cAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCACGAGGATCGCCCCGCtt                                  >  2:4828671/1‑151 (MQ=17)

                                                                                                        gAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTAt                                <  1:2803559/151‑1 (MQ=18)

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

                                                                                                         aaCCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATg                               <  1:2416612/151‑1 (MQ=18)

                                                                                                          aCCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGa                              <  2:4090766/151‑1 (MQ=255)

                                                                                                          aCCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGa                              <  2:1381787/151‑1 (MQ=255)

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

                                                                                                          aCCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGa                              <  1:856191/151‑1 (MQ=255)

                                                                                                          aCCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGCTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGa                              <  2:1597764/151‑1 (MQ=255)

                                                                                                           ccAGGGCAAGCAGCGCGGGGCAGGCGACGAAATCAAAACCAGGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTAGTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAg                             <  2:2270369/151‑1 (MQ=255)

                                                                                                           ccAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAg                             >  2:1215765/1‑151 (MQ=18)

                                                                                                           ccAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAg                             >  1:4678848/1‑151 (MQ=18)

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

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

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

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

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

                                                                                                                  aaGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTg                      <  1:3005224/151‑1 (MQ=255)

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

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

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

                                                                                                                   agcagcGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGa                     <  1:2884226/151‑1 (MQ=255)

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

                                                                                                                    gcagcGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGCACACCTGCTGCTCGGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGAc                    >  2:5014632/1‑151 (MQ=18)

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

                                                                                                                       gcgcGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCaa                 <  1:2323837/151‑1 (MQ=255)

                                                                                                                       gcgcGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCaa                 <  1:4779026/151‑1 (MQ=255)

                                                                                                                          cgAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTg              >  1:2722232/1‑151 (MQ=255)

                                                                                                                          cgAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCCAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTg              >  1:2996871/1‑151 (MQ=255)

                                                                                                                           caggcaggcGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTgg             <  1:4279759/150‑1 (MQ=255)

                                                                                                                            aggcaggcGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGCCCAATTGGc            >  1:1330315/1‑151 (MQ=18)

                                                                                                                             ggcaggcgACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCt           <  1:4286885/144‑1 (MQ=21)

                                                                                                                             ggcaggcgACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCt           >  1:619621/8‑151 (MQ=255)

                                                                                                                               caggcgACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGGTCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCa         >  2:3258392/6‑151 (MQ=255)

                                                                                                                                aggcgACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGAt                                            <  1:1211345/111‑1 (MQ=255)

                                                                                                                                aggcgACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGAt                                            >  2:1211345/5‑115 (MQ=255)

                                                                                                                                aggcgACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCaa        >  2:2392399/5‑151 (MQ=21)

                                                                                                                                aggcgACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCaa        >  1:1169797/5‑151 (MQ=21)

                                                                                                                                 ggcgACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAg       <  1:1812855/148‑1 (MQ=21)

                                                                                                                                  gcgACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGa      <  1:4486157/149‑1 (MQ=255)

                                                                                                                                   cgACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGAt     <  2:74041/150‑1 (MQ=255)

                                                                                                                                   cgACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGAt     <  2:260480/150‑1 (MQ=255)

                                                                                                                                   cgACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGAt     <  1:2649930/150‑1 (MQ=255)

                                                                                                                                   cgACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGAt     >  2:389126/2‑151 (MQ=255)

                                                                                                                                   cgACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGAt     >  2:3593210/2‑151 (MQ=255)

                                                                                                                                     aCGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAACTGGCTCAAGAttt   >  1:4495904/1‑151 (MQ=255)

                                                                                                                                     aCGAAACCAAAACCATGAGCGCCGGCTCTTGGCCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGAttt   <  2:3882279/151‑1 (MQ=21)

                                                                                                                                      cGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTc  <  1:3630754/151‑1 (MQ=255)

GGTGCCAGCTCGTACATCTACGAAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCGGGTGACGAGATCAAGACCATGAGTGCTGGCGCCTGGGCCAAGCAGGGCGGCGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTC  >  CP009974/2738370‑2738654

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