BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,967,040	C→G	5.3%	T218S (ACC→AGC) ‡	PP_1762 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,967,040	0	C	G	5.3%	506.7 / 13.5	151	T218S (ACC→AGC) ‡	PP_1762	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (58/85); new base G (3/5); total (61/90)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

Only 100 of 154 total aligned reads displayed.

TGTGCTCGACCGACACCTATGCGCAGGGGTTCCCTGAAGAAGGGCGCAAGCGGGCGATGAAGTGGCTGAAGATGCATGCCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCG  >  NC_002947/1966893‑1967189

tgtgCTCGACCGACACCTATGCGCAGGGGTTCCCTGAAGAAGGGCGCAAGCGGGCGATGAAGTGGCTGAAGATGCATGCCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCAccc                                                                                                                                                     >  2:4312632/1‑150 (MQ=255)

tgtgCTCGACCGACACCTATGCGCAGGGGTTCCCTGAAGAAGGGCGCAAGCGGGCGATGAAGTGGCTGAAGATGCATGCCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCAccc                                                                                                                                                     <  2:4016066/150‑1 (MQ=255)

  tgCTCGACCGACACCTATGCGCAGGGGTTCCCTGAAGAAGGGCGCAAGCGGGCGATGAAGTGGCTGAAGATGCATGCCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTg                                                                                                                                                   <  1:4205334/150‑1 (MQ=255)

    cTCGACCGACACCTATGCGCAGGGGTTCCCTGAAGAAGGGCGCAAGCGGGCGATGAAGTGGCTGAAGATGCATGCCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTggg                                                                                                                                                 >  2:2569458/1‑150 (MQ=255)

      cgaccgacACCTATGCGCAGGGGTTCCCTGAAGAAGGGCGCAAGCGGGCGATGAAGTGGCTGAAGATGCATGCCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGt                                                                                                                                               <  2:4843129/150‑1 (MQ=255)

         ccgacAGCTATGCGCAGTGGTTCACTGAAGAAGGGCGCAAGCGGGCGATGAAGTGGCTGAAGATGCTTGGCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGTAAATCATCTGTACCCTGTCCGGCGAGAACCCCACCCTGGGGTTGc                                                                                                                                            <  1:3882627/150‑1 (MQ=255)

          cgacACCTATGCGCAGGGGTTCCCTGAAGAAGGGCGCAAGCGGGCGATGAAGTGGCTGAAGATGCATGCCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCAGGGTGGGGttctc                                                                                                                                           <  2:3532137/150‑4 (MQ=255)

          cgacACCTATGCGCAGGGGTTCCCTGAAGAAGGGCGCAAGCGGGCGATGAAGTGGCTGAAGATGCATGCCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCg                                                                                                                                           >  2:400469/1‑150 (MQ=255)

            acacCTATGCGCAGGGGTTCCCTGAAGAAGGGCGCAAGCGGGCGATGAAGTGGCTGAAGATGCATGCCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATGTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGa                                                                                                                                         >  2:3481284/1‑150 (MQ=255)

            acacCTATGCGCAGGGGTTCCCTGAAGAAGGGCGCAAGCGGGCGATGAAGTGGCTGAAGATGCATGCCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGa                                                                                                                                         <  1:1334507/150‑1 (MQ=255)

             cacCTATGCGCAGGGGGTCCCTGAAGAAGGGCGCAAGCGGGCGATGAAGTGGCTGAAGATGCATGCCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGAc                                                                                                                                        <  2:2307113/150‑1 (MQ=255)

              acCTATGCGCAGGGGTTCCCTGAAGAAGGGCGCAAGCGGGCGATGAAGTGGCTGAAGATGCATGCCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGAcc                                                                                                                                       <  1:1655403/150‑1 (MQ=255)

                cTATGCGCAGGGGTTCCCTGAAGAAGGGCGCAAGCGGGCGATGAAGTGGCTGAAGATGCATGCCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGa                                                                                                                                     <  1:4581286/150‑1 (MQ=255)

                cTATGCGCAGGGGTTCCCTGAAGAAGGGCGCAAGCGGGCGATGAAGTGGCTGAAGATGCATGCCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGa                                                                                                                                     >  2:3273754/1‑150 (MQ=255)

                   tGCGCAGGGGTTCCCTGAAGAAGGGCGCAAGCGGGCGATGAAGTGGCTGAAGATGCATGCCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGAcc                                                                                                                                       <  1:1361015/145‑1 (MQ=255)

                    gcgcAGGGGTTCCCTGAAGAAGGGCGCAAGCGGGCGATGAAGTGGCTGAAGATGCATGCCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTg                                                                                                                                 <  1:2749667/150‑1 (MQ=255)

                     cgcAGGGGTTCCCTGAAGAAGGGCGCAAGCGGGCGATGAAGTGGCTGAAGATGCATGCCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTGCCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTgc                                                                                                                                <  2:2103953/150‑1 (MQ=255)

                     cgcAGGGGTTCCCTGAAGAAGGGCGCAAGCGGGCGATGAAGTGGCTGAAGATGCATGCCCAGGATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTGCCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTgc                                                                                                                                <  2:2103797/150‑1 (MQ=255)

                        aGGGGTTCCCTGAAGAAGGGCGCAAGCGGGCGATGAAGTGGCTGAAGATGCATGCCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCa                                                                                                                             >  1:2370566/1‑150 (MQ=255)

                          gggTTCCCTGAAGAAGGGCGCAAGCGGGCGATGAAGTGGCTGAAGATGCATGCCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGttctc                                                                                                                                           <  1:3320810/134‑4 (MQ=255)

                          gggTTCCCTGAAGAAGGGCGCAAGCGGGCGATGAAGTGGCTGAAGATGCATGCCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAgg                                                                                                                           >  2:501785/1‑150 (MQ=255)

                          gggTTCCCTGAAGAAGGGCGCAAGCGGGCGATGAAGTGGCTGAAGATGCATGCCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAgg                                                                                                                           >  1:2949155/1‑150 (MQ=255)

                                ccTGAAGAAGGGCGCAAGCGGGCGATGAAGTGGCTGAAGATGCATGCCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATc                                                                                                                     <  1:2409226/150‑1 (MQ=255)

                                 cTGAAGAAGGGCGCAAGCGGGCGATGAAGTGGCTGAAGATGCATGCCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCt                                                                                                                    >  1:914177/1‑150 (MQ=255)

                                  tGAAGAAGGGCGCAAGCGGGCGATGAAGTGGCTGAAGATGCATGCCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTg                                                                                                                   >  1:391114/1‑150 (MQ=255)

                                  tGAAGAAGGGCGCAAGCGGGCGATGAAGTGGCTGAAGATGCATGCCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTg                                                                                                                   >  1:2948827/1‑150 (MQ=255)

                                   gaagaaGGGCGCAAGCGGGCGATGAAGTGGCTGAAGATGCATGCCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGc                                                                                                                  <  1:4681428/150‑1 (MQ=255)

                                   gaagaaGGGCGCAAGCGGGCGATGAAGTGGCTGAAGATGCATGCCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGc                                                                                                                  <  2:249723/150‑1 (MQ=255)

                                   gaagaaGGGCGCAAGCGGGCGATGAAGTGGCTGAAGATGCATGCCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGc                                                                                                                  <  1:602605/150‑1 (MQ=255)

                                     agaagGGCGCAAGCGGGCGATGAAGTGGCTGAAGATGCATGCCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCaa                                                                                                                <  2:2543664/150‑1 (MQ=255)

                                         gggCGCAAGCGGGCGATGAAGTGGCTGAAGATGCATGCCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGt                                                                                                            >  1:412406/1‑150 (MQ=255)

                                          ggCGCAAGCGGGCGATGAAGTGGCTGAAGATGCATGCCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGtt                                                                                                           <  1:1825144/150‑1 (MQ=255)

                                          ggCGCAAGCGGGCGATGAAGTGGCTGAAGATGCATGCCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGtt                                                                                                           <  1:5069471/150‑1 (MQ=255)

                                          ggCGCAAGCGGGCGATGAAGTGGCTGAAGATGCATGCCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGtt                                                                                                           <  1:266953/150‑1 (MQ=255)

                                           gcgcAAGCGGGCGATGAAGTGGCTGAAGATGCATGCCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTa                                                                                                          >  2:3559893/1‑150 (MQ=255)

                                           gcgcAAGCGGGCGATGAAGTGGCTGAAGATGCATGCCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTa                                                                                                          >  2:612301/1‑150 (MQ=255)

                                            cgcAAGCGGGCGATGAAGTGGCTGAAGATGCATGCCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTAc                                                                                                         >  2:3866885/1‑150 (MQ=255)

                                               aaGCGGGCGATGAAGTGGCTGAAGATGCATGCCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGAc                                                                                                      <  1:4775487/150‑1 (MQ=255)

                                               aaGCGGGCGATGAAGTGGCTGAAGATGCATGCCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGAc                                                                                                      >  1:1735775/1‑150 (MQ=255)

                                                  cGGGCGATGAAGTGGCTGAAGATGCATGCCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGc                                                                                                   <  2:2487043/150‑1 (MQ=255)

                                                        aTGAAGTGGCTGAAGATGCATGCCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCAtgttc                                                                                             <  1:3122068/150‑1 (MQ=255)

                                                               ggCTGAAGATGCATGCCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCAtgttcctgttcc                                                                                      <  2:1575537/150‑1 (MQ=255)

                                                                gCTGAAGATGCATGCCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATgttcctgttcct                                                                                     <  2:408927/150‑1 (MQ=255)

                                                                gCTGAAGATGCATGCCATGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGTTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGAATgttcctgttcct                                                                                     <  2:3404325/150‑1 (MQ=255)

                                                                  tGAAGATGCATGCCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTgg                                                                                   >  2:175508/1‑150 (MQ=255)

                                                                  tGAAGATGCATGCCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTgg                                                                                   >  2:1722670/1‑150 (MQ=255)

                                                                  tGAAGATGCATGCCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTgg                                                                                   >  1:1567509/1‑150 (MQ=255)

                                                                    aaGATGCATGCCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAg                                                                                 >  1:4754904/1‑150 (MQ=255)

                                                                        tgcatgcCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCt                                                                             <  1:2308484/150‑1 (MQ=255)

                                                                        tgcatgcCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCt                                                                             <  1:2153846/150‑1 (MQ=255)

                                                                           atgcCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCtgca                                                                          <  1:614040/150‑1 (MQ=255)

                                                                           atgcCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCtgca                                                                          >  1:3145951/1‑150 (MQ=255)

                                                                            tgcCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCtgcat                                                                         >  1:494488/1‑150 (MQ=255)

                                                                             gcCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCtgcatg                                                                        >  1:2124840/1‑150 (MQ=255)

                                                                             gcCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCtgcatg                                                                        >  2:4169069/1‑150 (MQ=255)

                                                                              cccAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCtgcatgc                                                                       <  2:3645344/150‑1 (MQ=255)

                                                                              cccAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCtgcatgc                                                                       <  2:3064982/150‑1 (MQ=255)

                                                                               ccAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCtgcatgca                                                                      >  1:757647/1‑150 (MQ=255)

                                                                               ccAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCtgcatgca                                                                      >  1:4353121/1‑150 (MQ=255)

                                                                                  ttATGACGATGCGCACCCGAGGGATGCGCTTGAAATCATCTGTACCCTGGCCGGCGAGAACCACACCCTGGGTTTGATGACCGAGTTTCGCAGTGCGATATGCAAGAGTTACTACTTCATGTTCCTGTTCCTGGAGATCTGCATGCATCt                                                                   <  2:3699514/149‑1 (MQ=255)

                                                                                      gACGATGTGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGaa                                                               <  1:1494684/150‑1 (MQ=255)

                                                                                      gACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGaa                                                               >  2:4882162/1‑150 (MQ=255)

                                                                                      gACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGaa                                                               >  1:4430828/1‑150 (MQ=255)

                                                                                        cGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAgg                                                             >  1:1199323/1‑150 (MQ=255)

                                                                                        cGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAgg                                                             <  1:3543647/150‑1 (MQ=255)

                                                                                          aTGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGgc                                                           <  1:1194323/150‑1 (MQ=255)

                                                                                            gcgcACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGgcgc                                                         <  1:3596120/150‑1 (MQ=255)

                                                                                              gcACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCAGGGTGGGGttct                                                                                                                                            >  2:4978998/1‑63 (MQ=255)

                                                                                                  ccGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCaa                                                   <  1:1559856/150‑1 (MQ=255)

                                                                                                    gTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAgg                                                 >  1:4124384/1‑150 (MQ=255)

                                                                                                     tGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGc                                                <  2:1344408/150‑1 (MQ=255)

                                                                                                     tGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGc                                                <  2:1535220/150‑1 (MQ=255)

                                                                                                               tGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGcggc                                      <  2:37469/150‑1 (MQ=255)

                                                                                                               tGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGcggc                                      >  1:4543458/1‑150 (MQ=255)

                                                                                                                ggAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGcggcc                                     <  2:2821368/150‑1 (MQ=255)

                                                                                                                      catcTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGtt                               <  1:3273754/150‑1 (MQ=255)

                                                                                                                      catcTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGtt                               <  2:4928604/150‑1 (MQ=255)

                                                                                                                      catcTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGtt                               <  1:4683502/150‑1 (MQ=255)

                                                                                                                       atcTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGCATAcggcctt                                   >  1:4852266/1‑143 (MQ=255)

                                                                                                                        tcTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTgg                             <  2:4542302/150‑1 (MQ=255)

                                                                                                                        tcTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTgg                             <  1:4526430/150‑1 (MQ=255)

                                                                                                                             aCCCTGTCCGGCGAGAACACCACCCTGGGGTTGCGGACCGAGTTGCGCAGTGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCATCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGTCTGCg                        <  2:2807044/150‑1 (MQ=255)

                                                                                                                             aCCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCg                        <  2:412406/150‑1 (MQ=255)

                                                                                                                             aCCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCg                        <  1:4312632/150‑1 (MQ=255)

                                                                                                                               ccTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTAt                                           >  1:2759527/1‑129 (MQ=255)

                                                                                                                               ccTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTAt                                           <  2:2759527/129‑1 (MQ=255)

                                                                                                                               ccTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTAt                                           <  2:2759457/129‑1 (MQ=255)

                                                                                                                                 tGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCt                    <  1:3753585/150‑1 (MQ=255)

                                                                                                                                   gCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGa                  <  2:228163/150‑1 (MQ=255)

                                                                                                                                     cGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATc                <  2:2228910/150‑1 (MQ=255)

                                                                                                                                         cgcaaCCCCAGGGTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCtt            >  2:1126716/4‑150 (MQ=255)

                                                                                                                                          gcaaCCCCAGGGTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCttt           >  1:360771/3‑150 (MQ=255)

                                                                                                                                           gaACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTc          <  1:2815463/150‑1 (MQ=255)

                                                                                                                                           gaACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTc          <  1:501785/150‑1 (MQ=255)

                                                                                                                                           gaACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTc          <  2:4353121/150‑1 (MQ=255)

                                                                                                                                           gaACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTc          >  2:4574084/1‑150 (MQ=255)

                                                                                                                                           gaACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTc          <  2:4634803/150‑1 (MQ=255)

                                                                                                                                                ccACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTa     <  2:3775388/150‑1 (MQ=255)

                                                                                                                                                ccACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTa     <  1:4032549/150‑1 (MQ=255)

                                                                                                                                                ccACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTa     <  2:1181576/150‑1 (MQ=255)

TGTGCTCGACCGACACCTATGCGCAGGGGTTCCCTGAAGAAGGGCGCAAGCGGGCGATGAAGTGGCTGAAGATGCATGCCCAGTATGACGATGCGCACCCGTGGGAGGCGCTGGAAATCATCTGTACCCTGGCCGGCGAGAACCCCACCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCG  >  NC_002947/1966893‑1967189

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