BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,339,900	C→G	6.7%	R376R (CGC→CGG)	pcnB →	poly(A) polymerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,339,900	0	C	G	6.7%	316.1 / 8.0	99	R376R (CGC→CGG)	pcnB	poly(A) polymerase
Reads supporting (aligned to +/- strand): ref base C (42/49); new base G (2/4); total (45/54)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.14e-01

Only 100 of 104 total aligned reads displayed.

AGAACCAGGGCGTGCCGCCGATCCCGGCCATGAACGGTGCGGCCCACGACCTGATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAGTGGTGGACCGATTACCAGGACGCAAACGACAGCGAGCGCC  >  NC_002947/5339755‑5340048

aGAACCAGGGCGTGCCGCCGATCCCGGCCATGAACGGTGCGGCCCACGACCTGATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGc                                                                                                                                                    <  6:409896/149‑1 (MQ=255)

 gAACCAGGGCGTGCCGCCGATCCCGGCCATGAACGGTGCGGCCCACGACCTGATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCg                                                                                                                                                   <  4:397350/149‑1 (MQ=255)

 gAACCAGGCCGGGCCGCCGATCCCGGCCATGAACGGCGCGGCCCACGTCCTGATCGCCGAGCAGGGCGCGCGCATTGCCATCCCCCCGGGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCg                                                                                                                                                   <  5:232181/149‑1 (MQ=255)

  aaCCAGGGCGTGCCGCCGATCCCGGCCATGAACGGTGCGGCCCACGACCTGATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCgg                                                                                                                                                  <  8:77348/149‑1 (MQ=255)

  aaCCAGGGCGTGCCGCCGATCCCGGCCATGAACGGTGCGGCCCACGACCTGATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCgg                                                                                                                                                  >  1:356839/1‑149 (MQ=255)

     cAGGGCGTGCCGCCGATCCCGGCCATGAACGGTGCGGCCCACGACCTGATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCTGCGGCGCGGCag                                                                                                                                               <  3:281934/149‑2 (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                                                                                                                                       >  7:84150/1‑148 (MQ=255)

             gccgccGATCCCGGCCATGAACGGTGCGGCCCACGACCTGATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑Gcc                                                                                                                                       >  2:359911/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                                                                                                                            <  7:190804/149‑1 (MQ=255)

                       ccGGCCATGAACGGTGCGGCCCACGACCTGATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATgctgct                                                                                                                            <  3:160267/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)

                                    gTGCGGCCCACGACCTGATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGTAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGt                                                                                                               <  5:131670/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                                                                                                          >  8:214705/1‑133 (MQ=255)

                                                         ccGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCcg                                                                                                          <  7:214705/133‑1 (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                                                                                                         <  4:247717/123‑1 (MQ=255)

                                                                    gcgcacATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCcgc                                                                                                         >  3:247717/1‑123 (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                                                                                <  5:315446/131‑1 (MQ=255)

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

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

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

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

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

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

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

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

                                                                                          gTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAAcc                                                         <  3:361196/149‑1 (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                                       >  2:346012/1‑149 (MQ=255)

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

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

                                                                                                            gcgAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAgtgg                                       <  6:276867/149‑1 (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)

AGAACCAGGGCGTGCCGCCGATCCCGGCCATGAACGGTGCGGCCCACGACCTGATCGCCGAGCAGTGCGCGCGCATTGCCATCCCCAAGCGTTTTACCCTGCCGATCCGCGAGA‑‑TCTGGGACATGCAGGAGCGCCTGCCGCGCCGCAGCGGCAAGCGC‑GCCGACATGCTGCTGGACAACCCGCGTTTCCGCGCCGGTTACGACTTCCTGCTGCTGCGTGAAAGCGCCGGGGAAGAAACCGACGAACTGGGCCAGTGGTGGACCGATTACCAGGACGCAAACGACAGCGAGCGCC  >  NC_002947/5339755‑5340048

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