BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,784,075	T→C	6.4%	Q15R (CAG→CGG) ‡	glnD ←	uridylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,784,075	0	T	C	6.4%	535.0 / 15.3	171	Q15R (CAG→CGG) ‡	glnD	uridylyltransferase
Reads supporting (aligned to +/- strand): ref base T (84/75); new base C (4/7); total (89/82)
Fisher's exact test for biased strand distribution p-value = 3.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 178 total aligned reads displayed.

GCAGGATATTGTCGACGAGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGATCTCGTAGCCAGTGGCGG  >  NC_002947/1783926‑1784224

gCAGGATATTGTCGACGAGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCt                                                                                                                                                       >  1:4595062/1‑150 (MQ=255)

 cAGGATATTGTCGACGAGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCGGTTCCGCCTg                                                                                                                                                      >  1:643852/1‑150 (MQ=255)

 cAGGATATTGTCGACGAGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTg                                                                                                                                                      >  2:2905235/1‑150 (MQ=255)

  aGGATATTGTCGACGAGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTgg                                                                                                                                                     >  2:3822104/1‑150 (MQ=255)

       atTGTCGACGAGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTg                                                                                                                                                <  2:652325/150‑1 (MQ=255)

        ttGTCGACGAGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCAGGCGGAACTgg                                                                                                                                               <  2:1138049/150‑1 (MQ=255)

        ttGTCGACGAGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCAGGCGGAACTgg                                                                                                                                               <  2:3855239/150‑1 (MQ=255)

          gTCGACGAGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCAGGCGGAACTGGcc                                                                                                                                             <  2:1145831/150‑1 (MQ=255)

           tCGACGAGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCa                                                                                                                                            >  1:4239439/1‑150 (MQ=255)

              acgaGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACgg                                                                                                                                         >  1:4876026/1‑150 (MQ=255)

                gaGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTc                                                                                                                                       >  2:3914009/1‑150 (MQ=255)

                  gCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGa                                                                                                                                     >  2:2606417/1‑150 (MQ=255)

                   ccAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAACCGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGaa                                                                                                                                    <  2:5419522/150‑1 (MQ=255)

                   ccAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGaa                                                                                                                                    <  2:4436573/150‑1 (MQ=255)

                   ccAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGAGGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGaa                                                                                                                                    >  2:4170566/1‑150 (MQ=255)

                    cAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAAc                                                                                                                                   <  2:1752522/150‑1 (MQ=255)

                     aGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGCGGGCTTGCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACa                                                                                                                                  >  1:5107967/1‑150 (MQ=255)

                     aGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACa                                                                                                                                  >  2:3702552/1‑150 (MQ=255)

                           gggCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCgg                                                                                                                            >  1:4357544/1‑150 (MQ=255)

                             gCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCgggg                                                                                                                          <  2:2495359/150‑1 (MQ=255)

                              ccTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGt                                                                                                                         >  2:5081470/1‑150 (MQ=255)

                               cTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTc                                                                                                                        >  1:2795659/1‑150 (MQ=255)

                                  gATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCAc                                                                                                                     >  2:3014946/1‑150 (MQ=255)

                                    tCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCt                                                                                                                   <  1:3717386/150‑1 (MQ=255)

                                     cAGCGGTGGGTTTTCGCTGCCGGCGCGAAAACGTTTGGCGAGCACCTCGGCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTg                                                                                                                  <  1:1380797/150‑1 (MQ=255)

                                      aGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTgg                                                                                                                 <  1:5347440/150‑1 (MQ=255)

                                      aGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTgg                                                                                                                 <  1:902336/150‑1 (MQ=255)

                                       gCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTggg                                                                                                                <  2:5019490/150‑1 (MQ=255)

                                                         cGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCt                                                                                              <  2:4840926/150‑1 (MQ=255)

                                                          ggCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTg                                                                                             >  1:3120353/1‑150 (MQ=255)

                                                           gcgcGAAACCGTGGCCCGAGCCCCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCGTGCGGGCCAGGTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGt                                                                                            <  2:4620589/150‑1 (MQ=255)

                                                             gcgAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCg                                                                                          >  1:2903645/1‑150 (MQ=255)

                                                                 aaaCGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTggg                                                                                      >  1:685999/1‑150 (MQ=255)

                                                                   aCGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGt                                                                                    <  2:1660626/150‑1 (MQ=255)

                                                                     gTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTgg                                                                                  >  1:4467396/1‑150 (MQ=255)

                                                                     gTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTgg                                                                                  >  1:4979677/1‑150 (MQ=255)

                                                                      tttGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGc                                                                                 >  1:1809219/1‑150 (MQ=255)

                                                                         gTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTgg                                                                                       <  1:1272229/141‑1 (MQ=255)

                                                                         gTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTgg                                                                                       >  2:1272229/1‑141 (MQ=255)

                                                                            gAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCg                                                                           >  2:799286/1‑150 (MQ=255)

                                                                            gAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCg                                                                           >  2:1532159/1‑150 (MQ=255)

                                                                            gAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCg                                                                           >  1:113847/1‑150 (MQ=255)

                                                                            gAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCg                                                                           >  1:4736932/1‑150 (MQ=255)

                                                                            gAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCAGGCGGAACTGGcc                                                                                                                                             <  2:2099074/84‑1 (MQ=38)

                                                                                 ccTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTg                                                                      >  1:1609199/1‑150 (MQ=255)

                                                                                  cTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGt                                                                                    >  1:5309031/1‑135 (MQ=255)

                                                                                  cTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGt                                                                                    <  2:5309031/135‑1 (MQ=255)

                                                                                  cTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGc                                                                     >  1:5151277/1‑150 (MQ=255)

                                                                                   tCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCg                                                                    <  1:291555/150‑1 (MQ=255)

                                                                                    cGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCgg                                                                   <  1:3822104/150‑1 (MQ=255)

                                                                                     gCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGt                                                                  >  2:2259140/1‑150 (MQ=255)

                                                                                        ggCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAg                                                               <  1:4816313/150‑1 (MQ=255)

                                                                                        ggCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAg                                                               >  1:2343135/1‑150 (MQ=255)

                                                                                              gCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGaa                                                         <  1:4152568/150‑1 (MQ=255)

                                                                                              gCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGaa                                                         >  2:4153240/1‑150 (MQ=255)

                                                                                               cGGCTGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGACCTGGCCACGGTCGACCAGCTCGGGGTCCACCTGGGGCATCGCGTCACGGTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAg                                                        <  2:4088904/150‑1 (MQ=255)

                                                                                                ggATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGt                                                       >  2:922627/1‑150 (MQ=255)

                                                                                                 gATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTg                                                      >  1:4008157/1‑150 (MQ=255)

                                                                                                    ggCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCggg                                                                         >  1:4376558/1‑128 (MQ=255)

                                                                                                    ggCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTgcgc                                                   <  2:1855453/150‑1 (MQ=255)

                                                                                                    ggCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTgcgc                                                   <  2:5328155/150‑1 (MQ=255)

                                                                                                       tttCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTc                                                                                                >  1:2614674/1‑102 (MQ=255)

                                                                                                       tttCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTc                                                                                                <  2:2614674/102‑1 (MQ=255)

                                                                                                        ttCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGa                                               <  1:4094901/150‑1 (MQ=255)

                                                                                                          cTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATg                                             <  2:4868367/150‑1 (MQ=255)

                                                                                                          cTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATg                                             <  2:4922134/150‑1 (MQ=255)

                                                                                                          cTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATg                                             <  2:4008157/150‑1 (MQ=255)

                                                                                                          cTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATg                                             <  1:3863528/150‑1 (MQ=255)

                                                                                                           ttGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCGTGTCGGGGTGCGGTCGGGCCGAAGTGCGCGGGAggg                                            >  2:3609402/1‑150 (MQ=255)

                                                                                                           ttGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGc                                                                     >  1:3210728/1‑125 (MQ=255)

                                                                                                           ttGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGc                                                                     <  2:3210728/125‑1 (MQ=255)

                                                                                                           ttGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATgg                                            >  2:3609413/1‑150 (MQ=255)

                                                                                                           ttGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATgg                                            >  2:1324198/1‑150 (MQ=255)

                                                                                                            tGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGgt                                           >  2:2688008/1‑150 (MQ=255)

                                                                                                               aGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTgtc                                        <  1:798691/150‑1 (MQ=255)

                                                                                                               aGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTgtc                                        <  1:2211257/150‑1 (MQ=255)

                                                                                                                 gcggcgATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTgtcgt                                      <  1:98551/150‑1 (MQ=255)

                                                                                                                   ggcgATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTgtcgtca                                    >  2:2954995/1‑149 (MQ=255)

                                                                                                                     cgATGGGGCTCGGCTTGAGGGCCAGTTCCGCGTGGAACTGGCCACGGTCGAACAGCTCGGGGTCGACCTGGGGCTTCGCGTCTCGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTGGTCCtt                                  <  1:1165398/150‑1 (MQ=255)

                                                                                                                     cgATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCtt                                  <  2:3388065/150‑1 (MQ=255)

                                                                                                                            gCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCTGGCCGAAGTGCGCGGGATGGTTTCGTCCTTGCGCAgg                           <  2:4538086/150‑1 (MQ=255)

                                                                                                                            gCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAgg                           <  1:2011098/150‑1 (MQ=255)

                                                                                                                            gCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAgg                           <  2:2373783/150‑1 (MQ=255)

                                                                                                                             cTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGt                                                                                    <  2:2983635/92‑1 (MQ=255)

                                                                                                                             cTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGt                                                                                    >  1:2983635/1‑92 (MQ=255)

                                                                                                                              tCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGt                         >  1:3144570/1‑150 (MQ=255)

                                                                                                                                gCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGa                       >  1:3844526/1‑150 (MQ=255)

                                                                                                                                gCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGa                       <  1:3633463/150‑1 (MQ=255)

                                                                                                                                 ccTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGaa                      >  1:1145244/1‑150 (MQ=255)

                                                                                                                                 ccTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGaa                      <  1:3914009/150‑1 (MQ=255)

                                                                                                                                    tGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGAt                   <  2:2543099/150‑1 (MQ=255)

                                                                                                                                     cgtggCCAGTTCCAGGCGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGAtc                  >  1:338306/4‑150 (MQ=255)

                                                                                                                                        ggCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGATCTCg               <  2:3882767/150‑1 (MQ=255)

                                                                                                                                        ggCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGATCTCg               <  2:3882329/150‑1 (MQ=255)

                                                                                                                                              ttCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGt                                                                                         <  2:599828/70‑1 (MQ=255)

                                                                                                                                              ttCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGt                                                                                         >  1:599828/1‑70 (MQ=255)

                                                                                                                                                    cTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGATCTCGTAGCCAGTGGCg   >  2:435210/1‑150 (MQ=255)

                                                                                                                                                    cTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGATCTCGTAGCCAGTGGCg   <  1:2606417/150‑1 (MQ=255)

                                                                                                                                                     tGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGATCTCGTAGCCAGTGGCgg  >  1:2305080/1‑150 (MQ=255)

                                                                                                                                                     tGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGACCGAGGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGATCTCGTAGCCAGTGGCgg  <  1:1333133/150‑1 (MQ=255)

GCAGGATATTGTCGACGAGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGATCTCGTAGCCAGTGGCGG  >  NC_002947/1783926‑1784224

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