BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,291,468	A→C	5.9%	K85Q (AAG→CAG)	estC →	carboxylesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,291,468	0	A	C	5.9%	397.6 / 13.5	137	K85Q (AAG→CAG)	estC	carboxylesterase
Reads supporting (aligned to +/- strand): ref base A (61/67); new base C (3/5); total (65/72)
Fisher's exact test for biased strand distribution p-value = 7.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.86e-01

Only 100 of 144 total aligned reads displayed.

GGTGGGTGGCGAAACCGTCGTTGACCTGTGGGCCGGCAGTGCCGACAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTATGACTGG  >  NC_002947/1291320‑1291617

ggtgggtgGCGAAACCGTCGTTGACCTGTGGGCCGGCAGTGCCGACAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTTGTCGGCGAGGGCaa                                                                                                                                                      >  2:333635/1‑150 (MQ=255)

    ggtggCGAAACCGTCGTTGACCTGTGGGCCGGCAGTGCCGACAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTg                                                                                                                                                  <  2:2342149/150‑1 (MQ=255)

    ggtggCGAAACCGTCGTTGACCTGTGGGCCGGCAGTGCCGACAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTg                                                                                                                                                  <  1:311705/150‑1 (MQ=255)

     gtggCGAAACCGTCGTTGACCTGTGGGCCGGCAGTGCCGACAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTgg                                                                                                                                                 <  2:564045/150‑1 (MQ=255)

      tggCGAAACCGTCGTTGACCTGTGGGCCGGCAGTGCCGACAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGc                                                                                                                                                >  2:2295350/1‑150 (MQ=255)

       ggCGAAACCGTCGTTGCCCTGTGGGCCGGCAGTGCCGACAAGCACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCAGGcc                                                                                                                                               <  2:1117075/150‑1 (MQ=255)

          taaaCCGTCGTTGACCTGTGGGCCGGCAGTGCCGACAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTc                                                                                                                                            <  1:439602/149‑1 (MQ=255)

          gAAACCGTCGTTGACCTGTGGGCCGGCAGTGCCGACAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCta                                                                                                                                            >  1:1453043/1‑149 (MQ=255)

          gAAACCGTCGTTGACCTGTGGGCCGGCAGTGCCGACAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTCCTCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTc                                                                                                                                            >  1:1729748/1‑150 (MQ=255)

           aaaCCGTCGTTGACCTGTGGGCCGGCAGTGCCGACAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCCAGCTTGCCCTCg                                                                                                                                           <  1:920093/150‑1 (MQ=255)

            aaCCGTCGTTGACCTGTGGGCCGGCAGTGCCGACAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCCAGCTTGCCCTCgc                                                                                                                                          <  1:2159843/150‑2 (MQ=255)

               cgtcgtTGACCTGTGGGCCGGCAGTGCCGACAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGc                                                                                                                                       <  1:2453316/150‑1 (MQ=255)

               cgtcgtTGACCTGTGGGCCGGCAGTGCCAACAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGc                                                                                                                                       >  2:1255970/1‑150 (MQ=255)

                 tcgtTGACCTGTGGGCCGGCAGTGCCGACAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGccc                                                                                                                                     <  1:502210/150‑1 (MQ=255)

                         cTGTGGGCCGGCAGTGCCGACAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGc                                                                                                                                                <  1:998061/131‑1 (MQ=255)

                          tgtgGGCCGGCAGTGCCGACAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGACAACGACGTTGCCCTCGATGCCCCGGTGGCCa                                                                                                                            >  1:1181464/1‑150 (MQ=255)

                                  ggCAGTGCCGACAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTgg                                                                                                                    <  2:2094571/150‑1 (MQ=255)

                                  ggCAGTGCCGACAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTgg                                                                                                                    <  2:116106/150‑1 (MQ=255)

                                  ggCAGTGCCGACAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTgg                                                                                                                    <  1:1630866/150‑1 (MQ=255)

                                     aGTGCCGACAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGccc                                                                                                                 <  1:2295350/150‑1 (MQ=255)

                                      gTGCCGACAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCTAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCg                                                                                                                >  2:1699582/1‑150 (MQ=255)

                                       tGCCGACAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGa                                                                                                               >  1:584551/1‑150 (MQ=255)

                                           gACAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTc                                                                                                           >  1:633775/1‑150 (MQ=255)

                                             cAAGGACGGCCGGCAGGCCGGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGc                                                                                                         <  2:1197082/150‑1 (MQ=255)

                                                ggACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCa                                                                                                      >  1:478442/1‑150 (MQ=255)

                                                ggACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCa                                                                                                      >  2:202767/1‑150 (MQ=255)

                                                  aCGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAgg                                                                                                    <  2:584551/150‑1 (MQ=255)

                                                   cGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGc                                                                                                   <  1:2082905/150‑1 (MQ=255)

                                                    ggCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACATGCCCGAGTTCGCCCAGGCg                                                                                                  >  1:722777/1‑150 (MQ=255)

                                                     gCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCTAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCgg                                                                                                 <  1:1699582/150‑1 (MQ=255)

                                                     gCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCgg                                                                                                 <  2:1162451/150‑1 (MQ=255)

                                                     gCCAGCACCCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCgg                                                                                                 <  1:702007/150‑1 (MQ=255)

                                                      ccAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGccc                                                                                                                 >  1:1855461/1‑133 (MQ=255)

                                                      ccAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGccc                                                                                                                 <  2:1855461/133‑1 (MQ=255)

                                                         gcagGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTaa                                                                                             >  1:2372865/1‑150 (MQ=255)

                                                            ggCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCa                                                                                          <  1:1255970/150‑1 (MQ=255)

                                                            ggCCTGGAACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCa                                                                                          >  2:832304/1‑150 (MQ=255)

                                                            ggCCTCGTACAGCGCTACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCAGGTCCTCGCCC‑GTAGGTACGACGCGAGGGCAAGCTGGCCATCGATGCCCATGTGCTCACCTAGTTGCCCGAGTTCGCCCAGGCGGGTAAGCa                                                                                          <  1:2387854/150‑1 (MQ=255)

                                                              ccTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTg                                                                                                                     >  2:1013250/1‑121 (MQ=255)

                                                              ccTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTg                                                                                                                     <  1:1013250/121‑1 (MQ=255)

                                                                 ggCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGa                                                                                     >  2:2277641/1‑150 (MQ=255)

                                                                 ggCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGa                                                                                     <  1:1136621/150‑1 (MQ=255)

                                                                 ggCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGa                                                                                     >  2:1022009/1‑150 (MQ=255)

                                                                    acaGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCa                                                                                  >  1:1821187/1‑150 (MQ=255)

                                                                        cGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCtt                                                                              <  2:478442/150‑1 (MQ=255)

                                                                         gATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTg                                                                             >  2:2099099/1‑150 (MQ=255)

                                                                            aCCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGt                                                                          <  1:2183027/150‑1 (MQ=255)

                                                                             ccATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGc                                                                            <  1:1315923/147‑1 (MQ=255)

                                                                              cATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTc                                                                                                           <  1:2242086/115‑1 (MQ=255)

                                                                              cATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTaa                                                                                             <  2:366898/129‑1 (MQ=255)

                                                                              cATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCa                                                                        >  1:202519/1‑150 (MQ=255)

                                                                              cATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCACTGCGCCCCTGCAAGCGCTGGTTCGCCCAGGCGGTGAAGCAGCCGATCACCTTGCGTcg                                                                        >  2:62688/1‑149 (MQ=255)

                                                                                    cAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGa                                                                                                               >  2:738218/1‑105 (MQ=255)

                                                                                     aaCCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTc                                                                 <  1:1775787/150‑1 (MQ=255)

                                                                                           ttCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGAt                                                                                    <  1:1587505/125‑1 (MQ=255)

                                                                                           ttCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGAt                                                                                    >  2:1587505/1‑125 (MQ=255)

                                                                                           ttCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCAc                                                           >  2:221646/1‑150 (MQ=255)

                                                                                            tctcCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCAcc                                                          >  2:2377169/1‑150 (MQ=255)

                                                                                                cTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCa                                                                                          <  2:758942/114‑1 (MQ=255)

                                                                                                cTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCa                                                                                          >  1:758942/1‑114 (MQ=255)

                                                                                                cTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGc                                                      >  1:487232/1‑150 (MQ=255)

                                                                                                cTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGc                                                      >  2:2019420/1‑150 (MQ=255)

                                                                                                cTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGATGGCCCTCTATGCACCGGTGGCCAACTACGGACGCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTAAACTGCACAGACAACGTGc                                                      >  2:1092721/1‑150 (MQ=255)

                                                                                                     ccAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCt                                                 >  2:1700521/1‑150 (MQ=255)

                                                                                                      cAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCa                                                                        >  2:143251/1‑126 (MQ=255)

                                                                                                      cAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCtt                                                <  2:452830/150‑1 (MQ=255)

                                                                                                       aaGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTg                                               >  2:2444312/1‑150 (MQ=255)

                                                                                                       aaGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTg                                               >  1:370692/1‑150 (MQ=255)

                                                                                                        aGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGc                                              <  1:2019420/150‑1 (MQ=255)

                                                                                                        aGACATTCGCGGCCGTTAC‑CGGCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGc                                              >  2:2368589/1‑150 (MQ=255)

                                                                                                           ccTTCGCGGCCGTTAC‑CGCCCTTCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCgg                                           >  1:1044738/1‑150 (MQ=255)

                                                                                                             ttCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCAAGCCGGca                                         >  1:1450074/1‑149 (MQ=255)

                                                                                                                gcgGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGAAAGCTGGCCATCGATGAACCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCAtt                                      <  2:1707466/150‑1 (MQ=255)

                                                                                                                  ggCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCg                                    >  2:294986/1‑150 (MQ=255)

                                                                                                                   gCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGt                                   <  2:487232/150‑1 (MQ=255)

                                                                                                                     cGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGa                                 <  1:786176/150‑1 (MQ=255)

                                                                                                                     cGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGa                                 <  2:633775/150‑1 (MQ=255)

                                                                                                                          c‑cGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAgctgc                            <  1:2327622/150‑1 (MQ=255)

                                                                                                                             gCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGctgctg                          >  2:2399621/1‑150 (MQ=255)

                                                                                                                               ccTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGcc                        <  2:1769156/150‑1 (MQ=255)

                                                                                                                                 tGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCgg                      <  1:294986/150‑1 (MQ=255)

                                                                                                                                    aGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTCAGCAATCTATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTg                   >  2:566769/1‑150 (MQ=255)

                                                                                                                                      ttGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGaa                 <  1:1022009/150‑1 (MQ=255)

                                                                                                                                         gTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCtgctgc                         <  1:2034535/139‑1 (MQ=255)

                                                                                                                                         gTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGcc              >  1:1327656/1‑150 (MQ=255)

                                                                                                                                           cGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCtt            <  1:333635/150‑1 (MQ=255)

                                                                                                                                            ggCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTg           <  1:2433447/150‑1 (MQ=255)

                                                                                                                                            atcGAGGGCCAGCTTGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTc                                                                         >  2:813567/3‑88 (MQ=255)

                                                                                                                                            atcGAGGGCCAGCTTGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTc                                                                         <  1:813567/86‑1 (MQ=255)

                                                                                                                                             tcGAGGGCCAGCTTGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGc               >  2:1747069/2‑145 (MQ=255)

                                                                                                                                             tcGAGGGCCAGCTTGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGc               <  1:1747069/144‑1 (MQ=255)

                                                                                                                                             gCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGt          <  1:19681/150‑1 (MQ=255)

                                                                                                                                                 gggCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTATGa      >  2:1167158/1‑150 (MQ=255)

                                                                                                                                                    cAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGaa                 >  1:1938079/1‑136 (MQ=255)

                                                                                                                                                    cAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGaa                 <  2:1938079/136‑1 (MQ=255)

                                                                                                                                                     aaGCTGGCCCTCGCTGCCCCCGTGCCCGACCAGTGGCCCGAGTTGCCCCAGCCGGGTAAGCATTCGATCACCTTTCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTATGACTgg  <  2:146335/150‑1 (MQ=255)

                                                                                                                                                     aaGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTATGACTgg  >  1:1867556/1‑150 (MQ=255)

                                                                                                                                                     aaGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTATGACTgg  <  1:1167158/150‑1 (MQ=255)

                                                                                                                                                     aaGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCAAGCCGGcac                                        >  2:1374754/1‑110 (MQ=255)

                                                                                                                                                     aaGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCAAGCCGGcac                                        <  1:1374754/112‑3 (MQ=255)

GGTGGGTGGCGAAACCGTCGTTGACCTGTGGGCCGGCAGTGCCGACAAGGACGGCCAGCAGGCCTGGCACAGCGATACCATCGCCAACCTGTTCTCCTGCACCAAGACCTTCGCGGCCGTTAC‑CGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTATGACTGG  >  NC_002947/1291320‑1291617

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