BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,884,405	G→A	5.5%	K69K (AAG→AAA)	PP_1686 →	glutathione peroxidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,884,405	0	G	A	5.5%	1111.9 / 32.6	346	K69K (AAG→AAA)	PP_1686	glutathione peroxidase
Reads supporting (aligned to +/- strand): ref base G (138/189); new base A (14/5); total (152/194)
Fisher's exact test for biased strand distribution p-value = 8.63e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 353 total aligned reads displayed.

ATGCGTGCACATTGGTTGACGGTGCCGGTGCTGGCCCTGCTTGCCAGTACATCGAGCTGGGCTGCAGATTGCCCTGCGTTGTTGCAGGGCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTGGAAATGCTCGGCGTGCCGTCCAATGACTTCAAGCAGGAAGACGCCGACAGTGA  >  NC_002947/1884256‑1884554

 tGCGTGCACATTGGTTGACGGTGCCGGTGCTGGCCCTGCTTGCCAGTACATCGAGCTGGGCTGCAGATTGCCCTGCGTTGTTGCAGGGCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGc                                                                                                                                                      <  1:2198081/150‑1 (MQ=255)

    gTGCACATTGGTTGACGGTGCCGGTGCTGGCCCTGCTTGCCAGTACATCGAGCTGGGCTGCAGATTGCCCTGCGTTGTTGCAGGGCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGc                                                                                                                                                   <  1:3575330/150‑1 (MQ=255)

      gCACATTGGTTGACGGTGCCGGTGCTGGCCCTGCTTGCCAGTACATCGAGCTGGGCTGCAGATTGCCCTGCGTTGTTGCAGGGCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTc                                                                                                                                                 <  1:1216448/150‑1 (MQ=255)

            tGGTTGACGGTGCCGGTGCTGGCCCTGCTTGCCAGTACATCGAGCTGGGCTGCAGATTGCCCTGCGTTGTTGCAGGGCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCgtggtg                                                                                                                                           <  2:7662627/150‑1 (MQ=255)

            tGGTTGACGGTGCCGGTGCTGGCCCTGCTTGCCAGTACATCGAGCTGGGCTGCAGATTGCCCTGCGTTGTTGCAGGGCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCgtggtg                                                                                                                                           <  1:2882818/150‑1 (MQ=255)

             ggTTGACGGTGCCGGTGCTGGCCCTGCTTGCCAGTACATCGAGCTGGGCTGCAGATTGCCCTGCGTTGTTGCAGGGCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCgtggtgg                                                                                                                                          <  2:5892832/150‑1 (MQ=255)

              gTTGACGGTGCCGGTGCTGGCCCTGCTTGCCAGTACATCGAGCTGGGCTGCAGATTGCCCTGCGTTGTTGCAGGGCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGGCTTGCCGGCAAAGCGCTCgcag                                                                                                                                             >  2:1688207/1‑143 (MQ=255)

              gTTGACGGTGCCGGTGCTGGCCCTGCTTGCCAGTACATCGAGCTGAGCTGCAGATTGCCCTGCGTTGTTGCAGGGCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCgtggtggt                                                                                                                                         >  2:430743/1‑150 (MQ=255)

                tGACGGTGCCGGTGCTGGCCCTGCTTGCCAGTACATCGAGCTGGGCTGCAGATTGCCCTGCGTTGTTGCAGGGCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGa                                                                                                                                       >  2:703910/1‑150 (MQ=255)

                 gACGGTGCCGGTGCTGGCCCTGCTTGCCAGTACATCGAGCTGGGCTGCAGATTGCCCTGCGTTGTTGCAGGGCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGaa                                                                                                                                      <  1:5224931/150‑1 (MQ=255)

                  aCGGTGCCGGTGCTGGCCCTGCTTGCCAGTACATCGAGCTGGGCTGCAGATTGCCCTGCGTTGTTGCAGGGCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAac                                                                                                                                     >  1:5442170/1‑150 (MQ=255)

                  aCGGTGCCGGTGCTGGCCCTGCATGCCAGTACATCGAGCTGGGCTGCAGATTGCCCTGCGTTGTTGCAGGGCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAac                                                                                                                                     >  1:4160788/1‑150 (MQ=255)

                     gtgccggtgcTGGCCCTGCTTGCCAGTACATCGAGCTGGGCTGCAGATTGCCCTGCGTTGTTGCAGGGCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACt                                                                                                                                  <  1:206964/150‑1 (MQ=255)

                     gtgccggtgcTGGCCCTGCTTGCCAGTACATCGAGCTGGGCTGCAGATTGCCCTGCGTTGTTGCAGGGCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACt                                                                                                                                  >  2:3613725/1‑150 (MQ=255)

                      tgccggtgcTGGCCCTGCTTGCCAGTACATCGAGCTGGGCTGCAGATTGCCCTGCGTTGTTGCAGGGCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTg                                                                                                                                 >  2:7786046/1‑150 (MQ=255)

                         cggtgcTGGCCCTGCTTGCCAGTACATCGAGCTGGGCTGCAGATTGCCCTGCGTTGTTGCAGGGCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCa                                                                                                                              <  1:7145362/150‑1 (MQ=255)

                         cggtgcTGGCCCTGCTTGCCAGTACATCGAGCTGGGCTGCAGATTGCCCTGCGTTGTTGCAGGGCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCa                                                                                                                              <  1:8204099/150‑1 (MQ=255)

                           gtgcTGGCCCTGCTTGCCAGTACATCGAGCTGGGCTGCAGATTGCCCTGCGTTGTTGCAGGGCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGt                                                                                                                            >  1:6593800/1‑150 (MQ=255)

                                ggCCCTGCTTGCCAGTACATCGAGCTGGGCTGCAGATTGCCCTGCGTTGTTGCAGGGCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTg                                                                                                                       <  2:2648037/150‑1 (MQ=255)

                                    cTGCTTGCCAGTACATCGAGCTGGGCTGCAGATTGCCCTGCGTTGTTGCAGGGCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGt                                                                                                                   <  2:5333704/150‑1 (MQ=255)

                                     tGCTTGCCAGTACATCGAGCTGGGCTGCAGATTGCCCTGCGTTGTTGCAGGGCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGtt                                                                                                                  >  1:368212/1‑150 (MQ=255)

                                     tGCTTGCCAGTACATCGAGCTGGGCTGCAGATTGCCCTGCGTTGTTGCAGGGCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGtt                                                                                                                  >  2:7900302/1‑150 (MQ=255)

                                          gCCAGTACATCGAGCTGGGCTGCAGATTGCCCTGCGTTGTTGCAGGGCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGcc                                                                                                             >  1:2850525/1‑150 (MQ=255)

                                            cAGTACATCGAGCTGGGCTGCAGATTGCCCTGCGTTGTTGCAGGGCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGcccc                                                                                                           >  1:7599812/1‑150 (MQ=255)

                                             aGTACATCGAGCTGGGCTGCAGATTGCCCTGCGTTGTTGCAGGGCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCg                                                                                                          <  1:6169665/150‑1 (MQ=255)

                                                 cATCGAGCTGGGCTGCAGATTGCCCTGCGTTGTTGCAGGGCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGt                                                                                                      <  2:8388748/150‑1 (MQ=255)

                                                          gggCTGCAGATTGCCCTGCGTTGTTGCAGGGCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGc                                                                                             >  1:3063255/1‑150 (MQ=255)

                                                              tGCAGATTGCCCTGCGTTGTTGCAGGGCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGGCTTGCCGGCAAAGCGCTCgcag                                                                                                                                             <  2:5216887/98‑4 (MQ=38)

                                                              tGCAGATTGCCCTGCGTTGTTGCAGGGCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGa                                                                                         <  1:6211575/150‑1 (MQ=255)

                                                              tGCAGATTGCCCTGCGTTGTTGCAGGGCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGa                                                                                         <  1:116765/150‑1 (MQ=255)

                                                               gCAGATTGCCCTGCGTTGTTGCAGGGCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAg                                                                                        <  2:2920795/150‑1 (MQ=255)

                                                                cAGATTGCCCTGCGTTGTTGCAGGGCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAgg                                                                                       <  2:5104960/150‑1 (MQ=255)

                                                                  gATTGCCCTGCGTTGTTGCAGGGCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCg                                                                                     >  2:1857217/1‑150 (MQ=255)

                                                                        ccTGCGTTGTTGCAGGGCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTAc                                                                               >  2:8002984/1‑150 (MQ=255)

                                                                        ccTGCGTTGTTGCAGGGCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTAc                                                                               <  1:6899124/150‑1 (MQ=255)

                                                                         cTGCGTTGTTGCAGGGCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACa                                                                              >  1:8812016/1‑150 (MQ=255)

                                                                          tGCGTTGTTGCAGGGCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACaa                                                                             >  1:4997272/1‑150 (MQ=255)

                                                                          tGCGTTGTTGCAGGGCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACaa                                                                             >  2:8303086/1‑150 (MQ=255)

                                                                            cGTTGTTGCAGGGCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGTTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAgg                                                                           <  2:6829402/150‑1 (MQ=255)

                                                                            cGTTGTTGCAGGGCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAgg                                                                           <  1:7339002/150‑1 (MQ=255)

                                                                             gttgttGCAGGGCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGa                                                                          <  1:589054/150‑1 (MQ=255)

                                                                               tgttgCAGGGCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGt                                                                        <  1:7924800/150‑1 (MQ=255)

                                                                                 ttgCAGGGCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCgg                                                                                                                    <  1:7643366/104‑1 (MQ=255)

                                                                                      gggCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGa                                                                 <  2:5467712/150‑1 (MQ=255)

                                                                                       ggCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAg                                                                <  2:2768982/150‑1 (MQ=255)

                                                                                       ggCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGGTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGTGGCGACCTACCAGGAGTATCAGGAg                                                                >  2:7299585/1‑150 (MQ=255)

                                                                                          aGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCaa                                                             <  1:2566391/150‑1 (MQ=255)

                                                                                          aGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCaa                                                             <  1:2262336/150‑1 (MQ=255)

                                                                                           gCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAg                                                            <  1:1414618/150‑1 (MQ=255)

                                                                                             cTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGc                                                          >  2:8714857/1‑150 (MQ=255)

                                                                                              tGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGcc                                                         >  2:4043653/1‑150 (MQ=255)

                                                                                                ccGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTg                                                       >  1:5132305/1‑150 (MQ=255)

                                                                                                ccGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTg                                                       <  2:2455131/150‑1 (MQ=255)

                                                                                                 cGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTgg                                                      >  2:4449656/1‑150 (MQ=255)

                                                                                                   gAGTTGCGCGGTTAGGGGCTGGTCGACCTCTGCGTGCGCTTTGCCGGCTTGGCGCTCGTGGTGGTGTTCTCTGCCTGTTTTTGCGGTTTTGCCCCGCAGTTCGTGGGGCTTGAGGCGACCTACTAGGTGTATCACGTGCAAGGCCTGGaa                                                    <  2:6517805/150‑1 (MQ=255)

                                                                                                      ttGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTGGAAATg                                                 <  1:6259255/150‑1 (MQ=255)

                                                                                                        gcgcGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTGGAAATGCt                                               <  1:695895/150‑1 (MQ=255)

                                                                                                         cgcgGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTGGAAATGCTc                                              <  1:1834475/150‑1 (MQ=255)

                                                                                                            ggTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTGGAAATGCTCGGc                                           <  2:5664272/150‑1 (MQ=255)

                                                                                                              tAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTGGAAATGCTCGGCGt                                         <  2:6089019/150‑1 (MQ=255)

                                                                                                                aGGGGCAGGTCGACCTCGGCGAGCCCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTGGAAATGCTCGGCGTGc                                       <  2:8538856/150‑1 (MQ=255)

                                                                                                                 ggggCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTGGAAATGCTCGGCGTGcc                                      >  1:1535192/1‑150 (MQ=255)

                                                                                                                  tggCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTGGAAATGCTCGGCGTGCCg                                     >  1:3480720/2‑150 (MQ=255)

                                                                                                                  gggCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTGTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTGGAAATGCTCGGCGTGCCg                                     <  2:8654705/150‑1 (MQ=255)

                                                                                                                    gCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTGGAAATGCTCGGCGTGCCGTc                                   >  2:896985/1‑150 (MQ=255)

                                                                                                                    gCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTGGAAATGCTCGGCGTGCCGTc                                   >  2:115319/1‑150 (MQ=255)

                                                                                                                    gCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTGGAAATGCTCGGCGTGCCGTc                                   <  1:8303384/150‑1 (MQ=255)

                                                                                                                        gTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTGGAAATGCTCGGCGTGCCGTCCAAt                               >  1:1150644/1‑150 (MQ=255)

                                                                                                                        gTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTGGAAATGCTCGGCGTGCCGTCCAAt                               <  1:5367141/150‑1 (MQ=255)

                                                                                                                        gTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTGGAAATGCTCGGCGTGCCGTCCAAt                               <  2:2095087/150‑1 (MQ=255)

                                                                                                                        gTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTGGAAATGCTCGGCGTGCCGTCCAAt                               <  2:3531262/150‑1 (MQ=255)

                                                                                                                           gACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTGGAAATGCTCGGCGTGCCGTCCAATGAc                            <  2:2046793/150‑1 (MQ=255)

                                                                                                                             ccaccacGAGCGGCTTGCCGGCAAAGCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTGGAAATGCTCGGCGTGCCGTc                                   >  1:5749829/7‑141 (MQ=255)

                                                                                                                              ctctGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTGGAAATGCTCGGCGTGCCGTCCAATGACTTc                         <  1:235563/150‑1 (MQ=255)

                                                                                                                               accacGAGCGGCTTGCCGGCAAAGCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTGGAAATGCTCGGCGTGCCGTCCAATGACTTCa                        >  2:3410903/5‑150 (MQ=255)

                                                                                                                                ctGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTGGAAATGCTCGGCGTGCCGTCCAATGACTTCaa                       >  2:3811152/1‑150 (MQ=255)

                                                                                                                                ctGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTGGAAATGCTCGGCGTGCCGTCCAATGACTTCaa                       >  2:3275711/1‑150 (MQ=255)

                                                                                                                                 tGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTGGAAATGCTCGGCGTGCCGTCCAATGACTTCAAg                      >  1:8224830/1‑150 (MQ=255)

                                                                                                                                 cacGAGCGGCTTGCCGGCAAAGCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTGGAAATGCTCGGCGTGCCGTCCAATGACTTCAAg                      >  2:6434716/3‑150 (MQ=255)

                                                                                                                                 cacGAGCGGCTTGCCGGCAAAGCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTGGAAATGCTCGGCGTGCCGTCCAATGACTTCAAg                      >  1:4399582/3‑150 (MQ=255)

                                                                                                                                  gCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTGGAAATGCTCGGCGTGCCGTCCAATGACTTCAAGc                     <  2:5707707/150‑1 (MQ=255)

                                                                                                                                  gCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTGGAAATGCTCGGCGTGCCGTCCAATGACTTCAAGc                     <  2:4646242/150‑1 (MQ=255)

                                                                                                                                  gCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTGGAAATGCTCGGCGTGCCGTCCAATGACTTCAAGc                     <  2:7462270/150‑1 (MQ=255)

                                                                                                                                    gAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTGGAAATGCTCGGCGTGCCGTCCAATGACTTCAAGCAg                   <  1:4207957/150‑1 (MQ=255)

                                                                                                                                     aGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTGGAAATGCTCGGCGTGCCGTCCAATGACTTCAAGCAgg                  <  2:372021/150‑1 (MQ=255)

                                                                                                                                       cgcTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTGGAAATGCTCGGCGTGCCGTCCAATGACTTCAAGCAGGaa                <  1:1755436/150‑1 (MQ=255)

                                                                                                                                       cgcTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTGGAAATGCTCGGCGTGCCGTCCAATGACTTCAAGCAGGaa                <  1:8002984/150‑1 (MQ=255)

                                                                                                                                       cgcTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGATTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTGGAAATGCTCGGCGTGCCGTCCAATGACTTCAAGCAGGaa                <  1:4607890/150‑1 (MQ=255)

                                                                                                                                         cTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTGGAAATGCTCGGCGTGCCGTCCaa                                >  1:5582083/1‑132 (MQ=255)

                                                                                                                                          tttGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTGGAAATGCTCGGCGTGCCGTCCAATGACTTCAAGCAGGAAGAc             >  2:8201922/1‑150 (MQ=255)

                                                                                                                                          tttGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTGGAAATGCTCGGCGTGCCGTCCAATGACTTCAAGCAGGAAGAc             >  1:148750/1‑150 (MQ=255)

                                                                                                                                            tGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTGGAAATGCTCGGCGTGCCGTCCAATGACTTCAAGCAGGAAGACGc           <  2:4164831/150‑1 (MQ=255)

                                                                                                                                              ccGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCAGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTGGAAATGCTCGGCGTGCCGTCCAATGACTTCAAGCAGGAAGACGCCg         >  2:2686063/1‑150 (MQ=255)

                                                                                                                                               cGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTGGAAATGCTCGGCGTGCCGTCCAATGACTTCAAGCAGGAAGACGCCGa        <  1:3343699/150‑1 (MQ=255)

                                                                                                                                                 gCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTGGAAATGCTCGGCGTGCCGTCCAATGACTTCAAGCAGGAAGACGCCGACa      >  1:3832406/1‑150 (MQ=255)

                                                                                                                                                 gCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTGGAAATGCTCGGCGTGCCGTCCAATGACTTCAAGCAGGAAGACGCCGACa      <  1:2654355/150‑1 (MQ=255)

                                                                                                                                                   aaGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTGGAAATGCTCGGCGTGCCGTCCAATGACTTCAAGCAGGAAGACGCCGACAGt    <  2:7502277/150‑1 (MQ=255)

                                                                                                                                                    aGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTGGAAATGCTCGGCGTGCCGTCCAATGACTTCAAGCAGGAAGACGCCGACAGTg   <  1:649580/150‑1 (MQ=255)

                                                                                                                                                     gCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTGGAAATGCTCGGCGTGCCGTCCAATGACTTCAAGCAGGAAGACGCCGACAGTga  <  2:4526073/150‑1 (MQ=255)

                                                                                                                                                     gCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTGGAAATGCTCGGCGTGCCGTCCAATGACTTCAAGCAGGAAGACGCCGACAGTga  >  1:3875542/1‑150 (MQ=255)

ATGCGTGCACATTGGTTGACGGTGCCGGTGCTGGCCCTGCTTGCCAGTACATCGAGCTGGGCTGCAGATTGCCCTGCGTTGTTGCAGGGCAGCCTGCCGGAGTTGCGCGGTAAGGGGCAGGTCGACCTCTGCGAGCGCTTTGCCGGCAAGCCGCTCGTGGTGGTGAACACTGCCAGTTATTGCGGTTTTGCCCCGCAGTTCGAGGGGCTTGAGGCGACCTACAAGGAGTATCACGAGCAAGGCCTGGAAATGCTCGGCGTGCCGTCCAATGACTTCAAGCAGGAAGACGCCGACAGTGA  >  NC_002947/1884256‑1884554

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