BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,288,188	C→G	7.8%	R54P (CGG→CCG)	PP_3759 ←	CheB‑like methylesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,288,188	0	C	G	7.8%	322.8 / 15.0	102	R54P (CGG→CCG)	PP_3759	CheB‑like methylesterase
Reads supporting (aligned to +/- strand): ref base C (46/48); new base G (5/3); total (51/51)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 111 total aligned reads displayed.

CGCGAGAAGTGCACGGGTTCTTCCTGGCTCAGCGACAGCGTCAGGTCGCGCTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACAATGCCGCCACGCCACCGGCCGAGGCGCCGATCACCACTGCGC  >  NC_002947/4288041‑4288335

cgcgAGAAGTGCACGGGTTCTTCCTGGCTCAGCGACAGCGTCAGGTCGCGCTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCcgccgc                                                                                                                                                    <  2:424534/150‑1 (MQ=255)

cgcgAGAAGTGCACGGGTTCTTCCTGGCTCAGCGACAGCGTCAGGTCGCGCTCCACCGACAAGTGGTAACCAACAGCCAACACATCGAAGCACCCGCGGCTGATGCTCTCCTGGTCGCGTGCCTCGCATACCGGCCGATGCAGCcgccgc                                                                                                                                                    <  1:1486631/150‑1 (MQ=255)

  cgAGAAGTGCACGGGTTCTTCCTGGCTCAGCGACAGCGTCAGGTCGCGCTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCACGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑g                                                                                                                                                 <  2:229331/150‑1 (MQ=255)

   gagaAGTGCACGGGTTCTTCCTGGCTCAGCGACAGCGTCAGGTCGCGCTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑gc                                                                                                                                                >  1:720906/1‑150 (MQ=255)

      aaGTGCACGGGTTCTTCCTGGCTCAGCGACAGCGTCAGGTCGCGCTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGACCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCGGCGGCT‑gcatc                                                                                                                                             <  2:1324498/150‑3 (MQ=255)

         tGCACGGGTTCTTCCTGGCTCAGCGACAGCGTCAGGTCGCGCTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACg                                                                                                                                          >  2:1825935/1‑150 (MQ=255)

           cACGGGTTCTTCCTGGCTCAGCGACAGCGTCAGGTCGCGCTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACgcc                                                                                                                                        <  1:575484/150‑1 (MQ=255)

              gggTTCTTCCTGGCTCAGCGACAGCGTCAGGTCGCGCTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACgccagc                                                                                                                                     >  2:1936366/1‑150 (MQ=255)

              gggTTCTTCCTGGCTCAGCGACAGCGTCAGGTCGCGCTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACgccagc                                                                                                                                     >  1:613002/1‑150 (MQ=255)

                       cTGGCTCAGCGACAGCGTCAGGTCGCGCTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCt                                                                                                                            <  1:1179103/150‑1 (MQ=255)

                         ggCTCAGCGACAGCGTCAGGTCGCGCTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGAAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCtgg                                                                                                                          >  1:1190528/1‑149 (MQ=255)

                           cTCAGCGACAGCGTCAGGTCGCGCTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTgg                                                                                                                        >  1:268447/1‑150 (MQ=255)

                                cGACAGCGTCAGGTCGCGCTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGgtc                                                                                                                   >  1:342899/1‑150 (MQ=255)

                                 gACAGCGTCAGGTCGCGCTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGgtcg                                                                                                                  >  2:1742218/1‑150 (MQ=255)

                                    aGCGTCAGGTCGCGCTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTcc                                                                                                               >  1:1455638/1‑150 (MQ=255)

                                       gTCAGGTCGCGCTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGc                                                                                                            >  1:708613/1‑150 (MQ=255)

                                        tCAGGTCGCGCTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCa                                                                                                           <  1:964607/150‑1 (MQ=255)

                                             tCGCGCTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTgc                                                                                                      <  2:1632792/150‑1 (MQ=255)

                                             tCGCGCTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTgc                                                                                                      <  2:431075/150‑1 (MQ=255)

                                                 gcTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑g                                                                                                                                                 <  1:1187587/103‑1 (MQ=255)

                                                 gcTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑g                                                                                                                                                 >  2:1187587/1‑103 (MQ=255)

                                                 gcTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTgcagca                                                                                                  >  2:1537428/1‑150 (MQ=255)

                                                  cTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCAc                                                                                                 <  1:474092/150‑1 (MQ=255)

                                                   tCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCGGCGGCT‑gca                                                                                                                                               >  1:2354095/1‑103 (MQ=255)

                                                   tCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCGGCGGCT‑gca                                                                                                                                               >  1:2354119/1‑103 (MQ=255)

                                                   tCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCGGCGGCT‑gca                                                                                                                                               <  2:2354119/103‑1 (MQ=255)

                                                   tCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACg                                                                                                >  2:1890638/1‑150 (MQ=255)

                                                      aCCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCag                                                                                             >  2:1534807/1‑150 (MQ=255)

                                                       ccGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCaga                                                                                            <  1:1532830/150‑1 (MQ=255)

                                                       ccGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCaga                                                                                            <  1:1367700/150‑1 (MQ=255)

                                                        cGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGTCTGTGGc                                                                                                                       <  1:210625/122‑1 (MQ=255)

                                                        cGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGc                                                                                                                       >  2:210625/1‑122 (MQ=255)

                                                              gTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTcc                                                                                                               <  2:2435124/124‑1 (MQ=255)

                                                              gTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTcc                                                                                                               >  1:2435124/1‑124 (MQ=255)

                                                              gTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCgg                                                                                     <  2:1361337/150‑1 (MQ=255)

                                                              gTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGACTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCgg                                                                                     >  2:1369387/1‑150 (MQ=255)

                                                               tGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGt                                                                                    >  1:1789687/1‑150 (MQ=255)

                                                                    aCCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGgtc                                                                                                                   >  1:1160637/1‑114 (MQ=255)

                                                                      ccGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGa                                                                             >  2:829223/1‑150 (MQ=255)

                                                                      ccGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGa                                                                             >  1:1141034/1‑150 (MQ=255)

                                                                      ccGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGa                                                                             >  2:496068/1‑150 (MQ=255)

                                                                            ccGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGt                                                                                                                                 <  1:1462221/92‑1 (MQ=255)

                                                                            ccGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGt                                                                                                                                 >  2:1462221/1‑92 (MQ=255)

                                                                             cGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGa                                                                      >  1:96002/1‑150 (MQ=255)

                                                                             cGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGa                                                                      >  1:1632723/1‑150 (MQ=255)

                                                                              gCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGAt                                                                     <  2:1141034/150‑1 (MQ=255)

                                                                               ccACATAGATGTGCCCTCTTCTGTTGCTCTCCTTGTCGCGTGCCACGCATACCGTCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATg                                                                    >  2:1910940/1‑150 (MQ=255)

                                                                                 acaTAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGccagccag                                                                                                                                  >  1:29950/1‑86 (MQ=255)

                                                                                 acaTAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGccagccag                                                                                                                                  <  2:29950/86‑1 (MQ=255)

                                                                                 acaTAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGc                                                                  >  2:871457/1‑150 (MQ=255)

                                                                                 acaTAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGc                                                                  <  1:1892790/150‑1 (MQ=255)

                                                                                    tAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGc                                                               <  2:342899/150‑1 (MQ=255)

                                                                                    tAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGc                                                               <  1:1537428/150‑1 (MQ=255)

                                                                                    tAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGc                                                               <  2:2120755/150‑1 (MQ=255)

                                                                                    tAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGc                                                               <  2:1789687/150‑1 (MQ=255)

                                                                                    tAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGc                                                               <  2:708613/150‑1 (MQ=255)

                                                                                       ctCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTg                                                            <  1:636952/149‑1 (MQ=255)

                                                                                         cTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGcc                                                          <  1:1890638/150‑1 (MQ=255)

                                                                                           gCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGa                                                        <  1:926015/150‑1 (MQ=255)

                                                                                               gctgctGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCaga                                                                                            >  2:1120580/1‑110 (MQ=255)

                                                                                               gctgctGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCaga                                                                                            <  1:1120580/110‑1 (MQ=255)

                                                                                                     gATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGaa                                              <  2:837357/150‑1 (MQ=255)

                                                                                                      aTGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGAGGTTGGCAGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAAc                                             >  1:1483431/1‑150 (MQ=255)

                                                                                                       tGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACa                                            >  2:399714/1‑150 (MQ=255)

                                                                                                       tGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACa                                            >  2:2125674/1‑150 (MQ=255)

                                                                                                        gCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACaa                                           >  1:2313958/1‑150 (MQ=255)

                                                                                                        gCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACaa                                           >  1:2207759/1‑150 (MQ=255)

                                                                                                         ctctGCTTGTCGCGTGCCTCGCATAGCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTAGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCCGCGTGCCGAGCACCTGGAACAAt                                          <  2:1190528/150‑1 (MQ=255)

                                                                                                         ctctCCTTGTCGCGTGCCTCGCATACCGGCCGCTGAAGCCCCCGGTAGAAGCGGGCTTG‑CAGGTGGCTGTGGCGGTCTTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACAAt                                          <  1:1742218/150‑1 (MQ=255)

                                                                                                         ctctCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACAAt                                          <  2:720906/150‑1 (MQ=255)

                                                                                                         ctctCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACAAt                                          <  2:1455638/150‑1 (MQ=255)

                                                                                                         ctctCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACAAt                                          <  1:256688/150‑1 (MQ=255)

                                                                                                            tcCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGa                                                        >  2:183309/1‑133 (MQ=255)

                                                                                                            tcCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGa                                                        <  1:183309/133‑1 (MQ=255)

                                                                                                              cTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACAAt                                          <  2:2190147/145‑1 (MQ=255)

                                                                                                              cTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACAAt                                          >  1:2190147/1‑145 (MQ=255)

                                                                                                                    gcgTGCCTCGCATACCGGCCGATGTAGCCGCCGCT‑GCAGCACGCCATCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCa                                                     <  2:1253787/128‑1 (MQ=255)

                                                                                                                    gcgTGCCTCGCATACCGGCCGATGTAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCa                                                     >  1:1253787/1‑128 (MQ=255)

                                                                                                                      gTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCgg                                                                                     >  2:2118373/1‑94 (MQ=255)

                                                                                                                      gTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCgg                                                                                     <  1:2118373/94‑1 (MQ=255)

                                                                                                                      gTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACAATGCCgccacgccac                             <  2:365543/150‑1 (MQ=255)

                                                                                                                      gTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACAATGCCgccacgccac                             <  1:829223/150‑1 (MQ=255)

                                                                                                                       tGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAACACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACAATGCCGCCACGCCAcc                            >  1:2305097/1‑150 (MQ=255)

                                                                                                                          cTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACAATGCCGCCACGCCACCGGc                         >  2:878393/1‑150 (MQ=255)

                                                                                                                           tCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACAATGCCGCCACGCCACCGGcc                        >  1:1140473/1‑150 (MQ=255)

                                                                                                                             gCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACAATGCCGCCACGCCACCGGCCGa                      <  1:1910940/150‑1 (MQ=255)

                                                                                                                                  ccGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACAATGCCGCCACGCCACCGGCCGAGgcgc                 <  1:1936366/150‑1 (MQ=255)

                                                                                                                                   cGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACAATGCCGCCACGCCACCGGCCGAGGCGcc                <  2:613002/150‑1 (MQ=255)

                                                                                                                                    ggtcGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAACACCTGGAACAATGCCGCCACGCCACCGGCCGAGGCGCCg               >  1:391987/4‑150 (MQ=255)

                                                                                                                                       tgctGCAGCGGCGGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACAATGCCGCCACGCCACCGGCCGAGGCGCCGATc            >  1:863522/4‑150 (MQ=255)

                                                                                                                                       tgctGCAGCGGCGGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACAATGCCGCCACGCCACCGGCCGAGGCGCCGATc            >  1:37790/4‑150 (MQ=255)

                                                                                                                                        gctGCAGCGGCGGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACAATGCCGCCACGCCACCGGCCGAGGCGCCGATca           >  2:2425811/3‑150 (MQ=255)

                                                                                                                                        gATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACAATGCCGCCACGCCACCGGCCGAGGCGCCGATca           >  1:960681/1‑150 (MQ=255)

                                                                                                                                           gCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGa                                                        <  2:1298494/102‑1 (MQ=255)

                                                                                                                                           gCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGa                                                        >  1:1298494/1‑102 (MQ=255)

                                                                                                                                           gCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACAATGCCGCCACGCCACCGGCCGAGGCGCCGATcacca        <  2:1632723/150‑1 (MQ=255)

                                                                                                                                             aGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGACCCATGCCGCCACGCTATAGGCCGTGCCGCTAATCTCCACt      <  2:1969532/150‑1 (MQ=255)

                                                                                                                                               ccgccgCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGc                                                                               >  2:768501/1‑75 (MQ=255)

                                                                                                                                                 gccgcT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACAATGCCGCCACGCCACCGGCCGAGGCGCCGATCACCACTgcgc  <  2:1140473/150‑1 (MQ=255)

                                                                                                                                                  ccgcT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGCACGCCGATGGCAGCGTGCCGATCACCTg                                                 >  1:491277/1‑102 (MQ=255)

CGCGAGAAGTGCACGGGTTCTTCCTGGCTCAGCGACAGCGTCAGGTCGCGCTCCACCGACAAGTGGTAACCCGGCCCCGCCACATAGATCTGCCCGCTGCTGATGCTCTCCTTGTCGCGTGCCTCGCATACCGGCCGATGCAGCCGCCGCT‑GCAGCACGCCAGCCAGTTGGCTGTGGCGGTCGTCCGGCAGGTGCAGCACGCAGAGCACCGGTATGGCGAACGCCGATGGCAGCGTGCCGAGCACCTGGAACAATGCCGCCACGCCACCGGCCGAGGCGCCGATCACCACTGCGC  >  NC_002947/4288041‑4288335

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