BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,992,345	C→T	11.7%	A100A (GCG→GCA)	ycaC‑III ←	putative hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,992,345	0	C	T	11.7%	610.7 / 50.8	207	A100A (GCG→GCA)	ycaC‑III	putative hydrolase
Reads supporting (aligned to +/- strand): ref base C (107/74); new base T (11/13); total (119/88)
Fisher's exact test for biased strand distribution p-value = 2.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01

Only 100 of 211 total aligned reads displayed.

TCACGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAA  >  NC_002947/5992196‑5992494

tCACGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCAcc                                                                                                                                                       >  2:1640659/1‑150 (MQ=255)

    gggTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCt                                                                                                                                                   <  2:574332/150‑1 (MQ=255)

     ggTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCtt                                                                                                                                                  <  1:4355885/150‑1 (MQ=255)

        catttcattGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGAc                                                                                                                                               <  1:1113664/150‑1 (MQ=255)

          tttcattGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGa                                                                                                                                             >  2:1686847/1‑150 (MQ=255)

          tttcattGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGa                                                                                                                                             >  2:1667887/1‑150 (MQ=255)

          tttcattGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGa                                                                                                                                             >  1:4138956/1‑150 (MQ=255)

          tttcattGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGa                                                                                                                                             >  2:5113077/1‑150 (MQ=255)

              attGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATc                                                                                                                                         >  1:2515799/1‑150 (MQ=255)

               ttGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATcc                                                                                                                                        <  1:2808228/150‑1 (MQ=255)

               ttGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATcc                                                                                                                                        >  1:2609349/1‑150 (MQ=255)

               ttGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATcc                                                                                                                                        >  1:70695/1‑150 (MQ=255)

               ttGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTGCTTGCCGGTCGCCTTCACCGCCTTGACGaagccc                                                                                                                                        >  2:3712819/1‑146 (MQ=255)

                 gAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTc                                                                                                                                      <  2:2412027/150‑1 (MQ=255)

                  aaGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCg                                                                                                                                     >  1:2042226/1‑150 (MQ=255)

                  aaGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCg                                                                                                                                     >  2:1712688/1‑150 (MQ=255)

                     gTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTg                                                                                                                                  >  1:3556120/1‑150 (MQ=255)

                     gTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTg                                                                                                                                  >  1:3318519/1‑150 (MQ=255)

                     gTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTg                                                                                                                                  >  2:3627082/1‑150 (MQ=255)

                         cGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTccc                                                                                                                              <  1:2871071/150‑1 (MQ=255)

                           gAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCATCGCCTTGACGAAATCCTCGTTGTCCCAg                                                                                                                            <  2:2630553/150‑1 (MQ=255)

                           gAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAg                                                                                                                            >  2:1187953/1‑150 (MQ=255)

                                aTCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCCAGGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCggg                                                                                                                       >  2:4287184/1‑148 (MQ=255)

                                aTCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGtt                                                                                                                       >  1:3397109/1‑150 (MQ=255)

                                aTCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGtt                                                                                                                       >  1:2899247/1‑150 (MQ=255)

                                      caccacGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTg                                                                                                                 <  2:726579/150‑1 (MQ=255)

                                      caccacGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTg                                                                                                                 >  1:3172886/1‑150 (MQ=255)

                                          acGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCa                                                                                                             <  1:661103/150‑1 (MQ=255)

                                                acCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAgggcgg                                                                                                       <  2:3423474/150‑1 (MQ=255)

                                                  cTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAgggcgggc                                                                                                     >  1:3109164/1‑150 (MQ=255)

                                                  cTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAACGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAgggcgggc                                                                                                     >  1:4261693/1‑150 (MQ=255)

                                                   tCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCa                                                                                                    >  1:1858969/1‑150 (MQ=255)

                                                   tCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCa                                                                                                    >  1:623630/1‑150 (MQ=255)

                                                   tCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCa                                                                                                    >  2:2483376/1‑150 (MQ=255)

                                                   tCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCa                                                                                                    >  1:3749577/1‑150 (MQ=255)

                                                     aaaTTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAAt                                                                                                  >  2:4921796/1‑150 (MQ=255)

                                                        ttcttcCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTa                                                                                               <  1:2249904/150‑1 (MQ=255)

                                                         tcttcCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTAc                                                                                              <  1:4272810/150‑1 (MQ=255)

                                                           ttcCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACgg                                                                                            >  1:4827129/1‑150 (MQ=255)

                                                            tcctccAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGgc                                                                                           <  1:4681515/150‑1 (MQ=255)

                                                            tcctccAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGgc                                                                                           >  2:4964537/1‑150 (MQ=255)

                                                               tccAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCa                                                                                        >  2:1211879/1‑150 (MQ=255)

                                                               tccAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCa                                                                                        >  1:3083723/1‑150 (MQ=255)

                                                               tccAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCa                                                                                        >  2:2255613/1‑150 (MQ=255)

                                                                 cAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATc                                                                                      >  1:2278039/1‑150 (MQ=255)

                                                                  aaGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATcc                                                                                     >  2:4522841/1‑150 (MQ=255)

                                                                   aGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCg                                                                                    <  2:4824123/150‑1 (MQ=255)

                                                                       ggacaatgcCGGGAAAGCCACGCAGACTTCGGTCACAACACCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGaa                                                                                >  1:133503/8‑150 (MQ=14)

                                                                        gACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAAc                                                                               >  1:461573/1‑150 (MQ=255)

                                                                              gcCGGGAAAGCCACGCAGACTTCGGTCACAACACCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGaatgtttct                                                                         <  1:4916040/150‑8 (MQ=14)

                                                                                cGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCttt                                                                       <  2:4001876/150‑1 (MQ=255)

                                                                                  ggAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGACCAGGGCTTTCa                                                                     >  2:565909/1‑150 (MQ=255)

                                                                                       gCCCCGCTAACCTCGGTCACATAACCTGCAATGATGAGCTGCTTCTTGCAGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGGGTGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTcc                                                                <  2:1132587/150‑1 (MQ=18)

                                                                                       gCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTcc                                                                <  1:254494/150‑1 (MQ=255)

                                                                                       gCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTcc                                                                >  2:1384170/1‑150 (MQ=255)

                                                                                        ccACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCg                                                               >  2:188262/1‑150 (MQ=255)

                                                                                        ccACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCg                                                               <  1:3292405/150‑1 (MQ=255)

                                                                                            gCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCAc                                                           >  1:3659136/1‑150 (MQ=255)

                                                                                                aCCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGc                                                       >  2:4616406/1‑150 (MQ=255)

                                                                                                aCCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGc                                                       >  1:1212033/1‑150 (MQ=255)

                                                                                                   tCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTGGACGAAATCCTCGTTGTCCCAGGCGTGGATCGGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTATCAGGTCCGGCCCCAGCggg                                                    >  2:2395636/1‑150 (MQ=18)

                                                                                                        cacaACACCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGAATGTTTCTTTCAGCTCCGGCACCAATGGGCCAtt                                               <  1:1617095/150‑1 (MQ=9)

                                                                                                        cacaACACCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGAATGTTTCTTTCAGCTCCGGCACCAATGGGCCAtt                                               <  1:323039/150‑1 (MQ=9)

                                                                                                        cacaACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGtt                                               <  2:2828244/150‑1 (MQ=255)

                                                                                                         acaacaCCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCa                                                                                                             <  1:2742802/87‑1 (MQ=255)

                                                                                                         acaacaCCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCa                                                                                                             >  2:2742802/1‑87 (MQ=255)

                                                                                                         acaacaCCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTg                                              >  1:1454788/1‑150 (MQ=255)

                                                                                                           aacaCCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCCGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTggg                                            >  2:3559801/1‑150 (MQ=255)

                                                                                                            acacCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGAATGTTTCTTTCAGCTCCGGCACCAATGGGCCATTcggg                                           >  2:2697911/1‑150 (MQ=2)

                                                                                                             cacCTGAAATGATGAGCTGCTTCTTGACGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCTGGGCGGGCAATGTACGGCGCATCCGGGAAGAGGGCTTTCATTTCCGGCACCAGCGGGCCGTTGGGGc                                          <  1:610032/150‑1 (MQ=18)

                                                                                                              acCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGAATGTTTCTTTCAGCTCCGGCACCAATGGGCCATTCGGGcc                                         >  1:3855155/1‑150 (MQ=2)

                                                                                                              acCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGAATGTTTCTTTCAGCTCCGGCACCAATGGGCCATTCGGGcc                                         <  2:2286157/150‑1 (MQ=2)

                                                                                                                cTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGAATGTTTCTTTCAGCTCCGGCACCAATGGGCCATTCGGGccgt                                       >  2:683525/1‑148 (MQ=2)

                                                                                                                cTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCtt                                       <  2:2399510/150‑1 (MQ=255)

                                                                                                                 tGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTg                                      >  2:2884319/1‑150 (MQ=255)

                                                                                                                  gCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCCCCAACGGCCCGCTGGTGccggaa                                     <  1:3075910/150‑5 (MQ=17)

                                                                                                                  gCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGc                                     <  2:2004318/150‑1 (MQ=255)

                                                                                                                    aaTGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTc                                   >  2:1447982/1‑150 (MQ=255)

                                                                                                                      tGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGa                                 <  2:3000120/150‑1 (MQ=255)

                                                                                                                       gATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGaa                                <  1:313222/150‑1 (MQ=255)

                                                                                                                       gATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGaa                                <  2:621005/150‑1 (MQ=255)

                                                                                                                        aTCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAg                               >  1:322346/1‑150 (MQ=255)

                                                                                                                        aTCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAg                               >  1:5102251/1‑150 (MQ=255)

                                                                                                                              tgctTCTTGCCGGTCGCCTTCACCGCCTTGACGAACTCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCtggtg                         <  2:3786388/150‑1 (MQ=255)

                                                                                                                              tgctTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCtggtg                         <  2:3481206/150‑1 (MQ=255)

                                                                                                                              tgctTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCtggtg                         <  1:2634234/150‑1 (MQ=255)

                                                                                                                              tgctTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCtggtg                         <  2:2842764/150‑1 (MQ=255)

                                                                                                                                  tcttGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGa                     >  1:930801/1‑150 (MQ=255)

                                                                                                                                  tcttGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGa                     <  2:2338047/150‑1 (MQ=255)

                                                                                                                                   cttGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAg                    <  1:1689583/150‑1 (MQ=255)

                                                                                                                                             gCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGc          >  1:774222/1‑150 (MQ=255)

                                                                                                                                              ccTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGcc         >  1:4502564/1‑149 (MQ=255)

                                                                                                                                               cTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAg        >  2:3927089/1‑150 (MQ=255)

                                                                                                                                               cTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAg        >  1:1643772/1‑150 (MQ=255)

                                                                                                                                               cTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAg        >  1:1546497/1‑150 (MQ=255)

                                                                                                                                               cTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAg        >  2:293477/1‑150 (MQ=255)

                                                                                                                                                ttCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAgg       >  2:2098331/1‑150 (MQ=255)

                                                                                                                                                ttCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAgg       >  1:1925121/1‑150 (MQ=255)

                                                                                                                                                 tCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGt      >  1:2186863/1‑150 (MQ=255)

                                                                                                                                                    ccGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTga   >  2:2000345/1‑150 (MQ=255)

TCACGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAA  >  NC_002947/5992196‑5992494

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