BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,502,992	C→T	5.5%	T96M (ACG→ATG)	speB →	agmatinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,502,992	0	C	T	5.5%	338.6 / 5.3	109	T96M (ACG→ATG)	speB	agmatinase
Reads supporting (aligned to +/- strand): ref base C (53/50); new base T (2/4); total (55/54)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.35e-01

Only 100 of 111 total aligned reads displayed.

CGCCCCGGCGCGCGCTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGCCCGCTGGCGCTGA  >  NC_002947/2502850‑2503138

cGCCCCGGCGCGCGCTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCAcgccg                                                                                                                                                 >  8:291075/1‑147 (MQ=255)

 gCCCCGGCGCGCGCTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAg                                                                                                                                              <  6:152792/149‑1 (MQ=255)

  ccccGGCGCGCGCTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACTGCGGCGTGCCGCTGt                                                                                                                                             <  7:257626/149‑1 (MQ=255)

  ccccGGCGCGCGCTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGt                                                                                                                                             >  6:80675/1‑149 (MQ=255)

   cccGGCGCGCGCTTCGGCCCGCGGCCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTAGACAGCGGCACCCCGCAgac                                                                                                                                            >  6:75089/1‑147 (MQ=255)

    gcggcggGGGCTTCGGCCGGCGGGCGATCCGTGCTGCCTCGGTCCAGCTGGCCTGGGCCCGGCTCTGGCCGTGGGCGTTCGACCCGTTCGACCACGTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCTCGCCGCAGTCg                                                                                                                                           <  2:308008/148‑1 (MQ=255)

      ggCGCGCGCTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACTGCGGCGTGCCGCTGTCgaa                                                                                                                                         <  4:319731/149‑3 (MQ=255)

       gcgcgcgcTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACTGCGGCGTGCCGCTGTCgaag                                                                                                                                        <  4:270241/149‑4 (MQ=255)

       gcgcgcgcTTCGGCCCGCGGGCGATCCGCGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACTGCGGCGTGCCGCTGTCgaag                                                                                                                                        <  2:92561/149‑4 (MQ=255)

        cgcgcgcTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTAc                                                                                                                                       <  8:161881/149‑1 (MQ=255)

          cgcgcTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCg                                                                                                                                     <  2:11286/149‑1 (MQ=255)

          cgcgcTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCg                                                                                                                                     <  4:71465/149‑1 (MQ=255)

                 cGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGc                                                                                                                               <  2:11667/148‑1 (MQ=255)

                    cccGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGa                                                                                                                                   <  6:163985/141‑1 (MQ=255)

                    cccGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATc                                                                                                                            >  2:184589/1‑148 (MQ=255)

                      cgcgGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAg                                                                                                                         >  8:130059/1‑149 (MQ=255)

                       gcgGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTTGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAgg                                                                                                                        <  8:224383/149‑1 (MQ=255)

                       gcgGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAgg                                                                                                                        <  6:145149/149‑1 (MQ=255)

                              aTCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGc                                                                                                                 >  7:280476/1‑149 (MQ=255)

                              aTCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGc                                                                                                                 >  3:231733/1‑149 (MQ=255)

                                 cGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGTCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGa                                                                                                              <  7:46654/149‑1 (MQ=255)

                                 cGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGa                                                                                                              <  6:11990/149‑1 (MQ=255)

                                   tgctgcCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGc                                                                                                            <  4:231733/149‑1 (MQ=255)

                                      tgcCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGcaca                                                                                                         >  1:117566/1‑149 (MQ=255)

                                        ccTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCGGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACAt                                                                                                        >  8:2600/1‑148 (MQ=255)

                                          tCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCt                                                                                                     >  4:295935/1‑149 (MQ=255)

                                          tCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCt                                                                                                     >  1:111721/1‑149 (MQ=255)

                                           cGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCCCGCCGCAGTCGCTACCGGACAGCATCCAGGACCATGACGACCACATTCtc                                                                                                    >  5:18889/1‑148 (MQ=255)

                                           cGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTg                                                                                                    >  5:96666/1‑149 (MQ=255)

                                            ggTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGTCCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCt                                                                                                     <  8:225279/147‑1 (MQ=255)

                                            ggTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTg                                                                                                    >  5:334121/1‑148 (MQ=255)

                                             gTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTg                                                                                                    >  5:285874/1‑147 (MQ=255)

                                                  gcagGCCTGGGCCCGCCACAGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCg                                                                                                                                           <  7:285461/103‑1 (MQ=255)

                                                  gcagGCCTGGGCCCGCCACAGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCg                                                                                                                                           >  8:285461/1‑103 (MQ=255)

                                                      gCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTg                                                                                        >  1:299599/1‑149 (MQ=255)

                                                       ccTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTgc                                                                                       <  6:96666/149‑1 (MQ=255)

                                                        cTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGCGCACATTCTGCAGGC‑GGGTTgcg                                                                                      >  8:134643/1‑149 (MQ=255)

                                                        cTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTgcg                                                                                      <  1:230007/149‑1 (MQ=255)

                                                        cTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTgcg                                                                                      >  4:309345/1‑149 (MQ=255)

                                                          gggCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGcc                                                                                    <  8:100480/149‑1 (MQ=255)

                                                            gCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATg                                                                                                                  >  1:275593/1‑118 (MQ=255)

                                                            gCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATg                                                                                                                  <  2:275593/118‑1 (MQ=255)

                                                              ccGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGTGCACATTCTGCAGGC‑GGGTTGCGCCATGc                                                                                <  7:134643/149‑1 (MQ=255)

                                                                 ccACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGa                                                                             <  2:299599/149‑1 (MQ=255)

                                                                     tGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGc                                                                                                   >  8:305638/1‑124 (MQ=255)

                                                                     tGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGc                                                                                                   <  7:305638/124‑1 (MQ=255)

                                                                       gCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCtt                                                                        <  3:160871/148‑1 (MQ=255)

                                                                         cGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTtggt                                                                     <  3:309345/149‑1 (MQ=255)

                                                                          gTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTtggtg                                                                    <  2:223669/149‑1 (MQ=255)

                                                                            gggCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTtgg                                                                      >  3:180776/1‑145 (MQ=255)

                                                                            gggCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTtgg                                                                      <  4:180776/145‑1 (MQ=255)

                                                                                  tcgacccgttcgaccACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGc                                                                                               <  3:188584/115‑1 (MQ=255)

                                                                                  tcgacccgttcgaccACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGc                                                                                               >  4:188584/1‑115 (MQ=255)

                                                                                   cgacccgttcgaccACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTTGGTGGCGACCACt                                                           >  3:69495/1‑149 (MQ=255)

                                                                                    gacccgttcgaccACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTtggtg                                                                    <  5:55090/139‑1 (MQ=255)

                                                                                    gacccgttcgaccACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTtgggg                                                                    >  6:55090/1‑137 (MQ=255)

                                                                                    gacccgttcgaccACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTTGGTGGCGACCACt                                                           >  4:253980/1‑148 (MQ=255)

                                                                                     acccgttcgaccACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCCCATTCTGCAGGC‑GGGTTGCGCCATGCTGAGCCTTGGTGGCGACCACtt                                                          >  2:46325/1‑148 (MQ=255)

                                                                                     acccgttcgaccACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTtc                                                         <  5:80675/149‑1 (MQ=255)

                                                                                     acccgttcgaccACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTtc                                                         <  1:184589/149‑1 (MQ=255)

                                                                                     acccgttcgaccACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCAGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTTGGTGGCGACCACtt                                                          >  5:228124/1‑148 (MQ=255)

                                                                                      ccagttcgtcaACCTGGAGGTCATCGACTATGTCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTtca                                                        <  6:228124/146‑1 (MQ=255)

                                                                                        cgttcgaccACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTtcatc                                                      <  7:2600/149‑1 (MQ=255)

                                                                                         gttcgaccACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTtcat                                                       >  2:57574/1‑147 (MQ=255)

                                                                                         gttcgaccACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTtc                                                         >  7:230550/1‑145 (MQ=255)

                                                                                         gttcgaccACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTtc                                                         <  8:230550/145‑1 (MQ=255)

                                                                                             gaccACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCATGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTa                                                 >  3:142023/1‑149 (MQ=255)

                                                                                             gaccACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTa                                                 >  4:275111/1‑149 (MQ=255)

                                                                                                caccTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTTGGTGGCGa                                                                <  4:120652/131‑1 (MQ=255)

                                                                                                caccTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTTGGTGGCGa                                                                >  3:120652/1‑131 (MQ=255)

                                                                                                 accTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCg                                             <  7:306549/149‑1 (MQ=255)

                                                                                                 accTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCg                                             >  8:306549/1‑149 (MQ=255)

                                                                                                  ccTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCgc                                            >  6:257734/1‑149 (MQ=255)

                                                                                                      gCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCgctgct                                        >  6:308718/1‑149 (MQ=255)

                                                                                                             tCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCt                                 >  5:185968/1‑149 (MQ=255)

                                                                                                              cGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCt                                 >  3:62758/1‑148 (MQ=255)

                                                                                                               gACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTc                                <  2:111721/148‑1 (MQ=255)

                                                                                                               gACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTc                                <  8:280476/148‑1 (MQ=255)

                                                                                                                aCTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCAc                              <  3:253980/149‑1 (MQ=255)

                                                                                                                aCTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCAc                              >  8:173099/1‑149 (MQ=255)

                                                                                                                aCTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCAc                              >  6:284412/1‑149 (MQ=255)

                                                                                                                aCTATGGCGACTATCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCAc                              >  4:301609/1‑149 (MQ=255)

                                                                                                                  tATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTtggtg                                                                    <  6:112726/109‑1 (MQ=255)

                                                                                                                  tATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTtggtg                                                                    >  5:112726/1‑109 (MQ=255)

                                                                                                                       cGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCAcgcccgcc                       >  7:35412/1‑149 (MQ=255)

                                                                                                                             tCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATg                 >  2:125255/1‑149 (MQ=255)

                                                                                                                              ggtaccGACTGCGGCGTGCCGCTGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATgg                >  2:162814/6‑149 (MQ=255)

                                                                                                                               ccTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATgg                <  3:275111/148‑1 (MQ=255)

                                                                                                                               ccTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGc               >  7:101926/1‑149 (MQ=255)

                                                                                                                                cTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAg                                                    >  5:251270/1‑111 (MQ=255)

                                                                                                                                cTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAg                                                    <  6:251270/111‑1 (MQ=255)

                                                                                                                                cTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCAcgccc                          <  5:140414/137‑1 (MQ=255)

                                                                                                                                cTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCAcgccc                          >  6:140414/1‑137 (MQ=255)

                                                                                                                                cTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGcc              >  7:155262/1‑149 (MQ=255)

                                                                                                                                cTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCACACGCCCGCCGCCATGGcc              >  2:246587/1‑149 (MQ=255)

                                                                                                                                  tCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGccc             <  6:18889/148‑1 (MQ=255)

                                                                                                                                        gCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGCCCGCTGgcg      <  5:194833/149‑1 (MQ=255)

                                                                                                                                        gCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGCCCGCTGgcg      <  7:130059/149‑1 (MQ=255)

                                                                                                                                           gCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGTCTCACGACCGCCGCCATGGCCCGATGGCGCt    >  5:292367/1‑148 (MQ=255)

                                                                                                                                             aCGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCTGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGCCCGCTGGCGCTGa  >  1:90786/1‑149 (MQ=255)

CGCCCCGGCGCGCGCTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGC‑GGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGCCCGCTGGCGCTGA  >  NC_002947/2502850‑2503138

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