BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,339,908	A→G	5.0%	K379R (AAG→AGG)	pcnB →	poly(A) polymerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,339,908	0	A	G	5.0%	295.5 / 3.5	99	K379R (AAG→AGG)	pcnB	poly(A) polymerase
Reads supporting (aligned to +/- strand): ref base A (45/49); new base G (2/3); total (47/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

Only 100 of 102 total aligned reads displayed.

CAGGGCGTGCCGCCGATCCCGGCCATGAACGGTGCGGCCCACGACCTGATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAGTGGTGGACCGATTACCAGGACGCAAACGACAGCGAGCGCCGCGAGAT  >  NC_002947/5339760‑5340055

cAGGGCGTGCCGCCGATCCCGGCCATGAACGGTGCGGCCCACGACCTGATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCaa                                                                                                                                                      <  3:215438/149‑1 (MQ=255)

 aGGGCGTGCCGCCGATCCCGGCCATGAACGGTGCGGCCCACGACCTGATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAg                                                                                                                                                     <  2:356839/149‑1 (MQ=255)

    gCGTGCCGCCGATCCCGGCCATGAACGGTGCGGCCCACGACCTAATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCTTGCCGCTGCGGCGCGGCAGgcgc                                                                                                                                                  <  6:263888/149‑1 (MQ=255)

     cGTGCCGCCGATCCCGGCCATGAACGGTGCGGCCCACGACCTGATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAgcgc‑g                                                                                                                                                >  1:431645/1‑149 (MQ=255)

        gccgccGATCCCGGCCATGAACGGTGCGGCCCACGACCTGATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑Gcc                                                                                                                                              >  2:359911/1‑148 (MQ=255)

        gccgccGATCCCGGCCATGAACGGTGCGGCCCACGACCTGATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑Gcc                                                                                                                                              >  7:84150/1‑148 (MQ=255)

         ccgccgATCCCGGCCATGAACGGTGCGGCCCACGACCTGATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGa                                                                                                                                            >  6:14549/1‑149 (MQ=255)

         ccgccgATCCCGGCCATGAACGGTGCGGCCCACGACCTGATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTCCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCTTGCCGCTGCGGCGCGGCAGGCGC‑Tcct                                                                                                                                             <  4:152096/148‑2 (MQ=255)

                  ccGGCCATGAACGGTGCGGCCCACGACCTGATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATgctgct                                                                                                                                   <  3:160267/149‑1 (MQ=255)

                  ccGGCCATGAACGGTGCGGCCCACGACCTGATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATgctgct                                                                                                                                   <  7:190804/149‑1 (MQ=255)

                   cGGCCATGAACGGTGCGGCCCACGACCTGATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGTACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGctgctg                                                                                                                                  >  2:143268/1‑149 (MQ=255)

                   cGGCCATGAACGGTGCGGCCCACGACCTGATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGAAACGGAAGGCGCGCCTGCCGCTCAGCAGCGGCAAGCGC‑GCCGACAtgctccgg                                                                                                                                  >  4:364030/1‑145 (MQ=255)

                   cGGCCATGAACGGTGCGGCCCACGACCTGATCGCCGAGCAGTGCGCGCGCATTGCCATCACCAAGCGATTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGAGCCTGCCGCTCCGCAGCTGCCAGCGC‑CCCGACATgctcctc                                                                                                                                  >  6:253706/1‑148 (MQ=255)

                    ggCCATGAACGGTGCGGCCCACGACCTGATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTgg                                                                                                                                 <  7:348333/149‑1 (MQ=255)

                       cATGAACGGTGCGGCCCACGACCTGATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGAc                                                                                                                               >  6:135301/1‑148 (MQ=255)

                        aTGAACGGTGCGGCCCACGACCTGATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACaa                                                                                                                             >  4:306865/1‑149 (MQ=255)

                         tGAACGGTGCGGCCCACGACCTGATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAAc                                                                                                                            <  4:120866/149‑1 (MQ=255)

                             cGGTGCGGCCCACGACCTGATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCcgc                                                                                                                        <  6:305385/149‑1 (MQ=255)

                                       cACGACCTGATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCgcgc                                                                                                              >  5:333694/1‑149 (MQ=255)

                                        aCGACCTGATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGcc                                                                                                             <  7:406048/149‑1 (MQ=255)

                                             cTGATCGCTGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGATA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTa                                                                                                        <  1:161728/149‑1 (MQ=255)

                                                  cGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACt                                                                                                   >  6:254579/1‑149 (MQ=255)

                                                   gCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACt                                                                                                   <  6:234004/148‑1 (MQ=255)

                                                    ccGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCcg                                                                                                                 <  7:214705/133‑1 (MQ=255)

                                                    ccGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCcg                                                                                                                 >  8:214705/1‑133 (MQ=255)

                                                       aGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCctg                                                                                              >  3:150604/1‑149 (MQ=255)

                                                         cAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCctgct                                                                                            >  5:302985/1‑149 (MQ=255)

                                                         cAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCctgct                                                                                            >  8:175228/1‑149 (MQ=255)

                                                          aGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCctgctg                                                                                           >  5:276867/1‑149 (MQ=255)

                                                            tGCGCGCGCATTGCCATCCCCAAGCGGTTTACCCTGCCGATCCGCGAGATCTCTGCTACATGTACCAATAGCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCctgctg                                                                                           <  1:359911/149‑1 (MQ=255)

                                                              cgcgcgcATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGATA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCctgctgctg                                                                                        <  1:143268/148‑1 (MQ=255)

                                                               gcgcacATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCcgc                                                                                                                >  3:247717/1‑123 (MQ=255)

                                                               gcgcacATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCcgc                                                                                                                <  4:247717/123‑1 (MQ=255)

                                                                   gcATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGaa                                                                                  <  5:135301/149‑1 (MQ=255)

                                                                    cATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGaaa                                                                                 <  2:340864/149‑1 (MQ=255)

                                                                    cATTGCCATCCCCAAGCGGTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGaaa                                                                                 <  4:28672/149‑1 (MQ=255)

                                                                      ttGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑CCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAgc                                                                               >  5:237193/1‑149 (MQ=255)

                                                                          cATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCg                                                                           >  3:167931/1‑149 (MQ=255)

                                                                           aTCCCCAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTTGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCgg                                                                          <  6:433779/149‑1 (MQ=255)

                                                                              cccAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCgggg                                                                        <  7:69988/148‑1 (MQ=255)

                                                                               ccAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTAc                                                                                                       >  8:337721/1‑116 (MQ=255)

                                                                               ccAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTAc                                                                                                       <  7:337721/116‑1 (MQ=255)

                                                                                cAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCtgctgctgc                                                                                       >  6:315446/1‑131 (MQ=255)

                                                                                cAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCtgctgctgc                                                                                       <  5:315446/131‑1 (MQ=255)

                                                                                cAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCgggg                                                                        >  8:371586/1‑146 (MQ=255)

                                                                                 actcGTATTACCCTGCCGATCAGCGAGA‑‑GCTGTGACATGCAGTAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGgaaga                                                                    <  8:270175/146‑1 (MQ=255)

                                                                                 aaGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGa                                                                                                                                >  8:171856/1‑89 (MQ=255)

                                                                                 aaGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGa                                                                                                                                <  7:171856/89‑1 (MQ=255)

                                                                                 aaGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGgaaga                                                                    <  2:329362/149‑1 (MQ=255)

                                                                                   gCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGaaa                                                                  >  3:82639/1‑149 (MQ=255)

                                                                                   gCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGaaa                                                                  >  4:361196/1‑149 (MQ=255)

                                                                                     gTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAAcc                                                                <  3:361196/149‑1 (MQ=255)

                                                                                              tGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTg                                                       >  6:50473/1‑149 (MQ=255)

                                                                                               gCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGctgctg                                                                                                                                  >  3:380771/1‑73 (MQ=255)

                                                                                               gCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGctgctg                                                                                                                                  <  4:380771/73‑1 (MQ=255)

                                                                                                ccGATCCGCGAGA‑‑TCTGGGACATGCAGGAGAGCCTGCCGCGCCGAAGCGGCAAGAGC‑TCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTggg                                                     >  3:144486/1‑149 (MQ=255)

                                                                                                ccGATCCGCGAGA‑‑TCTAGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGTCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTgg                                                      <  7:229502/148‑1 (MQ=255)

                                                                                                 cGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGc                                                    >  7:125975/1‑149 (MQ=255)

                                                                                                  caTCTGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGc                                                    <  5:418735/147‑1 (MQ=255)

                                                                                                   aTCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGGGGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGGTGGACAACCCGCGTTTCCGCGGCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCa                                                  <  3:208579/149‑1 (MQ=255)

                                                                                                   aTCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCa                                                  <  5:254579/149‑1 (MQ=255)

                                                                                                      cgcgAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAgtg                                               <  4:82639/149‑1 (MQ=255)

                                                                                                       gcgAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAgtgg                                              >  7:98301/1‑149 (MQ=255)

                                                                                                       gcgAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAgtgg                                              <  6:276867/149‑1 (MQ=255)

                                                                                                       gcgAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAgtgg                                              >  2:346012/1‑149 (MQ=255)

                                                                                                       gcgAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAgtgg                                              >  1:74396/1‑149 (MQ=255)

                                                                                                               tCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAGTGGTGGAcc                                        >  2:134040/1‑149 (MQ=255)

                                                                                                               tCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAGTGGTGGAcc                                        >  7:267680/1‑149 (MQ=255)

                                                                                                                   ggACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAGTGGTGGACCGa                                      <  4:413619/147‑1 (MQ=255)

                                                                                                                   ggACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAGTGGTGGACCGAt                                     >  3:333403/1‑148 (MQ=255)

                                                                                                                   ggACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAGTGGTGGACCGAt                                     <  6:333694/148‑1 (MQ=255)

                                                                                                                    gACATGCAGGAGCGCCTGCCGGGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAGTGGTGGACCGAtt                                    <  6:302985/148‑1 (MQ=255)

                                                                                                                     aCATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑CCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAGTGGTGGACCGATTAc                                  <  6:237193/149‑1 (MQ=255)

                                                                                                                       aTGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAGTGGTGGACCGATTACCa                                <  7:175228/149‑1 (MQ=255)

                                                                                                                       aTGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAACAAACCGACGAACTGGGCCAGTGGTGGACCGATTACCa                                <  3:113469/149‑1 (MQ=255)

                                                                                                                         gCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGCACTGGGCCAGTGGTGGACCGATTACCAg                               >  2:310811/1‑148 (MQ=255)

                                                                                                                         gCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAgtggtg                                            <  3:247570/135‑1 (MQ=255)

                                                                                                                         gCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAgtggtg                                            >  4:247570/1‑135 (MQ=255)

                                                                                                                         gCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAGTGGTGGACCGATTACCAg                               >  6:282268/1‑148 (MQ=255)

                                                                                                                         gCAGGAGCGCCTGCCGCGCCGCAGAGGCAAGCGC‑GCAGACATGATGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCAGGGGAAGAACCCGACGAACTGGGCCAGTGGTGGACCGATTACCAgg                              >  7:253775/1‑149 (MQ=255)

                                                                                                                          cAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAGTGGTGGACCGATTACCAGGa                             <  7:371586/149‑1 (MQ=255)

                                                                                                                           aagAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAGTGGTGGACCGATTACCAGGAc                            <  3:427246/147‑1 (MQ=255)

                                                                                                                           aGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAGTGGTGGACCGATTACCAGGAc                            <  1:134040/149‑1 (MQ=255)

                                                                                                                            ggcgcgcTTGCCGCTGCGGCGCGGCAGGCGC‑TCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAGTGGTGGACCGATTACCAGGACg                           <  1:145729/146‑1 (MQ=255)

                                                                                                                             gAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAGTGGTGGACCGATTACCAGGACGc                          <  5:14549/149‑1 (MQ=255)

                                                                                                                             gAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAGTGGTGGACCGATTACCAGGACGc                          <  3:306865/149‑1 (MQ=255)

                                                                                                                               gcgcTTGCCGCTGCGGCGCGGCAGGCGCTGCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAGTGGTGGACCGATTACCAgg                              >  1:247947/1‑144 (MQ=255)

                                                                                                                                cgcTTGCCGCTGCGGCGCGGCAGGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAGTGGTGGACCGATTACCAGGACGCaaa                       >  1:161311/1‑149 (MQ=255)

                                                                                                                                  ccTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAGTGGTGGACCGATTACCAGGACGCAAAc                      >  5:338658/1‑148 (MQ=255)

                                                                                                                                   cTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAGTGGTGGACCGATTACCAGGACGCAAACGa                    >  7:123292/1‑149 (MQ=255)

                                                                                                                                   cTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATCCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGGGAAAGCGCCGGGGAAGAAACCGACGACCTGGGCCACTGGTGGACCGAGAACCAGGATGCAAACgc                    >  7:107547/1‑148 (MQ=255)

                                                                                                                                      ccgcgcctcAGCGTCAAGCGC‑TCCGACATGCTGCTGGACAAGCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGTGAAGAAACCGACTAACTGGGCCAGTGGTGGACCGATTACCAGGACGCAAACGACa                  <  8:253775/148‑1 (MQ=255)

                                                                                                                                      ccgcgccgcAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCTGTGGTGGACCGATTACCAGGACGCAAACGACa                  <  3:390401/148‑1 (MQ=255)

                                                                                                                                           ccgcAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCtgctgctgc                                                                                       <  4:4866/74‑1 (MQ=255)

                                                                                                                                           ccgcAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCtgctgctgc                                                                                       >  3:4866/1‑74 (MQ=255)

                                                                                                                                              cAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAGTGGTGGACCGATTACCAGGACGCAAACGACAGCGAGcgcc         >  2:43577/1‑149 (MQ=255)

                                                                                                                                                gCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAGTGGTGGACCGATTACCAGGACGCAAACGACAGCGAGcgccgc       <  6:338658/149‑1 (MQ=255)

                                                                                                                                                   gCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAGTGGTGGACCGATTACCAGGACGCAAACGACAGCGAGcgccgc       >  3:87957/1‑146 (MQ=255)

                                                                                                                                                     aaGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAGTGGTGGACCGATTACCAGGACGCAAACGACAGCGAGCGCCGCgaga   <  4:87957/148‑1 (MQ=255)

                                                                                                                                                     aaGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAGTGGTGGACCGATTACCAGGACGCAAACGACAGCGAGCGCCGCGAgat  >  2:271674/1‑149 (MQ=255)

CAGGGCGTGCCGCCGATCCCGGCCATGAACGGTGCGGCCCACGACCTGATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAGTGGTGGACCGATTACCAGGACGCAAACGACAGCGAGCGCCGCGAGAT  >  NC_002947/5339760‑5340055

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 21 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`