BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,154,439	C→A	7.3%	G126W (GGG→TGG)	PP_1909 ←	septum formation protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,154,439	0	C	A	7.3%	382.7 / 11.5	120	G126W (GGG→TGG)	PP_1909	septum formation protein
Reads supporting (aligned to +/- strand): ref base C (62/49); new base A (5/4); total (67/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 125 total aligned reads displayed.

TCTGGAACAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGCAACCTGGTCCGAGCCTATGATC  >  NC_002947/2154292‑2154587

tctGGAACAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCgg                                                                                                                                                    >  1:1349355/1‑150 (MQ=255)

 ctGGAACAAGCTCACCCACAGCCCCTCCGCCTTGGAACTCCCTTCTCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGt                                                                                                                                                   <  2:843229/150‑1 (MQ=255)

  tGGAACAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTg                                                                                                                                                  >  2:1277768/1‑150 (MQ=255)

  tGGAACAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTg                                                                                                                                                  >  2:562093/1‑150 (MQ=255)

   ggAACAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTgg                                                                                                                                                 <  1:562093/150‑1 (MQ=255)

     aaCAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGgcg                                                                                                                                               <  1:578207/150‑1 (MQ=255)

     aaCAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGgcg                                                                                                                                               >  2:242378/1‑150 (MQ=255)

     aaCAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCCCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGgcg                                                                                                                                               >  2:109484/1‑150 (MQ=255)

       ctaGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCGCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCt                                                                                                                                             <  1:1081423/148‑1 (MQ=255)

          gCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTTCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCGGGCAATGCCAGGt                                                                                                                                                   <  1:231769/141‑1 (MQ=255)

          gCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCCTCGCCAGCTGACGCTCtt                                                                                                                                          >  2:1081423/1‑150 (MQ=255)

          gCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCGGGCAATGCCAGGt                                                                                                                                                   >  2:231659/1‑141 (MQ=255)

          gCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCGGGCAATGCCAGGt                                                                                                                                                   <  1:231659/141‑1 (MQ=255)

            tCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTca                                                                                                                                        >  1:2209683/1‑150 (MQ=255)

             cACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCCAGGTCAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTcag                                                                                                                                       <  1:540957/150‑1 (MQ=255)

              aCCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTcagc                                                                                                                                      >  1:231085/1‑150 (MQ=255)

               cccccAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTcagca                                                                                                                                     <  1:1153960/150‑1 (MQ=255)

                 cccAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGt                                                                                                                                   <  2:2209683/150‑1 (MQ=255)

                   taGCCCCTCCTCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGc                                                                                                                                 <  1:1257991/149‑1 (MQ=255)

                     gCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCa                                                                                                                               <  1:16119/150‑1 (MQ=255)

                      ccccTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCCGCTCCCTGCTCAGTTCCCGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCACTGCCCGGTGGCACTGCTGAACAGCGCCac                                                                                                                              <  1:683820/150‑2 (MQ=255)

                        ccTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGcc                                                                                                                            >  1:907956/1‑150 (MQ=255)

                        ccTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCGCGGTTAATGCGCTGGCCCTGCAGGGGGTGGGCAGGCAACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCGGcc                                                                                                                            >  1:860718/1‑150 (MQ=255)

                         cTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGccc                                                                                                                           >  1:2195201/1‑150 (MQ=255)

                           ccGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCCCGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCgg                                                                                                                         <  1:728948/150‑1 (MQ=255)

                           ccGCCTTGAAACTCCCTGCGCAATACGGCTTCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCgg                                                                                                                         <  1:1170462/150‑1 (MQ=255)

                             gCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTca                                                                                                                                        <  1:2293896/133‑1 (MQ=255)

                             gCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTca                                                                                                                                        >  2:2293896/1‑133 (MQ=255)

                                ttGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGc                                                                                                                                 >  2:399533/1‑137 (MQ=255)

                                ttGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGc                                                                                                                                 <  1:399533/137‑1 (MQ=255)

                                  gAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGaaa                                                                                                                  >  2:1568216/1‑150 (MQ=255)

                                        cccTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCAc                                                                                                            <  2:1849732/150‑1 (MQ=255)

                                               cAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCGGCCCGGTAAGAAAGCTCACCtgctgc                                                                                                     <  2:278349/150‑1 (MQ=255)

                                               cAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCGGCCCGGTAAGAAAGCTCACCtgctgc                                                                                                     >  1:278349/1‑150 (MQ=255)

                                               cAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCtgctgc                                                                                                     <  1:1957322/150‑1 (MQ=255)

                                                        ggCTGTTCTGCTGCCACGTAACGCTCGACCGGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCTGTGGCGCTGTGCAGCCTTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTCCAc                                                                                            >  2:731278/1‑150 (MQ=255)

                                                           tGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTcc                                                                                         >  1:271804/1‑150 (MQ=255)

                                                            gTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCa                                                                                        <  1:1951818/150‑1 (MQ=255)

                                                            gTTCTGCTGCCACGTAACGCCCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCa                                                                                        <  1:923993/150‑1 (MQ=255)

                                                                tgctgcCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTc                                                                                          >  1:1207437/1‑144 (MQ=255)

                                                                tgctgcCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTc                                                                                          <  2:1207437/144‑1 (MQ=255)

                                                                  ctgcCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGTAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGtt                                                                                  >  1:2223090/1‑150 (MQ=255)

                                                                   tgcCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTg                                                                                 >  2:376812/1‑150 (MQ=255)

                                                                   tgcCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCACTGCACTCCAGTAGTTg                                                                                 >  1:333467/1‑150 (MQ=255)

                                                                       aCGTAACGCTCGACCTTCTCCCTGCTTAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAACGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCg                                                                             >  2:1060180/1‑150 (MQ=255)

                                                                         gTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGt                                                                                                                                   >  2:834411/1‑94 (MQ=255)

                                                                          tAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCaa                                                                          <  1:1015432/150‑1 (MQ=255)

                                                                                  gACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTc                                                                  >  1:1417454/1‑150 (MQ=255)

                                                                                  gACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTc                                                                  <  1:1568216/150‑1 (MQ=255)

                                                                                    cctgctccctgctcAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCGGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGa                                                                >  1:1041539/1‑150 (MQ=255)

                                                                                        ctccctgctcAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGt                                                            >  1:2200274/1‑150 (MQ=255)

                                                                                         tccctgctcAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGaaa                                                                                                                  <  2:1406295/95‑1 (MQ=255)

                                                                                         tccctgctcAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGaaa                                                                                                                  >  1:1406295/1‑95 (MQ=255)

                                                                                           cctgctcAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATg                                                         >  2:1218859/1‑150 (MQ=255)

                                                                                           cctgctcAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATg                                                         >  2:396644/1‑150 (MQ=255)

                                                                                           cctgctcAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATg                                                         >  1:1597062/1‑150 (MQ=255)

                                                                                               ctcAGTTCCCGCCGTGCCACGGCAAATGACTCGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCACTATTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGc                                                     <  1:514330/150‑1 (MQ=255)

                                                                                               ctcAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGc                                                     >  1:1750148/1‑150 (MQ=255)

                                                                                               ctcAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCCCGAACGTATGCGGc                                                     >  2:54916/1‑150 (MQ=255)

                                                                                                tcAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGtt                                                                                  <  1:747408/120‑1 (MQ=255)

                                                                                                tcAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCt                                                    >  1:1402399/1‑150 (MQ=255)

                                                                                                tcAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCt                                                    >  2:1069376/1‑150 (MQ=255)

                                                                                                 cAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCtt                                                   <  1:396644/150‑1 (MQ=255)

                                                                                                 cAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCtt                                                   <  1:1218859/150‑1 (MQ=255)

                                                                                                   gTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCATGCGCGCTCGAACGTATGCGGCTTGc                                                 >  2:1706104/1‑150 (MQ=255)

                                                                                                   gTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGc                                                 <  2:271804/150‑1 (MQ=255)

                                                                                                      cACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCa                                              <  2:1402399/150‑1 (MQ=255)

                                                                                                       aCGCAGTGTCACGGTCACAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAg                                             >  2:208852/1‑150 (MQ=255)

                                                                                                       aCGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAg                                             >  2:1410950/1‑150 (MQ=255)

                                                                                                           aGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCAcc                                         <  1:1620406/150‑1 (MQ=255)

                                                                                                           aGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCAcc                                         >  1:1006075/1‑150 (MQ=255)

                                                                                                              gtCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGc                                      <  2:1417454/150‑1 (MQ=255)

                                                                                                                cACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTc                                    <  1:376812/150‑1 (MQ=255)

                                                                                                                   ggTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCTCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGcc                                 <  2:1349355/150‑1 (MQ=255)

                                                                                                                   ggTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCACGGGCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGcc                                 >  1:651202/1‑150 (MQ=255)

                                                                                                                      aaaaGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCaa                              >  1:1692316/1‑150 (MQ=255)

                                                                                                                          ggCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATAcc                          <  2:2200274/150‑1 (MQ=255)

                                                                                                                            cACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCa                                                                                        <  2:491870/86‑1 (MQ=255)

                                                                                                                            cACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGc                        <  2:1692316/150‑1 (MQ=255)

                                                                                                                             aCGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTg                                                  <  2:1130298/123‑1 (MQ=255)

                                                                                                                             aCGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTg                                                  >  1:1130298/1‑123 (MQ=255)

                                                                                                                             aCGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGCa                       >  2:2305209/1‑150 (MQ=255)

                                                                                                                             aCGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGACCACAACAGTCCAGTAGCTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGCa                       >  2:1134920/1‑150 (MQ=255)

                                                                                                                              cGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCAc                                          >  2:940521/1‑130 (MQ=255)

                                                                                                                              cGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCAc                                          <  1:940521/130‑1 (MQ=255)

                                                                                                                                 aaTCCACCTGGCATTTCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCt                                                                                                >  2:338562/1‑73 (MQ=255)

                                                                                                                                 aaTCCACCTGGCATTTCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCt                                                                                                <  1:338562/73‑1 (MQ=255)

                                                                                                                                 aaTCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGCAACCt                   >  2:2279702/1‑150 (MQ=255)

                                                                                                                                 aaTCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGCAACCt                   >  2:177244/1‑150 (MQ=255)

                                                                                                                                 aaTCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGCAACCt                   >  2:1109974/1‑150 (MQ=255)

                                                                                                                                  aTCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGCAACCTg                  <  2:2223090/150‑1 (MQ=255)

                                                                                                                                   gccACCGGGCAATGCCAGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTc                                                                              >  2:1520763/2‑89 (MQ=255)

                                                                                                                                     cACCGGGCAATGCCAGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCa                                                                                        >  2:409641/1‑77 (MQ=255)

                                                                                                                                     cACCGGGCAATGCCAGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCa                                                                                        <  1:409641/77‑1 (MQ=255)

                                                                                                                                      aCCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGCAACCTGGTcc              >  2:1758612/1‑150 (MQ=255)

                                                                                                                                        cTGGCATTGCCCGGTGGCGCGGTTCAGCAGTGCCATCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGCAACCTGGTCCGa            <  2:1643997/150‑1 (MQ=255)

                                                                                                                                             aTTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGTCTTGCCCAGCACCTGCTCGCCCAATACCGCAACCTGGTCCGAGCCTa       <  2:907956/150‑1 (MQ=255)

                                                                                                                                                  ccGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTc                                    >  1:114292/1‑116 (MQ=255)

                                                                                                                                                  ccGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTc                                    <  2:114292/116‑1 (MQ=255)

                                                                                                                                                  ccGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGCAACCTGGTCCGAGCCTATGATc  <  1:110702/150‑1 (MQ=255)

TCTGGAACAAGCTCACCCCCAGCCCCTCCGCCTTGAAACTCCCTGCGCAATCCAGCGGCTGTTCTGCTGCCACGTAACGCTCGACCTGCTCCCTGCTCAGTTCACGCAGTGTCACGGTAAAAGGCACGCAATCCACCTGGCATTGCCCGGTGGCGCTGTTCAGCAGTGCCAGCCCGGTAAGAAAGCTCACCTGCTGCCCGCTGCACTCCAGTAGTTGCTCGCAAGCGCGCTCGAACGTATGCGGCTTGCCCAGCACCTGCTCGCCCAATACCGCAACCTGGTCCGAGCCTATGATC  >  NC_002947/2154292‑2154587

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