BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,841,649	T→C	6.4%	S71G (AGC→GGC) ‡	PP_5529 ←	epimerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,841,649	0	T	C	6.4%	582.5 / 18.4	186	S71G (AGC→GGC) ‡	PP_5529	epimerase
Reads supporting (aligned to +/- strand): ref base T (107/67); new base C (3/9); total (110/76)
Fisher's exact test for biased strand distribution p-value = 1.63e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

Only 100 of 188 total aligned reads displayed.

CAGTGCAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTTCTCGTGCGGCAGGCGCCTCCTGCGCGT  >  NC_002947/2841500‑2841798

  gttcAATGACATGCACAAGCCGAGATAAACGATCATAAAGACTTGAGCACGCGTTACATATAGGAGATTAGACGGAGCCTTCGATTAGAAAAACTAGGTAATACCCAACGTAGAGACCATGATATATCGCAGCATCGTAAGCCTGGCTgc                                                                                                                                                     >  2:2624581/4‑149 (MQ=255)

     cAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCAt                                                                                                                                                  >  1:3180504/1‑150 (MQ=255)

       aTGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTATGCCAGCCCGGCTTAc                                                                                                                                                <  1:5191481/150‑1 (MQ=255)

        tGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATAcc                                                                                                                                               >  1:936538/1‑150 (MQ=255)

         gACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCa                                                                                                                                              <  2:4293148/150‑1 (MQ=255)

          aCATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAg                                                                                                                                             >  2:975732/1‑150 (MQ=255)

            aTGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGc                                                                                                                                           >  2:2173832/1‑150 (MQ=255)

             tGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTATGCCAGCCCGGCTTACCAGGcg                                                                                                                                          <  2:373363/150‑2 (MQ=255)

                acaAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGa                                                                                                                                       <  1:1037684/150‑1 (MQ=255)

                acaAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGa                                                                                                                                       <  1:4465157/150‑1 (MQ=255)

                acaAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGa                                                                                                                                       >  1:2404946/1‑150 (MQ=255)

                    gCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCgg                                                                                                                                   <  2:5288490/150‑1 (MQ=255)

                     ccGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGc                                                                                                                                  >  1:4812968/1‑150 (MQ=255)

                      cGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCt                                                                                                                                 >  1:916768/1‑150 (MQ=255)

                      cGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCt                                                                                                                                 >  1:2117417/1‑150 (MQ=255)

                      cGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCt                                                                                                                                 >  1:2484112/1‑150 (MQ=255)

                          atatACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCa                                                                                                                             >  2:4802475/1‑150 (MQ=255)

                           tataCGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCat                                                                                                                            >  1:4390714/1‑150 (MQ=255)

                                 cTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTc                                                                                                                      >  1:1705217/1‑150 (MQ=255)

                                 cTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTc                                                                                                                      <  2:2448720/150‑1 (MQ=255)

                                 cTCATAACGACGTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTc                                                                                                                      <  2:5393955/150‑1 (MQ=255)

                                  tCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATTTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCAGCACCGGAGAGACGATGCTTTAGCGCCGACTGGTAAGCCGGGATGGCATACAAGGCTCGAGCGGCTGTCATATCTTCg                                                                                                                     >  1:1060299/1‑150 (MQ=255)

                                    aTAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGaa                                                                                                                   >  1:1549244/1‑150 (MQ=255)

                                         gACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTa                                                                                                              <  2:2127737/150‑1 (MQ=255)

                                            ttGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGaa                                                                                                           >  2:2607478/1‑150 (MQ=255)

                                            ttGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGaa                                                                                                           <  1:3853129/150‑1 (MQ=255)

                                               tgCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGa                                                                                                        <  1:2607478/150‑1 (MQ=255)

                                               tgCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGa                                                                                                        >  1:996288/1‑150 (MQ=255)

                                                gCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGAc                                                                                                       >  1:1009548/1‑150 (MQ=255)

                                                 cACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGAcc                                                                                                      >  1:3195287/1‑150 (MQ=255)

                                                 cACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGAcc                                                                                                      >  2:4450338/1‑150 (MQ=255)

                                                           aGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCtcga                                                                                            <  2:712303/150‑1 (MQ=255)

                                                            gTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCTGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCtcgat                                                                                           <  1:21001/150‑1 (MQ=255)

                                                            gTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCtcgat                                                                                           >  1:859988/1‑150 (MQ=255)

                                                            gTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCtcgat                                                                                           >  2:2501873/1‑150 (MQ=255)

                                                                ggTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTcgatcgat                                                                                       >  1:4062085/1‑150 (MQ=255)

                                                                ggTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTcgatcgat                                                                                       <  1:3269481/150‑1 (MQ=255)

                                                                ggTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTcgatcgat                                                                                       <  1:4585821/150‑1 (MQ=255)

                                                                   gATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGa                                                                                    <  1:3768105/150‑1 (MQ=255)

                                                                     tgtgTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGAcc                                                                                  >  1:71964/1‑150 (MQ=255)

                                                                         tCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTcc                                                                              >  1:4416666/1‑150 (MQ=255)

                                                                            gTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAac                                                                           >  1:3026465/1‑150 (MQ=255)

                                                                            gTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAac                                                                           >  1:5271994/1‑150 (MQ=255)

                                                                              gCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAacac                                                                         <  1:3480092/150‑1 (MQ=255)

                                                                                cTTCGATTAGGAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCt                                                                       >  2:2069651/1‑150 (MQ=255)

                                                                                 ttCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCTATCGATGGACCTTCCAACACCTc                                                                      >  2:4009044/1‑150 (MQ=255)

                                                                                 ttCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTc                                                                      >  1:3835498/1‑150 (MQ=255)

                                                                                 ttCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTc                                                                      >  2:3360042/1‑150 (MQ=255)

                                                                                 ttCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTc                                                                      <  1:1143597/150‑1 (MQ=255)

                                                                                    gATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAg                                                                   >  2:1847485/1‑150 (MQ=255)

                                                                                    gATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAg                                                                   >  2:3076155/1‑150 (MQ=255)

                                                                                    gATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAg                                                                   <  1:3301919/150‑1 (MQ=255)

                                                                                    gATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAg                                                                   >  1:2212958/1‑150 (MQ=255)

                                                                                    gATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAg                                                                   <  2:5065133/150‑1 (MQ=255)

                                                                                    caatAGAAAAACTCCGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCCATCGATGGACCTTCCAACACCTCCAg                                                                   <  2:5065040/147‑1 (MQ=255)

                                                                                      ttAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGtt                                                                 <  2:3292595/150‑1 (MQ=255)

                                                                                         gAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGAc                                                              >  1:3619641/1‑150 (MQ=255)

                                                                                         gAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGAc                                                              >  2:1053758/1‑150 (MQ=255)

                                                                                         gAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGAc                                                              <  2:5154731/150‑1 (MQ=255)

                                                                                          aaaaaCTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGTCTGGCATAACAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATTACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGaac                                                             >  1:5363330/1‑148 (MQ=255)

                                                                                          aaaaaCTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGAcc                                                             >  1:3185967/1‑150 (MQ=255)

                                                                                            aaaCTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAAAGACCGCTCCCGCGATCGATGGACCTTACAACACCTCCAGTTGACCGt                                                           >  2:2074903/1‑150 (MQ=255)

                                                                                             aaCTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTa                                                          <  1:620605/150‑1 (MQ=255)

                                                                                             aaCTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTa                                                          <  2:136392/150‑1 (MQ=255)

                                                                                             aaCTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTa                                                          <  1:3058760/150‑1 (MQ=255)

                                                                                               cTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGa                                                        <  2:1548451/150‑1 (MQ=255)

                                                                                               cTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGa                                                        >  2:3799769/1‑150 (MQ=255)

                                                                                                tCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGaa                                                       >  2:1224062/1‑150 (MQ=255)

                                                                                                tCGGAACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGCAGGACCTTCCAACACCTCCAGTTGACCGTAGaa                                                       >  1:5221631/1‑150 (MQ=255)

                                                                                                   gTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGc                                                    >  1:4325351/1‑150 (MQ=255)

                                                                                                   gTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGc                                                    <  1:4231221/150‑1 (MQ=255)

                                                                                                      cTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAg                                                 >  1:2061789/1‑150 (MQ=255)

                                                                                                       tCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCa                                                                                                                                              <  1:4095225/56‑1 (MQ=255)

                                                                                                       tCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGa                                                >  1:1288371/1‑150 (MQ=255)

                                                                                                       tCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGa                                                >  2:375739/1‑150 (MQ=255)

                                                                                                        ccGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGAc                                               <  1:4675568/150‑1 (MQ=255)

                                                                                                         cGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGGTGACCGTAGAAAGCCAGACg                                              >  2:2931588/1‑150 (MQ=255)

                                                                                                          gCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACgt                                             <  2:4416330/150‑1 (MQ=255)

                                                                                                          gCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCAGATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACgt                                             >  1:1588229/1‑150 (MQ=255)

                                                                                                            aCCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACgtgt                                           <  2:2910048/150‑1 (MQ=255)

                                                                                                            aCCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACgtgt                                           <  1:3360042/150‑1 (MQ=255)

                                                                                                                    gaCGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGAcgtcgt                                   <  1:1155770/150‑1 (MQ=255)

                                                                                                                     aCGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGACACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTg                                  >  2:2393391/1‑150 (MQ=255)

                                                                                                                           cTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTtc                            <  2:2700705/150‑1 (MQ=255)

                                                                                                                            tttAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTtct                           >  2:1989219/1‑150 (MQ=255)

                                                                                                                            tttAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGAGGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTtct                           >  2:1526012/1‑150 (MQ=255)

                                                                                                                                gcgcCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTTCTCGTg                       >  2:3868322/1‑150 (MQ=255)

                                                                                                                                gcgcCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTTCTCGTg                       >  2:1128362/1‑150 (MQ=255)

                                                                                                                                    agccTGGTATGCCAGCCCGGCTTACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTTCTCGTGCGGc                   >  1:1566940/2‑150 (MQ=255)

                                                                                                                                      ccTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTTCTCGTGCGGCAg                 >  2:2785802/1‑150 (MQ=255)

                                                                                                                                       cTGGTATGCCAGCCCGGCTTACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTTCTCGTGCGGCAgg                >  2:1663423/1‑150 (MQ=255)

                                                                                                                                        tGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTTCTCGTGCGGCAGgc               <  2:916768/150‑1 (MQ=255)

                                                                                                                                          tgatcacGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTTCTCGTGCGGCAGgcgc             >  2:5163902/7‑150 (MQ=255)

                                                                                                                                          gTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTTCTCGTGCGGCAGgcgc             >  1:4235208/1‑150 (MQ=255)

                                                                                                                                               ccGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTTCTCGTGCGGCAGGCGcctcct        >  1:1410721/1‑150 (MQ=255)

                                                                                                                                               ccGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTTCTCGTGCGGCAGGCGcctcct        >  2:3040743/1‑150 (MQ=255)

                                                                                                                                                   gCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTTCTCGTGCGGCAGGCGCCTCCTgcgc    <  2:4166833/150‑1 (MQ=255)

                                                                                                                                                   gCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTTCTCGTGCGGCAGGCGCCTCCTgcgc    <  1:3524202/150‑1 (MQ=255)

                                                                                                                                                   gCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTTCTCGTGCGGCAGGCGCCTCCTgcgc    >  1:3684956/1‑150 (MQ=255)

                                                                                                                                                     tGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTTCTCGTGCGGCAGGCGCCTCCTGCGCGt  >  1:4991116/1‑150 (MQ=255)

CAGTGCAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTTCTCGTGCGGCAGGCGCCTCCTGCGCGT  >  NC_002947/2841500‑2841798

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