BRESEQ :: Evidence

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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,530,455	0	C	T	5.1%	793.6 / 18.5	284	R82R (CGC→CGT)	PP_4858	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (101/153); new base T (10/3); total (127/157)
Fisher's exact test for biased strand distribution p-value = 1.74e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 296 total aligned reads displayed.

TCCAGCGCTGCCGAAATCGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCGGATCGTCTGATCGAACAAGC  >  NC_002947/5530306‑5530604

tCCAGCGCTGCCGAAATCGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGc                                                                                                                                                       <  2:632565/150‑1 (MQ=255)

tCCAGCGCTGCCGAAATCGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGc                                                                                                                                                       <  2:7087839/150‑1 (MQ=255)

  cAGCGCTGCCGAAATCGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCt                                                                                                                                                     <  1:6116580/150‑1 (MQ=255)

    gcgcTGCCGAAATCGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGa                                                                                                                                                   <  1:8098529/150‑1 (MQ=255)

    gcgcTGCCGAAATAGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGa                                                                                                                                                   <  2:6302485/150‑1 (MQ=255)

     cgcTGCCGAAATCGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGAt                                                                                                                                                  <  2:7847144/150‑1 (MQ=255)

     cgcTGCCGAAATCGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGAt                                                                                                                                                  <  2:1238366/150‑1 (MQ=255)

      gcTGCCGAAATCGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATc                                                                                                                                                 <  1:7848497/150‑1 (MQ=255)

         gCCGAAATCGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGACTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGaa                                                                                                                                              >  1:7319326/1‑150 (MQ=255)

          ccGAAATCGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCGGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAAc                                                                                                                                             <  2:2094585/150‑1 (MQ=255)

            gAAATCGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGt                                                                                                                                           >  1:7037085/1‑150 (MQ=255)

             aaaTCGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGtt                                                                                                                                          <  1:2403881/150‑1 (MQ=255)

              aaTCTCCCAGCCTCTGATCACCGTAGCTGTAGGCGGCTCCGATGGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTc                                                                                                                                         <  2:8041406/150‑1 (MQ=17)

                tCGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTg                                                                                                                                       >  2:4594901/1‑150 (MQ=255)

                  gCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGc                                                                                                                                     <  2:3842823/150‑1 (MQ=255)

                  gCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGACTGATCGAACGTTCTGGCCGCGTTGCAGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGc                                                                                                                                     >  1:2894187/1‑150 (MQ=255)

                    ccAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGAt                                                                                                                                                  <  1:4988178/135‑1 (MQ=255)

                    ccAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCg                                                                                                                                   <  2:3658687/150‑1 (MQ=255)

                          gCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGc                                                                                                                             >  2:4980671/1‑150 (MQ=255)

                          gCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGc                                                                                                                             >  1:1641466/1‑150 (MQ=255)

                          gCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCGGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGc                                                                                                                             <  1:1786023/150‑1 (MQ=255)

                            tGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTg                                                                                                                           <  1:5693591/150‑1 (MQ=255)

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

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

                                    gTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGc                                                                                                                   >  1:884487/1‑150 (MQ=255)

                                           caGGCGGCTCCGAGCGCCTGATCGATCGTTCGGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATc                                                                                                            <  2:1831447/149‑1 (MQ=255)

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

                                                ggcTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCGGAAGGTGGCTCCGATCGTCTg                                                                                                       >  2:7660883/1‑150 (MQ=37)

                                                   tCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGATCTGGCCGGGTTGCTGACGGTGGCTCCGATCGGCTGATc                                                                                                    >  1:5397316/1‑150 (MQ=255)

                                                     cGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGa                                                                                                  >  1:4820381/1‑150 (MQ=37)

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

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

                                                                   aaCGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCcgcg                                                                                    >  1:6628196/1‑150 (MQ=38)

                                                                   aaCGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCcgcg                                                                                    <  2:4130919/150‑1 (MQ=38)

                                                                    aCGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGt                                                                                   <  1:8589386/150‑1 (MQ=39)

                                                                           ggCCGCGTTGCTGAATGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGaa                                                                            <  2:1036615/150‑1 (MQ=255)

                                                                            gCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCGGATCGCCt                                                                                                                                                     <  2:8703826/76‑1 (MQ=255)

                                                                            gCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGCAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCGGacg                                                                           >  1:1057311/1‑148 (MQ=255)

                                                                            gCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAg                                                                           <  2:5329136/150‑1 (MQ=255)

                                                                             ccGAGTTGCTGCAGGCGGCTCCGATCGCCTGCTCGACCGGCCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAgg                                                                          <  2:8496323/150‑1 (MQ=11)

                                                                              cgcgTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGt                                                                         <  1:1365362/150‑1 (MQ=255)

                                                                                cgTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGTTGAAGGTgg                                                                       <  1:5540458/150‑1 (MQ=255)

                                                                                   tGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAgggacagg                                                                    <  2:4529933/150‑7 (MQ=255)

                                                                                   agCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTc                                                                    <  2:4173523/149‑1 (MQ=38)

                                                                                     cTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCg                                                                  <  2:5230003/150‑1 (MQ=21)

                                                                                      tGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGa                                                                 >  1:3066545/1‑150 (MQ=39)

                                                                                      tGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGACTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGGGTTGCTGAAGGTGGCTCCGa                                                                 >  2:3353350/1‑150 (MQ=255)

                                                                                        aaGGTGGCTCCGATCGCCTGATCGACCGGTCTCGCCGCGTTGCTGAAGGGGGCTCCGATCGCCTGATCGACCGGTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATc                                                               <  2:7431899/150‑1 (MQ=11)

                                                                                        aaGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATc                                                               <  1:7909809/150‑1 (MQ=39)

                                                                                        aaGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATc                                                               <  2:5697589/150‑1 (MQ=39)

                                                                                           gTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTc                                                            <  1:8437223/150‑1 (MQ=38)

                                                                                           gTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTc                                                            <  1:7556274/150‑1 (MQ=38)

                                                                                           gTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTc                                                            >  2:3985557/1‑150 (MQ=38)

                                                                                           gTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTc                                                            >  1:19824/1‑150 (MQ=38)

                                                                                            tGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGACTGATCGAACGTTCTGGCCGGGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCGGAAGGTGGCTCCGATCGTCt                                                           >  2:1730137/1‑150 (MQ=255)

                                                                                              gCTCCGATCGCCTGATCGAACGTTCTTGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGa                                                         <  1:5183601/150‑1 (MQ=37)

                                                                                              gCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGa                                                         <  2:707443/150‑1 (MQ=37)

                                                                                              gCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGa                                                         <  2:5557695/150‑1 (MQ=37)

                                                                                              gCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGa                                                         <  2:549593/150‑1 (MQ=37)

                                                                                                tCCGATCGCCTGTTCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATc                                                       <  2:2273263/150‑1 (MQ=17)

                                                                                                    aTCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGGGTGGCTGAAGGTGGCTCCGATCGTCTGATCGAAc                                                   >  2:2613188/1‑150 (MQ=21)

                                                                                                       gCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGtt                                                <  1:1567591/150‑1 (MQ=37)

                                                                                                        ccTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTc                                               >  1:4648123/1‑150 (MQ=37)

                                                                                                        ccTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTc                                               >  2:575088/1‑150 (MQ=37)

                                                                                                         cTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCt                                              >  2:2889835/1‑150 (MQ=37)

                                                                                                           gttCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGGGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTgg                                            <  2:675809/148‑1 (MQ=21)

                                                                                                            aTCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGACTGATAGAACGATCAGGc                                           >  1:7185918/1‑150 (MQ=255)

                                                                                                               gAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCgt                                        <  1:126338/150‑1 (MQ=38)

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

                                                                                                                aaCGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTc                                                                                                                 >  1:6380049/1‑76 (MQ=34)

                                                                                                                aaCGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCgtg                                       >  1:34027/1‑150 (MQ=38)

                                                                                                                 aCGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACTTACTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCgtgt                                      >  2:4296638/1‑150 (MQ=255)

                                                                                                                  cGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGtt                                     <  2:7213824/150‑1 (MQ=39)

                                                                                                                     tCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGcc                                  <  2:7037085/150‑1 (MQ=40)

                                                                                                                       tGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGa                                >  2:8096376/1‑150 (MQ=255)

                                                                                                                        ggCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGaa                               >  2:5914941/1‑150 (MQ=255)

                                                                                                                         gcagcgTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAg                              >  1:3532552/4‑150 (MQ=255)

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

                                                                                                                           cgcgTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTAGCCGAAGGt                            >  2:3549008/1‑150 (MQ=255)

                                                                                                                           cgcgTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCAGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGACGAAGGt                            >  2:3177955/1‑150 (MQ=255)

                                                                                                                            gcgTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCt                                                                               <  2:4447575/98‑1 (MQ=40)

                                                                                                                               ttGCTGAAGGTGGCTCCGATCGACTGTTCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCt                        <  2:1445976/150‑1 (MQ=255)

                                                                                                                                tGCTGTAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTc                       <  2:6622272/150‑1 (MQ=255)

                                                                                                                                  cTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCgg                     >  1:7423731/1‑150 (MQ=255)

                                                                                                                                  cTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGGGTTGCGGAAGGTGGCTCCGATCGTCGGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCgg                     >  2:8164488/1‑150 (MQ=255)

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

                                                                                                                                   tGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTAGCTGAAGGTGGCTCGGATCGTCTGGTCGAACTATCGGGCCTGGTAGCAGCCGGTGTGTCGGCTCGACTGATCGACGGTTCTCGCGGTGTTGCCGAAGGTGGCTCGGa                    >  2:8327243/1‑150 (MQ=9)

                                                                                                                                   tGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTAGCTGAAGGTGGCTCCGATCGTCTGATCGAACTTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGGTCTGGCCGTGTTGCCGAAGGTGGCTCgta                    >  1:828509/1‑148 (MQ=255)

                                                                                                                                    gAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACTTTCTGGACGCGTTGCTGAAGGTGGATCCGATCGTCTGATCGAACGTTCTGGGCGTGTTGAAGAAGGTGGCTCGGAt                   >  1:3453161/1‑150 (MQ=14)

                                                                                                                                     aaGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCGGATc                  <  1:2288015/150‑1 (MQ=21)

                                                                                                                                      aGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCGGATCg                 <  1:5920447/150‑1 (MQ=255)

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

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

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

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

                                                                                                                                               cGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGACTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCGGATCGTCTGATCGa        <  2:4385983/150‑1 (MQ=21)

                                                                                                                                                gATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCGGATCGTCTGATCGaa       <  2:6514179/150‑1 (MQ=255)

                                                                                                                                                     ccTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCGGATCGTCTGATCGAACAAgt  >  2:2391463/1‑149 (MQ=255)

                                                                                                                                                     ccTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCGGATCGTCTGATCGAACAAGc  >  1:4392198/1‑150 (MQ=255)

                                                                                                                                                     ccTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCAGAAGGTGGCTCGGATCGTCTGATCGAACAAGc  >  2:1379922/1‑150 (MQ=255)

TCCAGCGCTGCCGAAATCGCCCAGCCGCTGAACACCGTAGCTGAAGGCGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGCCTGATCGAACGTTCTGGCCGTGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGCGTTGCTGAAGGTGGCTCCGATCGTCTGATCGAACGTTCTGGCCGTGTTGCCGAAGGTGGCTCGGATCGTCTGATCGAACAAGC  >  NC_002947/5530306‑5530604

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