BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,826,376	C→A	5.4%	A50E (GCG→GAG)	PP_5102 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,826,376	0	C	A	5.4%	1166.1 / 33.8	345	A50E (GCG→GAG)	PP_5102	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (167/158); new base A (6/13); total (173/172)
Fisher's exact test for biased strand distribution p-value = 1.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 363 total aligned reads displayed.

CATGGATCTGGTACTTGATCTGCTCGCGACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTGTGTTACTGGTGGTGCTGGTATTGA  >  NC_002947/5826227‑5826525

 aTGGATCTGGTACTTGATCTGCTCGCGACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTcgcg                                                                                                                                                      >  1:2353785/1‑150 (MQ=255)

 aTGGATCTGGTACTTGATCTGCTCGCGACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTcgcg                                                                                                                                                      >  1:6182919/1‑150 (MQ=255)

         ggTACTTGATCTGCTCGCGACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGCCCGCGAGGCCGccca                                                                                                                                              <  2:240555/150‑2 (MQ=255)

         ggTACTTGATCTGCTCGCGACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGCCCGCGAGGCCGccca                                                                                                                                              <  2:2330158/150‑2 (MQ=255)

           gaCTTGATCTGCTCGCGACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTgc                                                                                                                                            <  2:661008/149‑1 (MQ=255)

            aCTTGATCTGCTCGCGACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTgcg                                                                                                                                           <  1:8666617/150‑1 (MQ=255)

             cTTGATCTGCTCGCGACGGTTTCCCGCTGGAGCCGCAGCAACCTGTTGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTgcgc                                                                                                                                          >  1:5145108/1‑150 (MQ=255)

             cTTGATCTGCTCGCGACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTgcgc                                                                                                                                          >  2:6775649/1‑150 (MQ=255)

              ttGATCTGCTCGCGACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGcgcg                                                                                                                                         >  2:1314617/1‑150 (MQ=255)

               tGATCTGCTCGCGACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGt                                                                                                                                        <  2:4956246/150‑1 (MQ=255)

                  tCTGCTCGCGACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGcc                                                                                                                                     <  1:7930275/150‑1 (MQ=255)

                    tGCTCGCGACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGt                                                                                                                                   >  1:1192606/1‑150 (MQ=255)

                    tGCTCGCGACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGt                                                                                                                                   >  1:465592/1‑150 (MQ=255)

                     gCTCGCGACGGTTACCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGtt                                                                                                                                  <  1:4510238/150‑1 (MQ=255)

                       tCGCGACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTct                                                                                                                                >  1:5492191/1‑149 (MQ=255)

                       tCGCGACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCa                                                                                                                                >  2:2878624/1‑150 (MQ=255)

                         gcgACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATg                                                                                                                              >  1:3057846/1‑150 (MQ=255)

                         gcgACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATg                                                                                                                              >  2:3385739/1‑150 (MQ=255)

                         gcgACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATg                                                                                                                              >  1:17355/1‑150 (MQ=255)

                         gcgACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATg                                                                                                                              <  1:4879871/150‑1 (MQ=255)

                         gcgACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATg                                                                                                                              >  2:8546763/1‑150 (MQ=255)

                           gACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGCCCGCGAGGCCGcccagg                                                                                                                                            <  2:2077481/134‑4 (MQ=255)

                           gACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGCCCGCGAGGCCGcccagg                                                                                                                                            >  1:2077669/1‑131 (MQ=255)

                                 ttCCCGCTGGGGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGAt                                                                                                                      >  2:5535829/1‑150 (MQ=255)

                                     cGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTc                                                                                                                  <  1:5692012/150‑1 (MQ=255)

                                      gCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCa                                                                                                                 <  2:8721905/150‑1 (MQ=255)

                                       cTGGAGCCGCAGCTACGTGTCGGAGATTTCACTGGCCTTGGTTGGCTGCCTGGTGGTTCTTTTCGGCAGCGATGTCAAAGCCTGGTCCTAGCAGCGCCTGGGCGGCCTCGGGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCAt                                                                                                                <  1:7713743/150‑1 (MQ=255)

                                       cTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCAt                                                                                                                <  1:594748/150‑1 (MQ=255)

                                         ggAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGa                                                                                                              <  1:8266158/150‑1 (MQ=255)

                                             ccGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCgg                                                                                                          <  1:8355802/150‑1 (MQ=255)

                                                          tCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTgcg                                                                                             >  1:7571405/1‑150 (MQ=255)

                                                            ggAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGTCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCt                                                                                           <  1:9151772/150‑1 (MQ=255)

                                                            ggAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCt                                                                                           <  1:2601346/150‑1 (MQ=255)

                                                            ggAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCt                                                                                           >  2:7738111/1‑150 (MQ=255)

                                                               gaTTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGAt                                                                                        <  1:1530569/150‑1 (MQ=255)

                                                                aTTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATc                                                                                       >  2:6029095/1‑150 (MQ=255)

                                                                   tCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAACGCCTGGGCCTAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTAt                                                                                    <  2:6981316/150‑1 (MQ=255)

                                                                   tCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTTCGCTGATCTAt                                                                                    <  1:285282/150‑1 (MQ=255)

                                                                     aCTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGc                                                                                  >  2:2111266/1‑150 (MQ=255)

                                                                     aCTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGc                                                                                  >  2:462422/1‑150 (MQ=255)

                                                                       tGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCca                                                                                >  2:5092117/1‑150 (MQ=255)

                                                                           cTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACAcc                                                                            >  2:1721668/1‑150 (MQ=255)

                                                                             tGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCgt                                                                          <  1:3718366/150‑1 (MQ=255)

                                                                              ggTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCgtg                                                                         <  2:7902873/150‑1 (MQ=255)

                                                                               gTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCgtgg                                                                        <  1:6891182/150‑1 (MQ=255)

                                                                               gTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCgtgg                                                                        <  2:548822/150‑1 (MQ=255)

                                                                               gTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCgtgg                                                                        <  2:6466771/150‑1 (MQ=255)

                                                                               gTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCgtgg                                                                        <  2:273109/150‑1 (MQ=255)

                                                                                tAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCgtggg                                                                       >  1:6962453/1‑150 (MQ=255)

                                                                                   gCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCgtgggtgg                                                                    <  2:3973275/150‑1 (MQ=255)

                                                                                    ctgcctgcTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGggtggt                                                                   <  1:4812244/150‑1 (MQ=255)

                                                                                     tgcctgcTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTAt                                                                                    >  2:4500239/1‑132 (MQ=255)

                                                                                        ctgctgGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAg                                                               <  1:4326519/150‑1 (MQ=255)

                                                                                         tgctgGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAgg                                                              <  2:5710210/150‑1 (MQ=255)

                                                                                         tgctgGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAgg                                                              <  1:6079162/150‑1 (MQ=255)

                                                                                         tgctgGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAgg                                                              <  1:4404994/150‑1 (MQ=255)

                                                                                          gctgGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGTGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCAGGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGGAggg                                                             >  1:7084343/1‑150 (MQ=255)

                                                                                             ggTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCt                                                          <  2:1386682/150‑1 (MQ=255)

                                                                                              gTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTc                                                         >  1:8210704/1‑150 (MQ=255)

                                                                                              gTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTc                                                         >  1:8511532/1‑150 (MQ=255)

                                                                                              gTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTc                                                         >  2:2666107/1‑150 (MQ=255)

                                                                                              gTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTc                                                         >  2:4043326/1‑150 (MQ=255)

                                                                                               tGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCTCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGCCTCa                                                        <  2:4908788/150‑1 (MQ=255)

                                                                                               tGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCa                                                        <  1:7172438/150‑1 (MQ=255)

                                                                                                   gTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCa                                                    <  2:3048893/150‑1 (MQ=255)

                                                                                                         cACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTcaa                                              <  1:8206782/150‑1 (MQ=255)

                                                                                                          aCCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTcaac                                             <  2:5182284/150‑1 (MQ=255)

                                                                                                             gATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCaacaac                                          >  2:2126822/1‑150 (MQ=255)

                                                                                                             gATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCaacaac                                          <  1:8097891/150‑1 (MQ=255)

                                                                                                                gTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTAt                                       >  2:439084/1‑150 (MQ=255)

                                                                                                                 tCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATg                                      >  1:5208415/1‑150 (MQ=255)

                                                                                                                   aaaGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATgcg                                    >  1:146615/1‑150 (MQ=255)

                                                                                                                   aaaGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATgcg                                    >  1:9061905/1‑150 (MQ=255)

                                                                                                                   aaaGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATgcg                                    >  2:6181183/1‑150 (MQ=255)

                                                                                                                     aGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCt                                  <  2:7487699/150‑1 (MQ=255)

                                                                                                                      gCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTg                                 >  1:3858864/1‑150 (MQ=255)

                                                                                                                        cTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGgc                               >  1:6330440/1‑150 (MQ=255)

                                                                                                                         tGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGgcg                              >  1:4559972/1‑150 (MQ=255)

                                                                                                                           ggCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGcc                            >  1:3309775/1‑150 (MQ=255)

                                                                                                                            gCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCt                           >  2:2396544/1‑150 (MQ=255)

                                                                                                                                agcagcGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTgtgt                       <  2:90161/150‑1 (MQ=255)

                                                                                                                                agcagcGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTgtgt                       <  1:367805/150‑1 (MQ=255)

                                                                                                                                   agcGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTGTGTTAc                    >  2:7086610/1‑150 (MQ=255)

                                                                                                                                     cgcCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTGTGTTACtg                  >  1:2001005/1‑150 (MQ=255)

                                                                                                                                     cgcCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTGTGTCACtg                  <  2:3787906/150‑1 (MQ=255)

                                                                                                                                        cagggCGCCCGCGAGGCCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTGTGTTACtggtg               >  2:2568528/3‑150 (MQ=255)

                                                                                                                                        cTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTGTGTTACtggtg               <  1:5299061/150‑1 (MQ=255)

                                                                                                                                         tGGGCGGCCTCACGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTGTGTTACtggtgg              >  1:78562/1‑150 (MQ=255)

                                                                                                                                          gggCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTGTGTTACTggtggt             <  2:1192606/150‑1 (MQ=255)

                                                                                                                                             cggcCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTGTGTTACTGGTGGTGCt          <  2:760670/150‑1 (MQ=255)

                                                                                                                                              ggcCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTGTGTTACTGGTGGTGc           <  1:7469489/148‑1 (MQ=255)

                                                                                                                                              ggcCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTGTGTTACTGGTGGTGCTg         <  2:4666338/150‑1 (MQ=255)

                                                                                                                                              ggcCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTGTGTTACTGGTGGTGCTg         <  2:8232167/150‑1 (MQ=255)

                                                                                                                                              ggcCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTGTGTTACTGGTGGTGCTg         <  2:1881888/150‑1 (MQ=255)

                                                                                                                                                 cTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTGTGTTACTGGTGGTGCTGGTa      >  2:7399509/1‑150 (MQ=255)

                                                                                                                                                 cTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTGTGTTACTGGTGGTGCTGGTa      <  1:3101820/150‑1 (MQ=255)

                                                                                                                                                 cTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTGTGTTACTGGTGGTGCTGGTa      >  2:5331176/1‑150 (MQ=255)

                                                                                                                                                  tCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTGTGTTACTGGTGGTGCTGGTat     >  2:7653632/1‑150 (MQ=255)

                                                                                                                                                  tCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTGTGTTACTGGTGGTGCTGGTat     >  2:7986799/1‑150 (MQ=255)

                                                                                                                                                    gcgGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTGTGTTACTGGTGGTGCTGGTattg   <  1:1037558/150‑1 (MQ=255)

CATGGATCTGGTACTTGATCTGCTCGCGACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTGTGTTACTGGTGGTGCTGGTATTGA  >  NC_002947/5826227‑5826525

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