BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,737,978	G→A	5.8%	G482G (GGC→GGT)	hutU ←	urocanate hydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,737,978	0	G	A	5.8%	875.0 / 27.0	276	G482G (GGC→GGT)	hutU	urocanate hydratase
Reads supporting (aligned to +/- strand): ref base G (111/149); new base A (7/9); total (118/158)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 302 total aligned reads displayed.

GCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTC  >  NC_002947/5737830‑5738127

gCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGgcgc                                                                                                                                                      >  2:8611656/1‑150 (MQ=255)

gCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGgcgc                                                                                                                                                      >  1:5236799/1‑150 (MQ=255)

gCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGgcgc                                                                                                                                                      >  1:8215390/1‑150 (MQ=255)

 cATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGcgcc                                                                                                                                                     <  1:8225311/150‑1 (MQ=255)

 cATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGcgcc                                                                                                                                                     <  1:3264682/150‑1 (MQ=255)

     aCGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCt                                                                                                                                                 <  1:5062659/150‑1 (MQ=255)

       gCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGc                                                                                                                                               <  1:734745/150‑1 (MQ=255)

       gCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGc                                                                                                                                               >  1:6321134/1‑150 (MQ=255)

       gCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGCGGCGCCACCTgg                                                                                                                                               <  1:2102979/150‑2 (MQ=255)

        ccAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCg                                                                                                                                              <  2:335895/150‑1 (MQ=255)

          aGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGgt                                                                                                                                            >  1:6766879/1‑150 (MQ=255)

          aGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGgt                                                                                                                                            <  2:3468979/150‑1 (MQ=255)

             cccTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGAGGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGtt                                                                                                                                         <  2:8914767/150‑1 (MQ=255)

              ccTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTc                                                                                                                                        <  1:1614689/150‑1 (MQ=255)

                 gggTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTcagc                                                                                                                                     <  1:4837222/150‑1 (MQ=255)

                  ggTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTcagca                                                                                                                                    <  2:543720/150‑1 (MQ=255)

                    tCGTTGGTCTGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAgg                                                                                                                                  <  2:1922227/150‑1 (MQ=255)

                    tCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAgg                                                                                                                                  <  2:2329525/150‑1 (MQ=255)

                     cGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAggg                                                                                                                                 <  1:2257951/150‑1 (MQ=255)

                      gTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGc                                                                                                                                <  2:3894310/150‑1 (MQ=255)

                      gTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGc                                                                                                                                >  1:3039777/1‑150 (MQ=255)

                          gTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTg                                                                                                                            <  1:2564506/150‑1 (MQ=255)

                          gTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTg                                                                                                                            <  1:5579122/150‑1 (MQ=255)

                           tCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCGTGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGa                                                                                                                           <  2:6134304/150‑1 (MQ=255)

                           tCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGa                                                                                                                           >  1:7657476/1‑150 (MQ=255)

                           tCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGa                                                                                                                           <  2:1148111/150‑1 (MQ=255)

                           tCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGa                                                                                                                           <  2:1323604/150‑1 (MQ=255)

                           tCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGa                                                                                                                           <  2:80401/150‑1 (MQ=255)

                           tCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGa                                                                                                                           <  2:1689175/150‑1 (MQ=255)

                           tCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGa                                                                                                                           <  2:7934016/150‑1 (MQ=255)

                                aCCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGagcagc                                                                                                                      <  2:8670348/150‑1 (MQ=255)

                                aCCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGagcagc                                                                                                                      <  2:7328920/150‑1 (MQ=255)

                                 cccGGGCGATGCGCTCGGCGGCCTCATCGGTCCCGTCGCAGCCGAGGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCAGGGTGCAGCGATACCCCGGTGGCGCCGCCGGCGGTGTTCAGCAGGGCGTTGAGCAGCg                                                                                                                     >  1:5246045/1‑150 (MQ=255)

                                  ccGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCgg                                                                                                                    <  2:267048/150‑1 (MQ=255)

                                  ccGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCgg                                                                                                                    <  2:2505344/150‑1 (MQ=255)

                                   cGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGc                                                                                                                   >  1:6207454/1‑150 (MQ=255)

                                     ggCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCa                                                                                                                 <  2:4783010/150‑1 (MQ=255)

                                      gCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGtt                                                                                                                                         <  1:5750213/125‑1 (MQ=255)

                                      gCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGtt                                                                                                                                         >  2:5750213/1‑125 (MQ=255)

                                        gATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTc                                                                                                              <  2:4719174/150‑1 (MQ=255)

                                        gATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTc                                                                                                              >  2:300479/1‑150 (MQ=255)

                                            cgcTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGaa                                                                                                          <  1:6625023/150‑1 (MQ=255)

                                             gcTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGaaa                                                                                                         <  2:7051159/150‑1 (MQ=255)

                                             gcTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGaaa                                                                                                         >  1:7975173/1‑150 (MQ=255)

                                             gcTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGaaa                                                                                                         <  2:3289090/150‑1 (MQ=255)

                                             gcTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGaaa                                                                                                         >  1:3796765/1‑150 (MQ=255)

                                             gcTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGaaa                                                                                                         <  2:3805655/150‑1 (MQ=255)

                                               tCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACa                                                                                                       >  1:4575131/1‑150 (MQ=255)

                                                  tcggcATCTTCCGAACCGTCGCAGCAGATGACCATGCCCGAATGTTGGGAGAAGCCCATGCCGACGCCGCCACCATGATGGAGCGATACCCATGTCGCACCCGATGCGGTGTTGAGCAACGCATTGAGCAGCGGCCAGTCGGAGACTGCg                                                                                                    >  2:2892170/2‑150 (MQ=255)

                                                    ggcCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGCTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCggc                                                                                                  >  2:3824794/1‑148 (MQ=255)

                                                       cTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGa                                                                                               >  1:2774908/1‑150 (MQ=255)

                                                        tCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGaa                                                                                              >  1:8855216/1‑150 (MQ=255)

                                                          ttccgaaCCGTCGCAGCAGATGACCATGCCCGAATGTTGGGAGAAGCCCATGCCGACGCCGCCACCATGATGGAGCGATACCCATGTCGCACCCGATGCGGTGTTGAGCAACGCATTGAGCAGCGGCCAGTCGGAGACTGCGTCCGAGcc                                                                                            >  1:1125367/7‑150 (MQ=255)

                                                             ggTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATc                                                                                         <  2:8776183/150‑1 (MQ=255)

                                                             ggTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATc                                                                                         <  1:4137359/150‑1 (MQ=255)

                                                              gaaCCGTCGCAGCAGATGACCATGCCCGAATGTTGGGAGAAGCCCATGCCGACGCCGCCACCATGATGGAGCGATACCCATGTCGCACCCGATGCGGTGTTGAGCAACGCATTGAGCAGCGGCCAGTCGGAGACTGCGTCCGAGCCGTct                                                                                        >  2:9014121/3‑149 (MQ=255)

                                                                   gTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCAt                                                                                   <  2:1799913/150‑1 (MQ=255)

                                                                    tCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATg                                                                                  <  1:3588422/150‑1 (MQ=255)

                                                                       cAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGcc                                                                               >  1:6063619/1‑150 (MQ=255)

                                                                         gACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTc                                                                             >  1:6716450/1‑150 (MQ=255)

                                                                             aTGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGtt                                                                         <  1:5414446/150‑1 (MQ=255)

                                                                              tGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGttt                                                                        >  2:4844578/1‑150 (MQ=255)

                                                                                 ccATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTcgc                                                                     <  2:5070405/150‑1 (MQ=255)

                                                                                  cATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTcgcg                                                                    >  1:8575759/1‑150 (MQ=255)

                                                                                   aTGCCCGAATGTTGGGAGAAGCCCATACCGACGCCGCCACCATGATGGAGCGATACCCATGTCGCACCCGATGCGGTGTTGAGCAACGCATTGAGCAGCGGCCAGTCGGAGACTGCGTCCGAGCCGTCTTTCATGGATTCGGTTTCACgg                                                                   <  2:6760592/150‑1 (MQ=255)

                                                                                       ccGAATGTTGGGAGAAGCCCATGCCGACGCCGCCACCATGATGGAGCGATACCCATGTCGCACCCGATGCGGTGTTGAGCAACGCATTGAGCAGCGGCCAGTCGGAGACTGCGTCCGAGCCGTCTTTCATGGATTCGGTTTCACGGtttg                                                               <  1:4253139/150‑3 (MQ=255)

                                                                                         gAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCggg                                                             <  2:3532745/150‑1 (MQ=255)

                                                                                           gTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCt                                                           >  2:1997641/1‑150 (MQ=255)

                                                                                               tgcGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGAt                                                       >  1:2649802/1‑150 (MQ=255)

                                                                                               tgcGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGAt                                                       >  1:1575028/1‑150 (MQ=255)

                                                                                                gcGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATa                                                      >  2:6278932/1‑150 (MQ=255)

                                                                                                 cGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATAc                                                     >  2:3273027/1‑150 (MQ=255)

                                                                                                 cGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATAc                                                     <  1:3698361/150‑1 (MQ=255)

                                                                                                 cGAGAAGCCAATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATAc                                                     <  1:2139568/150‑1 (MQ=255)

                                                                                                   agaAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCg                                                   >  1:2986664/1‑150 (MQ=255)

                                                                                                    gaAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGTGGTTCGGGCTGGATACCGa                                                  >  2:5624239/1‑150 (MQ=255)

                                                                                                    gaAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGa                                                  >  2:1599766/1‑150 (MQ=255)

                                                                                                       gCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGcc                                               <  1:3422476/150‑1 (MQ=255)

                                                                                                        cccATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCg                                              >  2:8298542/1‑150 (MQ=255)

                                                                                                           aTGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAg                                           >  1:7801196/1‑150 (MQ=255)

                                                                                                           aTGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAg                                           <  2:2711990/150‑1 (MQ=255)

                                                                                                               ccACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGa                                       >  1:6405057/1‑150 (MQ=255)

                                                                                                                cACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAg                                      <  2:5846422/150‑1 (MQ=255)

                                                                                                                    ccACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAggtgg                                  >  1:7008654/1‑150 (MQ=255)

                                                                                                                    ccACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAggtgg                                  <  1:7557249/150‑1 (MQ=255)

                                                                                                                      aCCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTc                                >  2:5220116/1‑150 (MQ=255)

                                                                                                                        cGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCAc                              >  2:6896218/1‑150 (MQ=255)

                                                                                                                           cATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGc                           >  1:4167685/1‑150 (MQ=255)

                                                                                                                               gTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGAt                       <  1:907215/150‑1 (MQ=255)

                                                                                                                                     cGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATcaccac                 >  2:4296333/1‑150 (MQ=255)

                                                                                                                                           ccAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTc                                <  2:7553899/129‑1 (MQ=255)

                                                                                                                                           ccAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGc           >  2:4116355/1‑150 (MQ=255)

                                                                                                                                           ccAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGc           >  1:3995982/1‑150 (MQ=255)

                                                                                                                                             aGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCg         >  1:1556159/1‑150 (MQ=255)

                                                                                                                                              ggTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGa        <  2:3711725/150‑1 (MQ=255)

                                                                                                                                                tGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACa      <  2:6062384/150‑1 (MQ=255)

                                                                                                                                                tGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACa      <  1:5885824/150‑1 (MQ=255)

                                                                                                                                                    gccgccTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTc  >  2:7504058/1‑150 (MQ=255)

                                                                                                                                                    gccgccTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTc  <  2:3969308/150‑1 (MQ=255)

                                                                                                                                                    gccgccTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTc  >  2:4458289/1‑150 (MQ=255)

GCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTC  >  NC_002947/5737830‑5738127

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