BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,299,727	T→C	48.1%	Y132C (TAC→TGC)	sbcC ←	exonuclease subunit SbcC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,299,727	0	T	C	48.1%	0.6 / 142.8	104	Y132C (TAC→TGC)	sbcC	exonuclease subunit SbcC
Reads supporting (aligned to +/- strand): ref base T (33/21); new base C (26/24); total (59/45)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-01

Only 100 of 105 total aligned reads displayed.

TCGTCGGCCTTGAGGAACGCGCCGAACTCGCTCTGGGCCAGCATCACCGCACGGGTGAACTGTTCGAAGTTCAGGCCCAGACGGGCTTCGACCAGTTGCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGTAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCCGACGAAATCGACCTCGGCAAAACCGCTGCCGGTGCCACGGCGCAGCAGGTTGCG  >  NC_002947/2299580‑2299875

tcgtcgGCCTTGAGGAACGCGCCTAACTCGCTCTGGGCCAGCATCACCGCACGGGTGAACTGTTCTAAGTTCAGGCCCAGACGTGCTTCGACCAGTTGCTTGTTTTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGTa                                                                                                                                                     <  7:121844/149‑1 (MQ=255)

  gtcgGCCTTGAGGAACGCGCCGAACTCGCTCTGGGCCAGCATCACCGCACGGGTGAACTGTTCGAAGTTCAGGCCCAGACGGGCTTCGACCAGTTGCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGCAGa                                                                                                                                                   <  1:149199/149‑1 (MQ=255)

  gtcgGCCTTGAGGAACGCGCCGAACTCGCTCTGGGCCAGCATCACCGCACGGGTGAACTGTTCGAAGTTCAGGCCCAGACGGGCTTCGACCAGTTGCTTGTATTCGTTCTTGCCACTCCCGAGCACTTGCTCGCTGTCGAGGTCGTAg                                                                                                                                                    <  1:30455/148‑1 (MQ=255)

  gtcgGCCTTGAGGAACGCGCCGAACTCGCTCTGGGCCAGCATCACCGCACGGGTGAACTGTTCGAAGTTCAAGCCCAGACGGGCTTCGACCAGTTGCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGTAg                                                                                                                                                    >  4:143407/1‑148 (MQ=255)

    cgGCCTTGAGGAACGCGCCGAACTCGCTCTGGGCCAGCATCACCGCACGGGTGAACTGTTCGAAGTTCAGGCCCAGACGGGCTTCGACCAGTTGCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGCAGAAg                                                                                                                                                 >  7:287786/1‑149 (MQ=255)

           gAGGAACGCGCCGAACTCGCTCTGGGCCAGCATCACCGCACGGGTGAACTGTTCGAAGTTCAGGCCCAGACGGGCTTCGACCAGTTGCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGTAGAAGCTTTGa                                                                                                                                           <  6:32612/148‑1 (MQ=255)

           gAGGAACGCGCCGAACTCGCTCTGGGCCAGCATCACCGCACGGGTGAACTGTTCGAAGTTCAGGCCCAGACGGGCTTCGACCAGATGCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCCCTGTCGAGGTCGTAGAAGCTTTGAc                                                                                                                                          >  2:250769/1‑149 (MQ=255)

            aGGAACGCGCCGAACTCGCTCTGGGCCAGCATCACCGCACGGGTGAACTGTTCGAAGTTCAGGCCCAGACGGGCCTCGACCAGTTGCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGCAGAAGCTTTGACg                                                                                                                                         <  4:262957/149‑1 (MQ=255)

               aaCGCGCCGAACTCGCTCTGGGCCAGCATCACCGCACGGGTGAACTGTTCGAAGTTCAGGCCCAGACGGGCTTCGACCAGTTGCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGCAGAAGCTTTGACGGCt                                                                                                                                      <  7:114679/149‑1 (MQ=255)

                  gcgcCGAACTCGCTCTGGGCCAGCATCACCGCACGGGTGAACTGTTCGAAGTTCAGGCCCAGACGGGCTTCGACCAGTTGCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGTAGAAGCTTTGACGGCtgtg                                                                                                                                   <  7:33510/149‑1 (MQ=255)

                   cgcCGAACTCGCTCTGGGCCAGCATCACCGCACGGGTGAACTGTTCGAAGTTCAGGCCCAGACGGGCTTCGACCAGTTGCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGCAGAAGCTTTGACGGCTGtgc                                                                                                                                  <  6:293551/149‑1 (MQ=255)

                    gcCGAACTCGCTCTGGGCCAGCATCACCGCACGGGTGAACTGTTCGAAGTTCAGGCCCAGACGGGCTTCGACCAGTTGCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGCAGAAGCTTTGACGGCTGtgct                                                                                                                                 >  8:19031/1‑149 (MQ=255)

                       gAACTCGCTCTGGGCCAGCATCACCGCACGGGTGAACTGTTCGAAGTTCAGGCCCAGACGGGCTTCGACCAGTTGCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGTAGAAGCTTTGACGGCTGtgctgc                                                                                                                               <  3:59232/148‑1 (MQ=255)

                       gAACTCGCTCTGGGCCAGCATCACCGCACGGGTGAACTGTTCGAAGTTCAGGCCCAGACGGGCTTCGACCAGTTGCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGTAGAAGCTTTGACGGCTGtgctgc                                                                                                                               >  4:59232/1‑148 (MQ=255)

                           tCGCTCTGGGCCAGCATCACCGCACGGGTGAACTGTTCGAAGTTCAGGCCCAGACGGGCTTCGACCAGTTGCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGCAGAAGCTTTGACGGCTGTGCTGCAGttt                                                                                                                          >  4:229818/1‑149 (MQ=255)

                             gCTCTGGGCCAGCATCACCGCACGGGTGAACTGTTCGAAGTTCAGGCCCAGACGGGCTTCGACCAGTTGCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGCAGAAGCTTTGACGGCTGTGCTGCAGttt                                                                                                                          <  8:287786/147‑1 (MQ=255)

                              ctctGGGCCAGCATCACCGCACGGGTGAACTGTTCGAAGTTCAGGCCCAGACGGGCTTCGACCAGTTGCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGATGTCGAGGTCGTAGAAGCTTTGACGGCTGTGCTGCAGTTTGca                                                                                                                       <  4:145786/149‑2 (MQ=255)

                               tctGGGCCAGCATCACCGCACGGGTGAACTGTTCGAAGTTCAGGCCCAGACGGGCTTCGACCAGTTGCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGCAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCg                                                                                                                      >  1:219650/1‑149 (MQ=255)

                                   ggCCAGCATCACCGCACGGGTGAACTGTTCGAAGTTCAGGCCCAGACGGGCTTCGACCAGTTGCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGTAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGt                                                                                                                     <  1:266237/146‑1 (MQ=255)

                                             agcGCACGGGTGAACTGTTCGAAGTTCAGGCCTAGACGGGCTTCGACCAGTTGCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGTAGaa                                                                                                                                                  <  4:191984/105‑1 (MQ=255)

                                             aCCGCACGGGTGAACTGTTCGAAGTTCAGGCCCAGACGGGCTTCGACCAGTTGCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGTAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTcgcg                                                                                                        >  4:48583/1‑149 (MQ=255)

                                             aCCGCACGGGTGAACTGTTCGAAGTTCAGGCCAAGACGGGCTTCGACCAGTTGCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGTAGaa                                                                                                                                                  >  3:191984/1‑107 (MQ=255)

                                                          aCTGTTCGAAGTTCAGGCCCAGACGGGCTTCGACCAGTTGCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGTAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCt                                                                                           >  5:327506/1‑149 (MQ=255)

                                                              ttCGAAGTTCAGGCCCAGACGGGCTTCGACCAGTTGCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGCAGAAGCTTTGACGGCtgt                                                                                                                                    >  2:210876/1‑104 (MQ=255)

                                                              ttCGAAGTTCAGGCCCAGACGGGCTTCGACCAGTTGCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGCAGAAGCTTTGACGGCtgt                                                                                                                                    <  1:210876/104‑1 (MQ=255)

                                                                 gAAGTTCAGGCCCAGACGGGCTTCGACCAGTTGCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGCAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGc                                                                                    >  6:75003/1‑149 (MQ=255)

                                                                 gAAGTTCAGGCCCAGACGGGCTTCGACCAGTTGCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGCAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGc                                                                                    >  7:85072/1‑149 (MQ=255)

                                                                 gAAGTTCAGGCCCAGACGGGCTTCGACCAGTTGCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGCAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGATGGCTTCCCAGc                                                                                    >  8:209467/1‑149 (MQ=255)

                                                                   aGTTCAGGCCCAGACGGGCTTCGACCAGTTGCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGCAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCgg                                                                                  <  4:74449/149‑1 (MQ=255)

                                                                    gTTCAGGCCCAGACGGGCTTCGACCAGTTGCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGTAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCggg                                                                                 >  2:37550/1‑149 (MQ=255)

                                                                          gCCCAGACGGGCTTCGACCAGTTGCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGCAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGtt                                                                                               >  7:17140/1‑129 (MQ=255)

                                                                          gCCCAGACGGGCTTCGACCAGTTGCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGCAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGtt                                                                                               <  8:17140/129‑1 (MQ=255)

                                                                           cccAGACGGGCTTCGACCAGTTGCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGTAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTa                                                                          <  3:48583/149‑1 (MQ=255)

                                                                            ccAGACGGGCTTCGACCAGTTGCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGTAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAg                                                                         >  8:295945/1‑149 (MQ=255)

                                                                               gACGGGCTTCGACCAGTTGCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGCAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAgcgg                                                                      >  2:311671/1‑149 (MQ=255)

                                                                                     cTTCGACCAGTTGCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGTAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCg                                                                >  6:26676/1‑149 (MQ=255)

                                                                                     cTTCGACCAGTTGCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGCAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCg                                                                >  4:185136/1‑149 (MQ=255)

                                                                                      ttCGACCAGTTGCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGCAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCgg                                                                                  <  8:80832/130‑1 (MQ=255)

                                                                                      ttCGACCAGTTGCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGCAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCgg                                                                                  >  7:80832/1‑130 (MQ=255)

                                                                                        cGAACAGTTGCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGTAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTc                                                              <  3:182079/148‑1 (MQ=255)

                                                                                          aCCAGTTGCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGCAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGAt                                                           <  3:282177/149‑1 (MQ=255)

                                                                                          aCCAGTTGCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGCAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGAt                                                           <  1:311671/149‑1 (MQ=255)

                                                                                          aCCAGTTGCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGCAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGAt                                                           <  1:4476/149‑1 (MQ=255)

                                                                                              gTTGCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGTAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGcc                                                        >  8:160440/1‑148 (MQ=255)

                                                                                              gTTGCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGTAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGcc                                                        >  2:249594/1‑148 (MQ=255)

                                                                                              gTTGCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGCAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGcc                                                        >  1:254362/1‑148 (MQ=255)

                                                                                               ttGCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGTAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCcga                                                      >  4:187386/1‑149 (MQ=255)

                                                                                                 gCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGCAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCcgac                                                     <  5:75003/148‑1 (MQ=255)

                                                                                                  cTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGCAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCcgacg                                                    >  5:272623/1‑148 (MQ=255)

                                                                                                   ttGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGTAGAAGCTTTGACGGCTGTGCTGCAGt                                                                                                                            >  5:175352/1‑75 (MQ=255)

                                                                                                   ttGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGTAGAAGCTTTGACGGCTGTGCTGCAGt                                                                                                                            <  6:175352/75‑1 (MQ=255)

                                                                                                   ttGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGTAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCcgacga                                                   >  7:149448/1‑148 (MQ=255)

                                                                                                    tGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGCAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCTGACGaaa                                                 >  4:163472/1‑149 (MQ=255)

                                                                                                     gTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGTAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCCGACGAAAt                                                >  2:183385/1‑149 (MQ=255)

                                                                                                      tATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGTAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCCGACGAAATC                                               <  3:143407‑M1/149‑1 (MQ=255)

                                                                                                       aTTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGCAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCCGACGAAATC                                               >  8:294968‑M1/1‑148 (MQ=255)

                                                                                                        ttCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGTAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCCGACGAAATCGA                                             >  4:270306‑M1/1‑149 (MQ=255)

                                                                                                        ttCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGTAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCCGACGAAATCGA                                             >  1:130731‑M1/1‑149 (MQ=255)

                                                                                                        ttCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGTAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCCGACGAAATCGA                                             >  5:17394‑M1/1‑149 (MQ=255)

                                                                                                        ttCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGCAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCCGACGAAATCGA                                             >  1:160000‑M1/1‑149 (MQ=255)

                                                                                                        ttCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGCAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCCGACGAAATCGA                                             <  1:109752‑M1/149‑1 (MQ=255)

                                                                                                        ttCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGACGAGGTCGTAGAAGCTATGACGGCTGTGCTGCAGTTTGCCGTTGGCCATGTCGCGGGCGCGGTTCGCTTCCCAGCGGGCGCGGTAGCGGCGGcc                                                                 >  1:183250/1‑129 (MQ=255)

                                                                                                        ttCGTTCTTGCCACTGCCGAGCACTTGCGCGCTGTCGAGGTCGTAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCCGACGAAATCGA                                             <  7:295945‑M1/149‑1 (MQ=255)

                                                                                                           gTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGTAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCCGACGAAATCGACCt                                          >  1:310478/1‑149 (MQ=255)

                                                                                                           gTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGTAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCCGACGAAATCGACC                                           <  1:249594‑M1/148‑1 (MQ=255)

                                                                                                            ttCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGTAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTg                                                                                              <  7:63418/96‑1 (MQ=255)

                                                                                                            ttCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGTAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTg                                                                                              >  8:63418/1‑96 (MQ=255)

                                                                                                            ttCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGCAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTAGATGCCGACGAAATCGACCTc                                         >  6:285594/1‑149 (MQ=255)

                                                                                                              cTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGCAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCTCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCCGACGAAATCGACCTCg                                        <  1:168134/148‑1 (MQ=255)

                                                                                                               ttGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGTAGAGGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCCGACGAAATCGACCTCgg                                       >  6:36334/1‑148 (MQ=255)

                                                                                                               ttGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGCAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCCGACGAAATCGACCTCgg                                       >  3:12413/1‑148 (MQ=255)

                                                                                                                 gCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGTAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCCGACGAAATCGACCTCGGc                                      <  1:37550/147‑1 (MQ=255)

                                                                                                                 gCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGCAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCCGACGAAATCGACCTCgg                                       <  7:209467/146‑1 (MQ=255)

                                                                                                                  ccACTGCCGAGCACTTGCTCGCTGTCGAGGTCGTAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGt                                                                           <  2:48086/109‑1 (MQ=255)

                                                                                                                   cACTGCCGAGCACTTGCTCGCTGTCGAGGTCGTAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCCGACGAAATCGACCTCGGCaaaa                                  <  2:310478/149‑1 (MQ=255)

                                                                                                                   cACTGCCGAGCACTTGCTCGCTGTCGAGGTCGTAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCCGACGAAATCGACCTCGGCaaaa                                  >  2:145497/1‑149 (MQ=255)

                                                                                                                      tGCCGAGCACTTGCTCGCTGTCGAGGTCGTAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCCGACGAAATCGACCTCGGCAAAACCg                               >  5:91640/1‑149 (MQ=255)

                                                                                                                        ccGAGCACTTGCTCGCTGTCGAGGTCGTAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCCGACGAAATCGACCTCGGCAAAACCGCt                             >  5:219412/1‑149 (MQ=255)

                                                                                                                           aGCACTTGCTCGCTGTCGAGGTCGCAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCCGACGAAATCGACCTCGGCAAAACCGCTGcc                          >  6:54323/1‑149 (MQ=255)

                                                                                                                           aGCACTTGCTCGCTGTCGAGGTCGCAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCCGACGAAATCGACCTCGGCAAAACCGCTGcc                          >  6:49473/1‑149 (MQ=255)

                                                                                                                             cACTTGCTCGCTGTCGAGGTCGTAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCCGACGAAATCGACCTCGGCAAAACCGCTGCCgg                        >  4:505/1‑149 (MQ=255)

                                                                                                                              aCTTGCTCGCTGTCGAGGTCGTAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCCGACGAAATCGACCTCGGCAAAACCGCTGCCGGt                       <  6:327506/149‑1 (MQ=255)

                                                                                                                              aCTTGCTCGCTGTCGAGGTCGTAGAAGCTTTGAAGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCCGACGAAATCGACCTCGGCAAAACCGCTGCCGGt                       <  7:287825/149‑1 (MQ=255)

                                                                                                                              aCTTGCTCGCTGTCGAGGTCGCAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCCGACGAAATCGACCTCGGCAAAACCGCTGCCGGt                       <  4:12413/149‑1 (MQ=255)

                                                                                                                               cTTGCTCGCTGTCGAGGTCGCAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCCGACGAAATCGACCTCGGCAAAACCGCTGCCGGTg                      <  5:49473/149‑1 (MQ=255)

                                                                                                                               cTTGCTCGCTGTCGAGGTCGCAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCCGACGAAATCGACCTCGGCAAAACCGCTGCCGGTg                      <  5:54323/149‑1 (MQ=255)

                                                                                                                                  gCTCGCTGTCGAGGTCGCAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCTGACGAAATCGACCTCGGCAAAACCGCTGCCGGTGc                     <  3:163472/147‑1 (MQ=255)

                                                                                                                                   cTCGCTGTCGAGGTCGTAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCCGACGAAATCGACCTCGGCAAAACCGCTGCCGGTGCCAc                  >  4:169289/1‑149 (MQ=255)

                                                                                                                                   cTCGCTGTCGAGGTCGCAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCCGACGAAATCGACCTCGGCAAAACCGCTGCCGGTGCCAc                  >  1:264808/1‑149 (MQ=255)

                                                                                                                                      gCTGTCGAGGTCGTAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCCGACGAAATCGACCTCGGCAAAACCGCTGCCGGTGCCACGgc               >  7:7216/1‑149 (MQ=255)

                                                                                                                                          tCGAGGTCGTAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCCGACGAAATCGACCTCGGCAAAACCGCTGCCGGTGCCACGGTgcag           >  7:67818/1‑149 (MQ=255)

                                                                                                                                           cGAGGTCGCAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCCGACGAAATCGACCTCGGCAAAACCGCTGCCGGTGCCACGGCgcagc          <  3:229818/149‑1 (MQ=255)

                                                                                                                                            gAGGTCGCAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCCGACGAAATCGACCTCGGCAAAACCGCTGCCGGTGCCACGGCgcag           >  2:199923/1‑147 (MQ=255)

                                                                                                                                            gAGGTCGCAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCCGACGAAATCGACCTCGGCAAAACCGCTGCCGGTGCCACGGCgcag           >  1:142684/1‑147 (MQ=255)

                                                                                                                                             aGGTCGCAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCCGACGAAATCGACCTCGGCAAAACCGCTGCCGGTGCCACGGCGcagcag        <  1:53410/149‑1 (MQ=255)

                                                                                                                                             aGGTCGCAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCCGACGAAATCGACCTCGGCAAAACCGCTGCCGGTGCCACGGCGcagcag        >  3:187006/1‑149 (MQ=255)

                                                                                                                                                tCGTAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCCGACGAATACGACCTCGGCAAAACCGCTGCCGTTGCCACGGCGCAGCAGGtt     >  7:127930/1‑149 (MQ=255)

                                                                                                                                                tCGTAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCCGACGAAATCGACCTCGGCAAAACCGCTGCCGGTGCCACGGCGCAGCAGGtt     >  8:56460/1‑149 (MQ=255)

                                                                                                                                                 cGTAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCCGACGAAATCGACCTCGGCAAAACCGCTGCCGGTGCCACGGCGCAGCAGGTTg    >  5:110404/1‑149 (MQ=255)

                                                                                                                                                   tAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCCGACGAAATCGACCTCGGCAAAACCGCTGCCGGTGCCACGGCGCAGCAGGTTgcg  >  3:175827/1‑149 (MQ=255)

TCGTCGGCCTTGAGGAACGCGCCGAACTCGCTCTGGGCCAGCATCACCGCACGGGTGAACTGTTCGAAGTTCAGGCCCAGACGGGCTTCGACCAGTTGCTTGTATTCGTTCTTGCCACTGCCGAGCACTTGCTCGCTGTCGAGGTCGTAGAAGCTTTGACGGCTGTGCTGCAGTTTGCCGTTGGCCTTGTCGCGGGCGCGGTTGGCTTCCCAGCGGGCGCGGTAGCGGCGGCCGTCGATGCCGACGAAATCGACCTCGGCAAAACCGCTGCCGGTGCCACGGCGCAGCAGGTTGCG  >  NC_002947/2299580‑2299875

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