BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,826,370	G→C	5.5%	G48A (GGC→GCC)	PP_5102 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,826,370	0	G	C	5.5%	1174.6 / 40.4	347	G48A (GGC→GCC)	PP_5102	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (181/146); new base C (6/13); total (187/160)
Fisher's exact test for biased strand distribution p-value = 5.72e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 363 total aligned reads displayed.

TGACCCATGGATCTGGTACTTGATCTGCTCGCGACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTGTGTTACTGGTGGTGCTGG  >  NC_002947/5826222‑5826519

  aCCCATGGATCTGGTACTTGATCTGCTCGCGACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCt                                                                                                                                                    >  1:5683629/1‑150 (MQ=255)

    ccATGGATCTGGTACTTGATCTGCTCGCGACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTcg                                                                                                                                                  >  1:8223417/1‑150 (MQ=255)

    ccATGGATCTGGTACTTGATCTGCTCGCGACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTcg                                                                                                                                                  >  1:2653430/1‑150 (MQ=255)

     cATGGATCTGGTACTTGATCTGCTCGCGACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTcgc                                                                                                                                                 >  2:3072881/1‑150 (MQ=255)

     cATGGATCTGGTACTTGATCTGCTCGCGACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTcgc                                                                                                                                                 <  1:854567/150‑1 (MQ=255)

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

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

      aTGGATCTGGTACTTGATCTGCTCGCGACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACAGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTcgcg                                                                                                                                                <  2:4005628/150‑1 (MQ=255)

             tGGTACTTGATCTGCTCGCGACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGccc                                                                                                                                         >  1:5039060/1‑150 (MQ=255)

             tGGTACTTGATCTGCTCGCGACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGccc                                                                                                                                         >  2:6613315/1‑150 (MQ=255)

              ggTACTTGATCTGCTCGCGACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGCTGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGCCCGCGAGGCCGccca                                                                                                                                        <  2:451326/150‑2 (MQ=255)

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

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

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

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

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

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

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

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

                             cgcgACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCAt                                                                                                                         <  1:6071335/150‑1 (MQ=255)

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

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

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

                                   ggTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCt                                                                                                                   >  1:7333777/1‑150 (MQ=255)

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

                                      ttCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCAGGGCGCCCGCGAGGCCGcccagg                                                                                                                                      >  1:4833011/1‑125 (MQ=255)

                                       tCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTTAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATg                                                                                                               <  2:8029182/150‑1 (MQ=255)

                                       tCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATg                                                                                                               <  1:4001002/150‑1 (MQ=255)

                                       tCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATg                                                                                                               >  1:5417657/1‑150 (MQ=255)

                                         ccGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGt                                                                                                             >  1:789811/1‑150 (MQ=255)

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

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

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

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

                                                   cGCAGCAACCTGTCGGAGATTTCACTGGCCGTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGt                                                                                                   >  2:8805960/1‑150 (MQ=255)

                                                    gcagcaACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTa                                                                                                  <  1:4367105/150‑1 (MQ=255)

                                                       gcaACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCg                                                                                               >  1:1584424/1‑150 (MQ=255)

                                                            cTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGtgct                                                                                          >  2:5672043/1‑150 (MQ=255)

                                                            cTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGGGCCGGTCATGGCGCTGATGGTCATGATCGGTAGCGGtgct                                                                                          >  2:7321714/1‑150 (MQ=255)

                                                                cGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTgcgc                                                                                      <  2:2450922/150‑1 (MQ=255)

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

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

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

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

                                                                        tCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTAt                                                                              <  1:1742991/150‑1 (MQ=255)

                                                                        tCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTAt                                                                              <  1:2780919/150‑1 (MQ=255)

                                                                        tCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTAt                                                                              <  1:4534186/150‑1 (MQ=255)

                                                                        tCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTAt                                                                              <  1:4733353/150‑1 (MQ=255)

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

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

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

                                                                          aCTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGc                                                                            >  1:7776071/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:7111773/150‑1 (MQ=255)

                                                                                    gTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCgtgg                                                                  <  2:548786/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)

                                                                                      aGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCgtgggt                                                                <  2:5125765/150‑1 (MQ=255)

                                                                                               gctgGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAggg                                                       >  1:7084278/1‑150 (MQ=255)

                                                                                                 tgGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGACGGGGc                                                     >  2:8656106/1‑150 (MQ=255)

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

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

                                                                                                   gTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGCTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTc                                                   >  1:8207471/1‑150 (MQ=255)

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

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

                                                                                                           cGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGtt                                           <  1:1119391/150‑1 (MQ=255)

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

                                                                                                               aCCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTcaac                                       >  1:8138963/1‑150 (MQ=255)

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

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

                                                                                                                      tCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATg                                <  1:3079830/150‑1 (MQ=255)

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

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

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

                                                                                                                        aaaGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGCAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATgcg                              <  2:7649542/150‑1 (MQ=255)

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

                                                                                                                               gggCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGgcgc                       <  1:4508596/150‑1 (MQ=255)

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

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

                                                                                                                                 gCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCt                     <  1:5535829/150‑1 (MQ=255)

                                                                                                                                  ccGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCtg                    >  2:1067922/1‑150 (MQ=255)

                                                                                                                                   cGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCtgt                   <  2:7587380/150‑1 (MQ=255)

                                                                                                                                    gAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCtgtg                  <  1:7946941/150‑1 (MQ=255)

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

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

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

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

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

                                                                                                                                           gcCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTGTGTTACtgg           <  1:6821679/150‑1 (MQ=255)

                                                                                                                                           gcCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTGTGTTACtgg           >  1:7539518/1‑150 (MQ=255)

                                                                                                                                            gcagggCGCCCGCGAGGCCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTGTGTTACtggg          >  2:1879550/4‑149 (MQ=255)

                                                                                                                                            ccTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGTGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTGTGTTACtggt          <  2:1410730/150‑1 (MQ=255)

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

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

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

                                                                                                                                                   ggcCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTGTGTTACTGGTGGTGCTg   >  1:2948058/1‑150 (MQ=255)

                                                                                                                                                    gcCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTGTGTTACTGGTGGTGCTgg  <  2:465592/150‑1 (MQ=255)

TGACCCATGGATCTGGTACTTGATCTGCTCGCGACGGTTTCCCGCTGGAGCCGCAGCAACCTGTCGGAGATTTCACTGGCCTTGGTAGGCTGCCTGCTGGTGCTGTTCGGCACCGATGTCAAAGCCTGGGCCGAGCAGCGCCTGGGCGGCCTCGCGGGCGCCCTGCGCGTGCCGTTCATGGCGCTGATGGTCATGATCGGTAGCGGTGCTGCGCTGATCTATGCCACACCGTGGGTGGTGAAGGGGCTCAGCCAGTTCAACAACTATGCGCTGGCGCCTGTGTTACTGGTGGTGCTGG  >  NC_002947/5826222‑5826519

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