BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	508,775	A→T	6.2%	V79D (GTC→GAC)	estP ←	esterase EstP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	508,775	0	A	T	6.2%	268.2 / 4.5	91	V79D (GTC→GAC)	estP	esterase EstP
Reads supporting (aligned to +/- strand): ref base A (42/42); new base T (4/2); total (47/44)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

ATTGGGGTTTCCTGGGGGGATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGG  >  NC_002947/508630‑508921

atgggggTTTCCTGGGGGGATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGGGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGc                                                                                                                                                  <  3:201075/146‑1 (MQ=255)

aTTGGGGTTTCCTGGGGGGATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTCCCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGc                                                                                                                                                  <  2:76383/149‑1 (MQ=255)

 ttGGGGTTTCCTGGGGGGATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGc‑c                                                                                                                                                <  2:179568/149‑1 (MQ=255)

 ttGGGGTTTCCTGGGGGGATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCAGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGc‑c                                                                                                                                                >  1:330511/1‑149 (MQ=255)

    gggTTTCCTGGGGGGATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAg                                                                                                                                             >  1:71034/1‑149 (MQ=255)

     ggTTTCCTGGGGGGATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAgg                                                                                                                                            <  5:95722/149‑1 (MQ=255)

      gTTTCCTGGGGGGATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAggc                                                                                                                                           <  5:149885/149‑1 (MQ=255)

       tttCCTGGGGGGATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAggcg                                                                                                                                          <  6:397897/149‑1 (MQ=255)

          ccTGGGGGGATGACGGTCTCGGACGTGGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAggcggc                                                                                                                                        <  7:153024/148‑1 (MQ=255)

           cTGGGGGGATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGCCGCCTGG‑GCGTCGCGCCAGGCGAcc                                                                                                                                       <  5:223138/148‑1 (MQ=255)

              gggggATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCAc                                                                                                                                     <  6:331125/147‑1 (MQ=255)

                gggATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAg                                                                                                                                 <  1:408036/149‑1 (MQ=255)

                  gATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGa                                                                                                                               >  7:359174/1‑149 (MQ=255)

                        ggTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCg                                                                                                                         >  6:307821/1‑149 (MQ=255)

                        ggTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACACGCCGACAGCCCAGTTAATACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAACGTCGCCTGGCGCGACGC‑CCAGGCCGCCACAGAGGATCATCg                                                                                                                         >  6:276283/1‑149 (MQ=255)

                             cGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGac                                                                                                                    <  7:180401/149‑1 (MQ=255)

                              ggACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGCGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGaca                                                                                                                   >  8:29297/1‑149 (MQ=255)

                                aCGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACAcc                                                                                                                 <  5:307821/149‑1 (MQ=255)

                                     gTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGc                                                                                                            <  4:223497/149‑1 (MQ=255)

                                        gTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCg                                                                                                           >  4:213516/1‑147 (MQ=255)

                                         tGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAg                                                                                                        <  2:71034/149‑1 (MQ=255)

                                         tGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACCGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCCGGGCGAAg                                                                                                        >  8:190971/1‑149 (MQ=255)

                                           aTCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGtt                                                                                                      >  6:65922/1‑149 (MQ=255)

                                            tCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTg                                                                                                     >  8:230244/1‑149 (MQ=255)

                                                  ccAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCg                                                                                               >  1:33840/1‑149 (MQ=255)

                                                      gATCTGCTGGGTGGTGTACCCGCAGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCg                                                                                            >  8:323300/1‑148 (MQ=255)

                                                          tgctgGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACg                                                                                       >  6:124457/1‑149 (MQ=255)

                                                          tgctgGGTGGTGTACCCGCCGACTGCCCAGTTATTACAATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACg                                                                                       >  2:70676/1‑149 (MQ=255)

                                                           gctgGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACgg                                                                                      >  6:421512/1‑149 (MQ=255)

                                                           gctgGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACgg                                                                                      >  5:46682/1‑149 (MQ=255)

                                                            ctgGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCAc                                                                                        <  3:213516/146‑1 (MQ=255)

                                                            ctgGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACgg                                                                                      <  5:65922/148‑1 (MQ=255)

                                                            ctgGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGt                                                                                     >  3:54179/1‑149 (MQ=255)

                                                               ggtggtGTACTAATCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGtt                                                                                    <  1:151014/147‑1 (MQ=255)

                                                                gtggtgTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGt                                                                                 <  7:230244/149‑1 (MQ=255)

                                                                 tggtgTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGt                                                                                 >  6:283232/1‑148 (MQ=255)

                                                                  ggtgTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTg                                                                                >  8:122980/1‑148 (MQ=255)

                                                                   gtgtACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGCCACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGa                                                                               >  3:131843/1‑148 (MQ=255)

                                                                    tgtACCTGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCTAGGGGTTAAGGTCGCCTGGCGCGACGC‑CGAGGCGGCCACCGAGGATCATCGGCGACAGGGGGGCGAAGTTGCGGTCGGCGTCACGGTTGGTGaaa                                                                             <  8:127411/149‑1 (MQ=255)

                                                                    tgtACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGaaa                                                                             <  8:359174/149‑1 (MQ=255)

                                                                     gtACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAAc                                                                            >  3:386236/1‑149 (MQ=255)

                                                                               aCTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGTCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTgg                                                                  <  4:259719/149‑1 (MQ=255)

                                                                               aCTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTgg                                                                  <  8:393097/149‑1 (MQ=255)

                                                                                cTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTggg                                                                 <  5:318985/149‑1 (MQ=255)

                                                                                    ccAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGCCGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTAc                                                             >  2:21811/1‑149 (MQ=255)

                                                                                        ttattaCCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGt                                                                                 >  8:197812/1‑125 (MQ=255)

                                                                                        ttattaCCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGt                                                                                 <  7:197812/125‑1 (MQ=255)

                                                                                           ttaCCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTTGTGAAACGCGCGCCTGGGGTACCGCcaa                                                       <  2:330511/148‑1 (MQ=255)

                                                                                           ttaCCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCcaa                                                       >  8:107694/1‑148 (MQ=255)

                                                                                           ttaCCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCcaa                                                       <  3:2693/148‑1 (MQ=255)

                                                                                              ccATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCcaacaa                                                    <  1:74332/148‑1 (MQ=255)

                                                                                               cATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCcaac                                                      >  8:249566/1‑145 (MQ=255)

                                                                                                 tCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATc                                                <  7:249566/149‑1 (MQ=255)

                                                                                                    ggCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCggg                                             <  7:323300/149‑1 (MQ=255)

                                                                                                       tGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGaaa                                          <  4:416533/149‑1 (MQ=255)

                                                                                                       tGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGaaa                                          >  5:126033/1‑149 (MQ=255)

                                                                                                         gATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTAcc                                                            <  4:89295/129‑1 (MQ=255)

                                                                                                         gATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTAcc                                                            >  3:89295/1‑129 (MQ=255)

                                                                                                              cTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCg                                   >  7:398265/1‑149 (MQ=255)

                                                                                                              cTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCg                                   >  5:417901/1‑149 (MQ=255)

                                                                                                              cTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCg                                   >  4:383215/1‑149 (MQ=255)

                                                                                                               tACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTAc                                                             >  1:393320/1‑122 (MQ=255)

                                                                                                               tACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTAc                                                             <  2:393320/122‑1 (MQ=255)

                                                                                                               tACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCgg                                  >  6:313798/1‑149 (MQ=255)

                                                                                                                  ggATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCAt                               >  4:398725/1‑149 (MQ=255)

                                                                                                                      gTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGc                           >  2:45293/1‑149 (MQ=255)

                                                                                                                       tCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGc                           <  6:126033/148‑1 (MQ=255)

                                                                                                                        cGACGGGTTAAGGTCTCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGAACAACGGCGACACCGGGGCGAAGATGCCGTCGGCGTCACCGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGCGAAACTCCCCGCCATCGGc                           >  5:235116/1‑147 (MQ=255)

                                                                                                                         gACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGaaa                                          <  4:111204/131‑1 (MQ=255)

                                                                                                                         gACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGaaa                                          >  3:111204/1‑131 (MQ=255)

                                                                                                                         gACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAg                        >  7:42229/1‑149 (MQ=255)

                                                                                                                            gggTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGc                                                          <  5:49865/112‑1 (MQ=255)

                                                                                                                            gggTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGc                                                          >  6:49865/1‑112 (MQ=255)

                                                                                                                               ttAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGACACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTc                  >  3:252716/1‑149 (MQ=255)

                                                                                                                               ttAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTc                  <  7:160244/149‑1 (MQ=255)

                                                                                                                                ggtggccGCCTGG‑GCGTCGCGCCAGGCGACCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCg                 >  3:138821/7‑149 (MQ=255)

                                                                                                                                 aaGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCg                 <  4:54179/148‑1 (MQ=255)

                                                                                                                                  aGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGc                <  2:33840/148‑1 (MQ=255)

                                                                                                                                     ccGCCTGG‑GCGTCGCGCCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCaaa            >  4:99651/2‑149 (MQ=255)

                                                                                                                                     ccGCCTGG‑GCGTCGCGCCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCaaa            >  2:400603/2‑149 (MQ=255)

                                                                                                                                     ccGCCTGG‑GCGTCGCGCCAGGCGACCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGa                                            <  2:160588/116‑1 (MQ=255)

                                                                                                                                     ccGCCTGG‑GCGTCGCGCCAGGCGACCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGa                                            >  1:160588/2‑117 (MQ=255)

                                                                                                                                        ccTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGGCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTAACGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAacac         <  1:327086/149‑1 (MQ=255)

                                                                                                                                         cTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGTGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCt                     <  4:54450/136‑1 (MQ=255)

                                                                                                                                         cTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCt                     >  3:54450/1‑136 (MQ=255)

                                                                                                                                         cTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAacac         >  1:221984/1‑148 (MQ=255)

                                                                                                                                           ggCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGAt      >  6:412353/1‑149 (MQ=255)

                                                                                                                                             cgcgACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCa    <  1:45293/149‑1 (MQ=255)

                                                                                                                                              gcgACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATc                                                <  4:135275/104‑1 (MQ=255)

                                                                                                                                              gcgACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATc                                                >  3:135275/1‑104 (MQ=255)

                                                                                                                                               cgACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAgg  <  5:124457/149‑1 (MQ=255)

ATTGGGGTTTCCTGGGGGGATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGG  >  NC_002947/508630‑508921

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