BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,488,065	T→A	5.2%	K89* (AAG→TAG)	degS ←	serine endoprotease DegS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,488,065	0	T	A	5.2%	630.4 / 13.9	189	K89* (AAG→TAG)	degS	serine endoprotease DegS
Reads supporting (aligned to +/- strand): ref base T (80/98); new base A (7/3); total (87/102)
Fisher's exact test for biased strand distribution p-value = 1.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 205 total aligned reads displayed.

GTTGTTGGTCAGCAGGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGAAAACGCTGGA  >  NC_002947/1487916‑1488214

gttgttGGTCAGCAGGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGTGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCtt                                                                                                                                                       >  1:3341863/1‑150 (MQ=255)

 ttgttgGTCAGCAGGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTtg                                                                                                                                                      >  1:3529642/1‑150 (MQ=255)

 ttgttgGTCAGCAGGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTtg                                                                                                                                                      >  1:3835945/1‑150 (MQ=255)

 ttgttgGTCAGCAGGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTtg                                                                                                                                                      >  1:4031952/1‑150 (MQ=255)

 ttgttgGTCAGCAGGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTtg                                                                                                                                                      <  2:4895754/150‑1 (MQ=255)

   gttgGTCAGCAGGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTtggt                                                                                                                                                    >  1:426435/1‑150 (MQ=255)

    ttgGTCAGCAGGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTtggtg                                                                                                                                                   <  2:577592/150‑1 (MQ=255)

       gTCAGCAGGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCGGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGgtggtg                                                                                                                                                <  1:972196/150‑1 (MQ=255)

        tCAGCAGGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGgtgt                                                                                                                                               >  1:2480367/1‑150 (MQ=255)

           gcagGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACa                                                                                                                                            <  2:5141723/150‑1 (MQ=255)

            cagGTAGCCCTCTGGGCTCATGAGTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAg                                                                                                                                           <  2:5346938/150‑1 (MQ=255)

             agGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAgg                                                                                                                                          <  1:1480298/150‑1 (MQ=255)

                tAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTg                                                                                                                                       <  1:1997142/150‑1 (MQ=255)

                tAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTg                                                                                                                                       <  2:1903839/150‑1 (MQ=255)

                     ctctGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCAc                                                                                                                                  >  1:5156103/1‑150 (MQ=255)

                     ctctGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCAc                                                                                                                                  >  1:3248404/1‑150 (MQ=255)

                           gCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCgg                                                                                                                            >  1:1072789/1‑150 (MQ=255)

                             tCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCgggg                                                                                                                          <  2:4720918/150‑1 (MQ=255)

                             tCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCgggg                                                                                                                          >  1:2117797/1‑150 (MQ=255)

                              cATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGggc                                                                                                                         <  2:3759074/150‑1 (MQ=255)

                              cATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGCTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGggc                                                                                                                         >  2:1712038/1‑150 (MQ=255)

                               aTGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGggcg                                                                                                                        <  2:836248/150‑1 (MQ=255)

                                tGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGggcgg                                                                                                                       >  1:2742604/1‑150 (MQ=255)

                                 gATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGggcggc                                                                                                                      >  2:256926/1‑150 (MQ=255)

                                   ttACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCa                                                                                                                    <  2:1913382/150‑1 (MQ=255)

                                   ttACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCa                                                                                                                    >  1:2882693/1‑150 (MQ=255)

                                    tACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAg                                                                                                                   <  1:2560541/150‑1 (MQ=255)

                                    tACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAg                                                                                                                   <  2:684441/150‑1 (MQ=255)

                                     aCCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCGCCCCGGCGGTGGCCaac                                                                                                                  <  1:5215576/150‑3 (MQ=255)

                                     aCCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAgg                                                                                                                  <  2:2235023/150‑1 (MQ=255)

                                       cgccgAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACTGGTTGGCCGCCCCGGCGGTGGCCaacct                                                                                                                <  2:2378608/150‑5 (MQ=255)

                                       cgccgAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGt                                                                                                                >  2:4915012/1‑150 (MQ=255)

                                            aaCCCAGGCTTGATTCCCAGCGGCGCGGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCAcgg                                                                                                           >  1:1415203/1‑150 (MQ=255)

                                              cccAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCAcggcg                                                                                                         >  1:1970029/1‑150 (MQ=255)

                                               ccAGGCTTGATTCCCAGCGGCGGTGTTTGGGCAGGTTATCGCCGAAGTAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCAcggcgt                                                                                                        >  2:3857891/1‑150 (MQ=255)

                                               ccAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCAcggcgt                                                                                                        <  1:5235985/150‑1 (MQ=255)

                                               ccAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCAcggcgt                                                                                                        <  1:487870/150‑1 (MQ=255)

                                                cAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGAACaa                                                                                                                                           >  1:3358168/1‑113 (MQ=255)

                                                cAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGAACaa                                                                                                                                           <  2:3358168/114‑2 (MQ=255)

                                                 aGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTtggtg                                                                                                                                                   >  2:931610/1‑105 (MQ=255)

                                                 aGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCAcggcgtcg                                                                                                      <  2:4443162/150‑1 (MQ=255)

                                                   gCTTGAGTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGCAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCAcggcgtcggc                                                                                                    <  2:5198865/150‑1 (MQ=255)

                                                      tGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTa                                                                                                 <  1:1608036/150‑1 (MQ=255)

                                                       gATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTAc                                                                                                <  2:426435/150‑1 (MQ=255)

                                                             cAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACAcc                                                                                          >  1:5046214/1‑150 (MQ=255)

                                                              aGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCg                                                                                         >  2:3102234/1‑150 (MQ=255)

                                                               gcggcgCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCgg                                                                                        <  2:2145049/150‑1 (MQ=255)

                                                                 ggcgCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTTCGACACCGGGc                                                                                      <  2:2882693/150‑1 (MQ=255)

                                                                 ggcgCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGc                                                                                      >  2:2783710/1‑150 (MQ=255)

                                                                       gTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCa                                                                                <  1:4224246/150‑1 (MQ=255)

                                                                        tttGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCAt                                                                               >  1:5385027/1‑150 (MQ=255)

                                                                            ggCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATc                                                                           >  2:4080703/1‑150 (MQ=255)

                                                                             gCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATcc                                                                          >  2:5043059/1‑150 (MQ=255)

                                                                             gCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATcc                                                                          >  2:3244549/1‑150 (MQ=255)

                                                                              cAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCg                                                                         >  2:1509176/1‑150 (MQ=255)

                                                                               aGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTATGACACCGGGCCTTGCATGATCCgg                                                                        <  1:569884/150‑1 (MQ=255)

                                                                               aGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCgg                                                                        <  2:1187847/150‑1 (MQ=255)

                                                                               aGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCGTTTGTTCACCACCTTGGTGGTGTACGGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCgg                                                                        <  2:4780212/150‑1 (MQ=255)

                                                                                  ttATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTa                                                                     <  1:3708526/150‑1 (MQ=255)

                                                                                     tCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTc                                                                  >  2:1491952/1‑150 (MQ=255)

                                                                                           aagaagCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGgc                                                            <  1:1215296/150‑1 (MQ=255)

                                                                                                    cggAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTtgctgc                                                   >  2:3544441/1‑150 (MQ=255)

                                                                                                     ggAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTa                                                                                                 <  2:1585235/103‑1 (MQ=255)

                                                                                                      gAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAg                                                 <  1:3752768/150‑1 (MQ=255)

                                                                                                      gAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAg                                                 <  2:4689044/150‑1 (MQ=255)

                                                                                                      gAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAg                                                 <  2:3980908/150‑1 (MQ=255)

                                                                                                      gAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAg                                                 <  1:2769686/150‑1 (MQ=255)

                                                                                                         cTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGtt                                              <  1:2783710/150‑1 (MQ=255)

                                                                                                          tGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTg                                             >  2:1245990/1‑150 (MQ=255)

                                                                                                           gCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCAccgccg                                                                                                                             >  1:1325855/1‑69 (MQ=255)

                                                                                                             gggTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACg                                          >  2:451268/1‑150 (MQ=255)

                                                                                                               gTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTc                                        <  1:456843/150‑1 (MQ=255)

                                                                                                                tCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTcc                                       >  1:1528155/1‑150 (MQ=255)

                                                                                                                  ttCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCctg                                     >  2:5317538/1‑150 (MQ=255)

                                                                                                                   tCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCctgg                                    >  1:2999707/1‑150 (MQ=255)

                                                                                                                   tCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCctgg                                    <  1:4557077/150‑1 (MQ=255)

                                                                                                                    cGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAg                                                 >  2:3802442/1‑136 (MQ=255)

                                                                                                                      aaCAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCctggctg                                 >  1:1354471/1‑150 (MQ=255)

                                                                                                                        cAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTggg                               <  1:4670093/150‑1 (MQ=255)

                                                                                                                          gCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCa                             >  1:924802/1‑150 (MQ=255)

                                                                                                                          gCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCa                             <  2:1354471/150‑1 (MQ=255)

                                                                                                                                gggCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGa                       >  1:3630566/1‑150 (MQ=255)

                                                                                                                                gggCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGa                       >  2:629968/1‑150 (MQ=255)

                                                                                                                                 ggCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGc                                                                                      <  2:3115232/86‑1 (MQ=255)

                                                                                                                                 ggCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGc                                                                                      >  1:3115232/1‑86 (MQ=255)

                                                                                                                                 ggCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGAc                      <  2:1975743/150‑1 (MQ=255)

                                                                                                                                  gcgcTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCct                                      <  1:2305455/133‑1 (MQ=255)

                                                                                                                                  gcgcTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGAcc                     >  2:3648733/1‑150 (MQ=255)

                                                                                                                                    gcTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGAc                                           >  1:4493376/1‑126 (MQ=255)

                                                                                                                                     cTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCAt                  >  1:1305522/1‑150 (MQ=255)

                                                                                                                                      ttttGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCAtt                 >  1:1038081/1‑150 (MQ=255)

                                                                                                                                        ctGTACACCACCAAGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGc               >  2:2690852/2‑150 (MQ=255)

                                                                                                                                        ctGTACACCACCAAGGTGGTGAACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGc               >  1:3518986/2‑150 (MQ=255)

                                                                                                                                           ttCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGa            >  2:5380152/1‑150 (MQ=255)

                                                                                                                                           tacaccacCAAGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGa            >  2:391419/3‑150 (MQ=255)

                                                                                                                                            tCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGgcg                                                           >  2:455739/1‑102 (MQ=255)

                                                                                                                                                  ccTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGAAAACGCt     <  2:4031952/150‑1 (MQ=255)

                                                                                                                                                  ccTGGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGAAAACGCt     <  2:3341863/150‑1 (MQ=255)

                                                                                                                                                    ttGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGAAAACGCTgg   <  1:1899957/150‑1 (MQ=255)

                                                                                                                                                     tggtggTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGAAAACGCTGga  >  1:4012979/1‑150 (MQ=255)

GTTGTTGGTCAGCAGGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGAAAACGCTGGA  >  NC_002947/1487916‑1488214

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