BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,784,074	C→G	6.1%	Q15H (CAG→CAC) ‡	glnD ←	uridylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,784,074	0	C	G	6.1%	558.9 / 16.9	170	Q15H (CAG→CAC) ‡	glnD	uridylyltransferase
Reads supporting (aligned to +/- strand): ref base C (85/74); new base G (4/7); total (89/81)
Fisher's exact test for biased strand distribution p-value = 3.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 173 total aligned reads displayed.

GCAGGATATTGTCGACGAGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGATCTCGTAGCCAGTGGCG  >  NC_002947/1783926‑1784223

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

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

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

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

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

       atTGTCGACGAGCCATGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTg                                                                                                                                               <  2:1800290/150‑1 (MQ=255)

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

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

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

                 aGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCg                                                                                                                                     >  2:4018439/1‑150 (MQ=255)

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

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

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

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

                           gggCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCgg                                                                                                                           <  1:4181701/150‑1 (MQ=255)

                             gCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCgggg                                                                                                                         <  2:2495359/150‑1 (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                                                                                                                <  2:3957139/150‑1 (MQ=255)

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

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

                                        cGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTgggg                                                                                                              >  1:3388065/1‑150 (MQ=255)

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

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

                                                              cgAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCgg                                                                                        >  2:3831469/1‑149 (MQ=255)

                                                                aaaaCGTTTGTCGAGCACCTCGGCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTgg                                                                                      >  1:3873696/1‑150 (MQ=255)

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

                                                                    cGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTg                                                                                  <  2:4449141/150‑1 (MQ=255)

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

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

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

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

                                                                            gAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTgg                                                                                 >  2:2398280/1‑143 (MQ=255)

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

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

                                                                            gAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCg                                                                          <  2:608893/150‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                                                                   tCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCAGGCGGAACTGGcc                                                                                                                                            <  2:1680273/77‑1 (MQ=38)

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

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

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

                                                                                            aGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCg                                                          >  2:4174010/1‑150 (MQ=255)

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

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

                                                                                              gCGGATGGCTGTCTTGAGGGCGGCGCTGGGGCTCGCCTTGAGGGCCGGTTGCGCCTGGCACTGGCCACGGTCGACCCGCTCGGGGTCCACCTGGGGCATCGCGTCACGTGCCTGTGGTGGGGTGGCGTGTGGGGGTGCGGTCGGGCCggg                                                        >  2:4153358/1‑148 (MQ=255)

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

                                                                                                    ggCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCggg                                                                        <  2:4376558/128‑1 (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)

                                                                                                       tttCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCggg                                               >  2:3863528/1‑150 (MQ=255)

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

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

                                                                                                          cTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGGTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATg                                            <  2:4868124/150‑1 (MQ=255)

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

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

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

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

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

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

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

                                                                                                            tGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGgt                                          <  1:4018439/150‑1 (MQ=255)

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

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

                                                                                                                  cggcgATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTgtcgtc                                    <  1:3014946/150‑1 (MQ=255)

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

                                                                                                                   ggcgATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTcc                                   >  2:3891140/1‑150 (MQ=255)

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

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

                                                                                                                           ggCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAg                           >  1:436699/1‑150 (MQ=255)

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

                                                                                                                            gCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAgg                          >  1:4420679/1‑150 (MQ=255)

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

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

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

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

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

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

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

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

                                                                                                                                       gggCCAGTTCCGCCGGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGAtctc               >  2:296858/1‑150 (MQ=255)

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

                                                                                                                                        ggCCAGTTCCGCCCGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGATCTCg              <  1:296858/150‑1 (MQ=255)

                                                                                                                                         gCCAGTTCCAGGCGGAACTGGCCACGGCCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGATCTCGt             >  1:3282899/1‑150 (MQ=255)

                                                                                                                                           cAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGATCTCGTAg           >  2:1086692/1‑150 (MQ=255)

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

                                                                                                                                                 cGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGATCTCGTAGCCAGTg     >  1:4549798/1‑150 (MQ=255)

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

GCAGGATATTGTCGACGAGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGATCTCGTAGCCAGTGGCG  >  NC_002947/1783926‑1784223

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