BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,403,831	C→T	5.5%	C69Y (TGC→TAC) ‡	PP_3883 ←	holin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,403,831	0	C	T	5.5%	284.2 / 3.0	91	C69Y (TGC→TAC) ‡	PP_3883	holin
Reads supporting (aligned to +/- strand): ref base C (44/42); new base T (3/2); total (47/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.94e-01

TTTCATGGGGCCCCGCAGCGCCCGGAGGCTTTCAGGAGGGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGCGGGTTCTTGCGCAACCATCCTTCACGAGC  >  NC_002947/4403685‑4403979

tttCATGGGGCCCCGCAGCGCCCGGAGGCTTTCAGGAGGGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGCCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGca                                                                                                                                                      >  8:240892/1‑148 (MQ=255)

 ttCATGGGGCCCCGCAGCGCCCGGAGGCTTTCAGGAGGGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGca                                                                                                                                                      >  7:184834/1‑147 (MQ=255)

  tCATGGGGCCCCGCAGCGCCCGGAGGCTTTCAGGAGGGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGcaca                                                                                                                                                    <  3:325224/148‑1 (MQ=255)

   cATGGGGCCCCGCAGCGCCCGGAGGCTTTCAGGAGGGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAg                                                                                                                                                   <  1:64547/148‑1 (MQ=255)

    aTGGGGCCCCGCAGCGCCCGGAGGCTTTCAGGAGGGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCg                                                                                                                                                 <  6:181756/149‑1 (MQ=255)

       gggCCCCGCAGCGCCCGGAGGCTTTCAGGAGGGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGAc                                                                                                                                               <  1:292785/148‑1 (MQ=255)

        ggCCCCGCAGCGCCCGGAGGCTTTCAGGAGGGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTg                                                                                                                                             >  4:165489/1‑149 (MQ=255)

         gCCCCGCAGCGCCCGGAGGCTTTCAGGAGGGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGc                                                                                                                                            <  6:294689/149‑1 (MQ=255)

          ccccGCAGCGCCCGGAGGCTTTCAGGAGGGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTCGCTGTGCACCAACTGcc                                                                                                                                           <  4:188593/149‑1 (MQ=255)

           cccGCAGCGCCCGGAGGCTTTCAGGAGGGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTCGCTGTGCACCAACTGCCa                                                                                                                                          <  4:226186/149‑1 (MQ=255)

            ccGCAGCGCCCGGAGGCTTTCAGGAGGGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAg                                                                                                                                         >  5:74950/1‑149 (MQ=255)

                     ccGGAGGCTTTCAGGAGGGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTcc                                                                                                                                  >  3:228005/1‑147 (MQ=255)

                        gAGGCTTTCAGGAGGGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGcc                                                                                                                             >  2:335319/1‑149 (MQ=255)

                          ggCTTTCAGGAGGGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCa                                                                                                                           <  4:183670/149‑1 (MQ=255)

                               tCAGGAGGGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCa                                                                                                                           <  2:351893/144‑1 (MQ=255)

                               tCAGGAGGGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCa                                                                                                                           >  1:351893/1‑144 (MQ=255)

                                cAGGAGGGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGaa                                                                                                                        >  2:71688/1‑146 (MQ=255)

                                 aggaggGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACaga                                                                                                                    <  4:372028/149‑1 (MQ=255)

                                   gaggGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACagag                                                                                                                   <  4:408091/148‑1 (MQ=255)

                                   cggggCCTATTCCAGTGGTGCCGGTGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACtga                                                                                                                    <  2:79879/145‑1 (MQ=255)

                                        ccTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGt                                                                                                             >  6:351364/1‑149 (MQ=255)

                                         cTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGTGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTc                                                                                                            >  5:227306/1‑149 (MQ=255)

                                         cTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGACCAGAGTGATGTc                                                                                                            >  2:424329/1‑149 (MQ=255)

                                          tATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCAGACTTCTCCAGCCGCTGCTTGAAGGAGTTGTGGCAGTAGGTGCACAGCGACAGCCACTTGGTCCGAGCCCAGAACAGAGTCATGTcg                                                                                                           >  4:330004/1‑149 (MQ=255)

                                              ccAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTcg‑ccg                                                                                                       <  1:335319/148‑1 (MQ=255)

                                               cAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTCGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAAc                                                                                                                       >  5:73666/1‑132 (MQ=255)

                                               cAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTCGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAAc                                                                                                                       <  6:73666/132‑1 (MQ=255)

                                                aGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCgg                                                                                                    >  5:429439/1‑149 (MQ=255)

                                                 gtggtgCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGt                                                                                                   >  2:306556/1‑149 (MQ=255)

                                                 gtggtgCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGt                                                                                                   >  8:21912/1‑149 (MQ=255)

                                                  tggtgCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTg                                                                                                  <  5:351364/149‑1 (MQ=255)

                                                    gtgCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCg                                                                                                >  4:323369/1‑149 (MQ=255)

                                                    gtgCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCg                                                                                                <  6:368851/149‑1 (MQ=255)

                                                    gtgCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCg                                                                                                >  1:362020/1‑149 (MQ=255)

                                                       ccGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGc                                                                                              <  8:324520/148‑1 (MQ=255)

                                                       ccGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCt                                                                                             <  3:165489/149‑1 (MQ=255)

                                                         gggggTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTg                                                                                            >  8:91443/1‑148 (MQ=255)

                                                         gggggTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGa                                                                                           <  6:272818/149‑1 (MQ=255)

                                                           gggTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAAt                                                                                         <  3:323369/149‑1 (MQ=255)

                                                             gTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAAtgt                                                                                       >  7:381669/1‑149 (MQ=255)

                                                               ttaGCGGCAGGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGa                                                                                     >  7:376544/3‑149 (MQ=255)

                                                                cAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGAt                                                                                    >  3:132747/1‑149 (MQ=255)

                                                                 aGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATc                                                                                   >  3:361617/1‑149 (MQ=255)

                                                                 aGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATc                                                                                   >  6:212824/1‑149 (MQ=255)

                                                                     gCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATcgac                                                                                >  1:91793/1‑148 (MQ=255)

                                                                         gCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAg                                                                           >  3:214931/1‑149 (MQ=255)

                                                                          ccGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACa                                                                            <  4:132747/147‑1 (MQ=255)

                                                                            gCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGa                                                                                     >  2:93195/1‑136 (MQ=255)

                                                                            gCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGa                                                                                     <  1:93195/136‑1 (MQ=255)

                                                                                aCGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCGTCCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCg                                                                      <  1:71688/148‑1 (MQ=255)

                                                                                 cGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGc                                                                   <  2:362020/149‑1 (MQ=255)

                                                                                  gTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCt                                                                  >  2:116141/1‑149 (MQ=255)

                                                                                  gTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCt                                                                  >  5:329029/1‑149 (MQ=255)

                                                                                   tCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGCATGGGAGCGACGCCAGTGGCGGGCtt                                                                 >  6:333481/1‑149 (MQ=255)

                                                                                    cGCAACCAGCCTCGCGCCCGGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGCGATCGACGACAGTGGCGGGCTTg                                                                <  8:184834/149‑1 (MQ=255)

                                                                                     gCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGc                                                               >  5:33829/1‑149 (MQ=255)

                                                                                      cAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGTGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGATAGTGGCGGGCTTGCt                                                              <  6:227306/149‑1 (MQ=255)

                                                                                      cAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCt                                                              >  1:431998/1‑149 (MQ=255)

                                                                                            gCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCtt                                                                 >  1:279678/1‑140 (MQ=255)

                                                                                            gCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCtt                                                                 <  2:279678/140‑1 (MQ=255)

                                                                                             ccTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCt                                                       <  5:205708/149‑1 (MQ=255)

                                                                                                 gcgcCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACt                                                   <  7:91443/149‑1 (MQ=255)

                                                                                                    cccTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGa                                                <  4:361617/149‑1 (MQ=255)

                                                                                                     ccTGACTTCTCCAGGCGCTGCTTGTCGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGGCCGGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAg                                               <  6:74950/149‑1 (MQ=255)

                                                                                                     ccTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAg                                               <  7:21912/149‑1 (MQ=255)

                                                                                                     ccTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAg                                               <  2:431998/149‑1 (MQ=255)

                                                                                                      cTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAgg                                              <  7:240892/149‑1 (MQ=255)

                                                                                                         aCTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGc                                             <  1:116141/147‑1 (MQ=255)

                                                                                                           ttCTCCAGGCGCTGCTTCTAGGAGTTGTGGCAGTTGGAGCACATCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCT‑CCGCTGTGCGGCTGAATGTGAACGACGACAGTGGCGGGCATCCTCAGGCCTCACAGGAGGCAGc                                          >  8:398608/1‑148 (MQ=255)

                                                                                                             ctcCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGaa                                       >  6:214732/1‑149 (MQ=255)

                                                                                                              tcCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGaa                                       >  4:228860/1‑148 (MQ=255)

                                                                                                                gaGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACTGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCTGGCTTGCTCAGGCCTGACTGGAGGCAGCGAaca                                     <  7:398608/147‑1 (MQ=255)

                                                                                                                cAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAaca                                     >  3:108827/1‑148 (MQ=255)

                                                                                                                     gctgctTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGc                               <  2:91793/149‑1 (MQ=255)

                                                                                                                      ctgctTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGCg                              >  5:202032/1‑149 (MQ=255)

                                                                                                                         ctTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGCgg                             <  3:228860/147‑1 (MQ=255)

                                                                                                                         ctTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGCGGGt                           <  5:212824/149‑1 (MQ=255)

                                                                                                                          ttGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGCGGGtt                          >  6:224689/1‑149 (MQ=255)

                                                                                                                                 ccaactgGCAGTCGCTGTGCACCAACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGCGAGGGTTCTTgcg                    >  7:51417/6‑148 (MQ=255)

                                                                                                                                 aGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGCGGGTTCTTgcgc                   <  7:64599/149‑1 (MQ=255)

                                                                                                                                 aGTTGTGGCAGTCGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGc                                                         <  1:16598/111‑1 (MQ=255)

                                                                                                                                 aGTTGTGGCAGTCGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGc                                                         >  2:16598/1‑111 (MQ=255)

                                                                                                                                  gTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGCGGGTTCTTgcgc                   >  8:363955/1‑148 (MQ=255)

                                                                                                                                    actgGCAGTCGCTGTGCACCAACTGCCAGTTGGTCTGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGTCGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGCACACAGAGCGGGTTCTTGCGCAAc                >  4:211219/3‑149 (MQ=255)

                                                                                                                                    actgGCAGTCGCTGTGCACCAACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGCGGGTTCTTGCGCAAc                >  6:126426/3‑149 (MQ=255)

                                                                                                                                     gtgGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGCGGGTTCTTGCGCAAc                >  6:80168/1‑148 (MQ=255)

                                                                                                                                      tgGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGCGGGTTCTTGCGCAAcc               <  6:202032/148‑1 (MQ=255)

                                                                                                                                        gCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGCGGGTTCTTGCGCAACCATc            >  5:385307/1‑149 (MQ=255)

                                                                                                                                         cAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCGTGACTGGAGGCAGCGAACACAGAGCGGGTTCTTGCGCAACCATc            <  4:371250/148‑1 (MQ=255)

                                                                                                                                         cAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGCGGGTTCTTGCGCAACCATc            >  3:371250/1‑148 (MQ=255)

                                                                                                                                                tGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGCGGGTTCTTGCGCAACCATCCTTCa       <  2:186949/146‑1 (MQ=255)

                                                                                                                                                tGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGCGGGTTCTTGCGCAACCATCCTTCa       >  1:186949/1‑146 (MQ=255)

                                                                                                                                                tGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGCGGGTTCTTGCGCAACCATCCTTCAc      >  5:185377/1‑147 (MQ=255)

                                                                                                                                                tGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGCGGGTTCTTGCGCAACCATCCTTCACGa    <  6:329029/149‑1 (MQ=255)

                                                                                                                                                  cacaGCGACTGCCAGTTGGTCCGAGTCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGCGGGTTCTTGCGCAACCATCCTTCACGAGc  >  3:10043/1‑149 (MQ=255)

                                                                                                                                                  cacaGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGCGGGTTCTTGCGCAACCATCCTTCACGAGc  <  1:306556/149‑1 (MQ=255)

                                                                                                                                                  cacaGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGCGGGTTCTTGCGCAACCATCCTTCACGAGc  >  4:106893/1‑149 (MQ=255)

TTTCATGGGGCCCCGCAGCGCCCGGAGGCTTTCAGGAGGGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGCGGGTTCTTGCGCAACCATCCTTCACGAGC  >  NC_002947/4403685‑4403979

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