BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,530,476	T→C	6.7%	R89R (CGT→CGC)	PP_4858 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,530,476	0	T	C	6.7%	671.9 / 19.9	247	R89R (CGT→CGC)	PP_4858	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (100/114); new base C (3/13); total (113/134)
Fisher's exact test for biased strand distribution p-value = 3.67e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 253 total aligned reads displayed.

CAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCGGATCGTCTGATCGAACAAGCTGGCCGTGTTGCCGAAGGCGG  >  NC_002947/5530327‑5530625

     gCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGc                                                                                                                                                  <  1:3530934/150‑1 (MQ=255)

      cTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCt                                                                                                                                                 <  1:6786909/150‑1 (MQ=255)

        gAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGCTCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGa                                                                                                                                               <  2:1094408/150‑1 (MQ=255)

         aaCACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGACCGAACGTTCTGGCCGTGTTGCTGaa                                                                                                                                              <  1:6862302/150‑1 (MQ=255)

           cacCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAgg                                                                                                                                            <  1:5499678/150‑1 (MQ=255)

            acCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAGCGATCTGGCCGCGGTGCTGAAGGAGGCTCCGAGCGCCTGATCGAACGGTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGGTCTGGCCGGGTTGcggcgggg                                                                                                                                           >  2:6207912/1‑143 (MQ=14)

            acCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGt                                                                                                                                           <  2:2889140/150‑1 (MQ=255)

            acCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGt                                                                                                                                           <  2:2839883/150‑1 (MQ=255)

              cGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTgg                                                                                                                                         <  1:5798942/150‑1 (MQ=255)

              cGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGACTTATCGAACGAGCTGGCGGCGTTGCGGAAGGTGGCTCCGGTCGACTGATCGAACGTTCTGGGCGGGTTGCTGAAGGTgg                                                                                                                                         >  2:687868/1‑150 (MQ=11)

               gTAGCTGTAGCCGGCTCCGATCCCATGATCGATCGTACTTGCAGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGc                                                                                                                                        <  1:1690968/150‑1 (MQ=11)

                tAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCt                                                                                                                                       <  2:5226305/150‑1 (MQ=40)

                 aGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTc                                                                                                                                      <  1:4355117/150‑1 (MQ=40)

                    tGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGACTGATCGAACGTTCTGGCCGGGTTGCTGAAGGTGGCTCCGa                                                                                                                                   >  1:6528168/1‑150 (MQ=21)

                    tGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGGGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGa                                                                                                                                   <  2:1613555/150‑1 (MQ=255)

                    tGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTGGCTGAAGGTGGCTCCGATCGCCTGATCGCACGGGCTGGCCGCGTTGCGGAAGGTGGCGCCGATCGCCTGATCGAACGTGCTGGCCCGGTTGCTGAAGGTGGCTCCGa                                                                                                                                   >  2:6786171/1‑150 (MQ=11)

                      aaGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATc                                                                                                                                 <  1:4876998/150‑1 (MQ=39)

                      aaGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATc                                                                                                                                 <  2:384604/150‑1 (MQ=39)

                         gcggcTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTc                                                                                                                              <  2:6386015/150‑1 (MQ=38)

                         gcggcTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTc                                                                                                                              <  1:4833768/150‑1 (MQ=38)

                             cTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGAt                                                                                                                          <  2:3674158/150‑1 (MQ=37)

                             cTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGAt                                                                                                                          <  1:2894017/150‑1 (MQ=37)

                              tCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATc                                                                                                                         <  2:3211642/150‑1 (MQ=37)

                                 gATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGaa                                                                                                                      <  1:2574730/150‑1 (MQ=37)

                                  aTCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAAc                                                                                                                     <  1:4688027/150‑1 (MQ=37)

                                    cGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTTTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGt                                                                                                                   <  1:3691380/150‑1 (MQ=37)

                                     gCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGGTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGGCTGATCGAACGGTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGACTGATCGAGCGtg                                                                                                                  >  2:4545875/1‑149 (MQ=11)

                                      ccTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTc                                                                                                                              <  2:2524165/137‑1 (MQ=37)

                                      ccTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTc                                                                                                                              >  1:2524165/1‑137 (MQ=37)

                                        tGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTg                                                                                                               <  1:98933/150‑1 (MQ=37)

                                          aTCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGc                                                                                                             <  1:2612529/150‑1 (MQ=37)

                                           tCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGGCTGATCGAAAGTTCTGGACGTGTTGCGGAATGTGGCTCCGATCGTCTGATGGAGCGGTCTGGcg                                                                                                            >  1:1496462/1‑149 (MQ=11)

                                             gAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCGGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCcgc                                                                                                          <  1:1229950/150‑1 (MQ=38)

                                                   tCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGAGGCTCCGATCGTCTGATCGAACGATCTGGCCGTGTAGCt                                                                                                    >  2:4978692/1‑150 (MQ=255)

                                                     tGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGa                                                                                                  >  1:2782511/1‑150 (MQ=255)

                                                      ggCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGaa                                                                                                 <  1:53777/150‑1 (MQ=255)

                                                      ggCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGaa                                                                                                 <  2:1486850/150‑1 (MQ=255)

                                                       gCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACTATCTGGCCGAGTTGCTGAAg                                                                                                >  1:543842/1‑150 (MQ=255)

                                                        ccGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGATCTGGCCGGGTTGCGGAAGGTGGCTCCGATCGTCTGATCGAACGAGCTGGCCGCGGTGCGGAAgg                                                                                               >  2:1476178/1‑150 (MQ=11)

                                                         cgcgTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGACTGATCGAACGTTCTGGCCGTGTTGCGGAAGGAGGCTCCGATCGTCTGATCGAACAATCTGGCCGAGTTGCTGAAgga                                                                                              >  2:6160307/1‑149 (MQ=11)

                                                              tGCGGAAGGTGGCTCCGCTCGCCTGATCGAACGTTCGGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTc                                                                                         <  2:4870149/150‑1 (MQ=255)

                                                                cTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGCCTggtt                                                                                                                         <  1:4111585/116‑4 (MQ=12)

                                                                cTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCAGATCGACTGATCGAACGTTCTGGCCGTGTTGCTGAAGGAGGCTCCGATCGTCTGATCGAACATTCTGGCCGCGTTGCTGAAGGAGGCTCCg                                                                                       >  1:312542/1‑150 (MQ=14)

                                                                 tGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGa                                                                                      <  2:6325333/150‑1 (MQ=21)

                                                                   aaGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGTTCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATc                                                                                    >  2:746295/1‑150 (MQ=255)

                                                                   aaGGTGGCTCCGAGCGCCTGATCGAACGTTCTTGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATc                                                                                    <  2:3243328/150‑1 (MQ=255)

                                                                      gcggCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTc                                                                                 <  2:1420047/148‑1 (MQ=25)

                                                                       cggCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCt                                                                                >  2:4917351/2‑150 (MQ=25)

                                                                        ggCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACTTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCGGAAGGTGGCTCCGATCGTCTg                                                                               >  2:390138/1‑150 (MQ=255)

                                                                          cTCCGATCGCCTGATCGAACGTCCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGAt                                                                             <  2:587923/150‑1 (MQ=37)

                                                                           tCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACTATCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATc                                                                            >  1:520415/1‑150 (MQ=255)

                                                                           tCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATc                                                                            >  1:2434513/1‑150 (MQ=37)

                                                                           tCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATc                                                                            >  1:6050405/1‑150 (MQ=37)

                                                                             cGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGAGCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGa                                                                          >  2:4973986/1‑150 (MQ=255)

                                                                                 cGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACga                                                                      >  1:1033145/1‑149 (MQ=37)

                                                                                  gCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGtt                                                                     >  2:1155247/1‑150 (MQ=37)

                                                                                   ccTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGATc                                                                    >  2:4734268/1‑150 (MQ=37)

                                                                                   ccTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGATc                                                                    >  2:4358847/1‑150 (MQ=37)

                                                                                      gttCGAACGGCCGGGCCGCGTTGCTGCAGGGGGCTCCGATCGCCTGATCGAACGTTCGGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTgg                                                                 <  2:2654734/148‑1 (MQ=11)

                                                                                       ttCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGc                                                                <  2:4046379/149‑1 (MQ=37)

                                                                                       aTCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACTGACAGGCCGCGTGGCTGAAGGTGGCTCCGATCGACAGATCGAACTATCTGGc                                                                >  2:657896/1‑150 (MQ=11)

                                                                                        tCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGcc                                                               <  2:6727610/150‑1 (MQ=37)

                                                                                          gAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCgt                                                             >  2:1593963/1‑150 (MQ=38)

                                                                                          gAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCgg                                                             >  2:4638634/1‑149 (MQ=37)

                                                                                             cgtccgggCCGAGTTGCTGACGGGGGCTCCGATCGCCTGCTCGACCGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGtt                                                          <  2:6031241/144‑1 (MQ=11)

                                                                                              gTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTg                                                         <  2:1312803/150‑1 (MQ=21)

                                                                                                  tGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGa                                                     <  2:908342/150‑1 (MQ=255)

                                                                                                  gggCCGCGTTGCTGCAGGGGGCGCCGATCGCCTGATCGAACGGCCGGGCCGCGTTGCTGAAGGGGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGa                                                     <  2:6436912/149‑1 (MQ=11)

                                                                                                   gcccGCGTAGCTGAAGGTGGCTCCGATCGCCTGCTCGAACGTTCTGGCCGAGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGaa                                                    <  2:3006095/148‑1 (MQ=255)

                                                                                                    gCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAg                                                   >  1:5082883/1‑150 (MQ=255)

                                                                                                    gCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAg                                                   <  2:4611699/150‑1 (MQ=255)

                                                                                                     ccGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAgg                                                  <  1:1368512/150‑1 (MQ=255)

                                                                                                      cgcgTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACTTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGAATGATCGAACGTTCTGGCCGTGTTGCCGAAGGt                                                 >  1:6490425/1‑150 (MQ=255)

                                                                                                      cgcgTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGt                                                 <  1:2272523/150‑1 (MQ=255)

                                                                                                           tGCTGAAGGGGGCTCCGATCGCCTGCTCGAACGGTCGGGCCGGGTTGCTGAAGGGGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGGTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTc                                            <  2:5274045/150‑1 (MQ=11)

                                                                                                              tGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCGGa                                         >  2:5328329/1‑150 (MQ=255)

                                                                                                                        tCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCGGATCGTCTGATc                               >  1:5619371/1‑150 (MQ=255)

                                                                                                                        tCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCGGATCGTCTGATc                               <  2:312061/150‑1 (MQ=37)

                                                                                                                          cGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCGGATCGTCTGATCGa                             >  2:19985/1‑150 (MQ=255)

                                                                                                                               gCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCGGATCGTCTGATCGAACAAg                        >  2:1818610/1‑150 (MQ=255)

                                                                                                                                ccTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACTATCTGGCCGTGTTGCCGAAGGTGGCTCGGATCGTCTGATCGAACAAGc                       >  2:2871528/1‑150 (MQ=255)

                                                                                                                                  tGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCGGATCGTCTGATCGAACAAGCTg                     <  1:798162/150‑1 (MQ=255)

                                                                                                                                     tCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCGGATCGTCTGATCGAACAAGCTGGcc                  >  1:1903395/1‑150 (MQ=255)

                                                                                                                                       gAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGATCTGGCCGCGTTGCTGAAGGAGGCTCCGAAAGTCTGATCGAACGATCAGGCCGTGTTGCCGAAGGTGGCTCGGATCGTCTGATCGAACAATCTGGCCgt                >  2:2558259/1‑150 (MQ=14)

                                                                                                                                        aaCGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCGGATCGTCTGATCGAACAAGCTGGCCgtg               >  1:5590044/1‑150 (MQ=255)

                                                                                                                                          cGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCGGATCGTCTGATCGAACAAGCTGGCCGTGtt             <  1:2365805/150‑1 (MQ=255)

                                                                                                                                          cGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCGGATCGTCTGATCGAACAAGCTGGCCGTGtt             >  2:6552250/1‑150 (MQ=255)

                                                                                                                                           gTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCGGATCGTCTGATCGAACAAGCTGGCCGTGTTg            >  2:3564579/1‑150 (MQ=255)

                                                                                                                                            ttCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCGGATCGTCTGATCGAACAAGCTGGCCGTGTTGc           >  1:327463/1‑150 (MQ=255)

                                                                                                                                            ttCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCGGATCGTCTGATCGAACAAGCTGGCCGTGTAGc           >  2:2767692/1‑150 (MQ=255)

                                                                                                                                              cTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCGGATCGTCTGATCGAACAATCTGGCCGTGTTGCCg         >  2:4287880/1‑150 (MQ=255)

                                                                                                                                              cTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCGGATCGTCTGATCGAACAAGCTGGCCGTGTTGCCg         >  2:5471244/1‑150 (MQ=255)

                                                                                                                                               tGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCGGATCGTCTGATCGAACAAGCTGGCCGTGTTGCCGa        >  2:5652960/1‑150 (MQ=255)

                                                                                                                                                ggCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCGGATCGTCTGATCGAACAAGCTGGCCGTGTTGCCGaa       <  1:4917351/150‑1 (MQ=255)

                                                                                                                                                 gCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCgt                                                             <  2:5948892/95‑1 (MQ=40)

                                                                                                                                                 gCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCgt                                                             >  1:5948892/1‑95 (MQ=40)

                                                                                                                                                 gCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCGGATCGTCTGATCGAACAAGCTGGCCGTGTTGCCGAAg      >  1:4597859/1‑150 (MQ=255)

                                                                                                                                                 gCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCGGATCGTCTGATCGAACAAGCTGGCCGTGTTGCCGAAg      >  1:3660892/1‑150 (MQ=255)

                                                                                                                                                    gtgtTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGACTCGGATCGTCTGATCGAACAAGCTGGCCGTGTTGCCGAAggcg   <  2:4216794/150‑1 (MQ=255)

                                                                                                                                                     tgtTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACTTTCTGGCCGTGTTGCCGAAGGTGGCTCGGATCGTCTGATCGAACTAGCTGGCCGTGTTGCCGAAggcgg  >  2:2244941/1‑150 (MQ=255)

CAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCGGATCGTCTGATCGAACAAGCTGGCCGTGTTGCCGAAGGCGG  >  NC_002947/5530327‑5530625

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