BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,926,936	A→G	5.4%	T54A (ACT→GCT)	PP_2565 →	UvrD/REP family helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,926,936	0	A	G	5.4%	332.9 / 5.6	111	T54A (ACT→GCT)	PP_2565	UvrD/REP family helicase
Reads supporting (aligned to +/- strand): ref base A (55/50); new base G (2/4); total (57/54)
Fisher's exact test for biased strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

Only 100 of 111 total aligned reads displayed.

TGGCCCCCGCCGACCTTACCCCTCCTGCCCAGTTGCCCTGGTTTCGCCGCCTGGCCGCCCGCCTGCTGGGGCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGGCGCCAGGTGCTGGTGATCCGCGATACGCGCCCCGACAGCGCCGCTGTGCCAGGCCAGGATGACAACTTGTTCGACGAC  >  NC_002947/2926796‑2927079

tGGCCCCCGCCGACCTTACCCCTCCTGCCCAGTTGCCCTGGTTTCGCCGCCTGGCCGCCCGCCTGCTGGGGCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGTGCGCTGGcctgt                                                                                                                                        <  1:955040/150‑4 (MQ=255)

tGGCCCCCGCCGACCTTACCCCTCCTGCCCAGTTGCCCTGGTTTCGCCGCCTGGCCGCCCGCCTGCTGGGGCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGTGCGCTGGcctgt                                                                                                                                        <  1:432617/150‑4 (MQ=255)

 ggCCCCCGCCGCCCTTACCCCTCCAGCCCAATGGCACTCGGTGGGCCTCCTGGCCGCCCGCCTGCTGGGGCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGc                                                                                                                                       <  2:857767/150‑1 (MQ=255)

 ggCCCCCGCCGACCTTACCCCTCCTGCCCAGTTGCCCTGGTTTCGCCGCCTGGCCGCCCGCCTGCTGGGGCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGTGCGCTGGcctgtg                                                                                                                                       <  1:2177287/150‑5 (MQ=255)

   cccccGCCGACCTTACCCCTCCTGCCCAGTTGCCCTGGTTTCGCCGCCTGGCCGCCCGCCTGCTGGGGCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCg                                                                                                                                     >  2:2292038/1‑150 (MQ=255)

    ccccGCCGACCTTACCCCTCCTGCCCAGTTGCCCTGGTTTCGCCGCCTAGCCGCCCGCCTGCTGGGGCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGAACTGGCCATGCCGa                                                                                                                                    >  2:1752867/1‑150 (MQ=255)

       cgccgACCTTACCCCTCCTGCCCAGTTGCCCTGGTTTCGCCGCCTGGCCGCCCGCCTGCTGGGGCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATgg                                                                                                                                 <  2:2281456/150‑1 (MQ=255)

        gccgACCTTACCCCTCCTGCCCAGTTGCCCTGGTTTCGCCGCCTGGCCGCCCGCCTGCTGGGGCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGgc                                                                                                                                <  1:357921/150‑1 (MQ=255)

              cTTACCCCTCCTGCCCAGTTGCCCTGGTTTCGCCGCCTGGCCGCCCGCCTGCTGGGGCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGgcgtgcgt                                                                                                                          <  2:1820371/150‑1 (MQ=255)

               ttACCCCTCCTGCCCAGTTGCCCTGGTTTCGCCGCCTGGCCGCCCGCCTGCTGGGGCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGcgtgcgtg                                                                                                                         >  2:586571/1‑150 (MQ=255)

                 aCCCCTCCTGCCCAGTTGCCCTGGTTTCGCCGCCTGGCCGCCCGCCTGCTGGGGCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGaa                                                                                                                       <  1:2161188/150‑1 (MQ=255)

                 aCCCCTCCTGCCCAGTTGCCCTGGTTTCGCCGCCTGGCCGCCCGCCTGCTGGGGCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGaa                                                                                                                       >  1:1188996/1‑150 (MQ=255)

                  ccccTCCTGCCCAGTTGCCCTGGTTTCGCCGCCTGGCCGCCCGCCTGCTGGGGCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGc                                                                                                                                       <  2:892323/133‑1 (MQ=255)

                  ccccTCCTGCCCAGTTGCCCTGGTTTCGCCGCCTGGCCGCCCGCCTGCTGGGGCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGc                                                                                                                                       >  1:892323/1‑133 (MQ=255)

                   cccTCCTGCCCAGTTGCCCTGGTTTCGCCGCCTGGCCGCCCGCCTGCTGGGGCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGTGCGCTGGcctgtggc                                                                                                                                     >  1:99624/1‑128 (MQ=255)

                   cccTCCTGCCCAGTTGCCCTGGTTTCGCCGCCTGGCCGCCCGCCTGCTGGGGCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAgg                                                                                                                     >  1:428904/1‑150 (MQ=255)

                     ctcctGCGCAGTTGCCCTGGTTTCGCCGCCTGGCCGCCCGCCTGCTGGGGCGTGGCCTGACCCGCCGGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAgggg                                                                                                                   >  1:790380/1‑149 (MQ=255)

                     ctcctGCCCAGTTGCCCTGGTTTCGCCGCCTGGCCGCCCGCCTGCTGGGGCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTGCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGt                                                                                                                   <  1:2172562/150‑1 (MQ=255)

                      tcctGCCCAGTTGCCCTGGTTTCGCCGCCTGGCCGCCCGCCTGCTGGGGCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGtt                                                                                                                  >  2:984335/1‑150 (MQ=255)

                         tGCCCAGTTGCCCTGGTTTCGCCGCCTGGCCGCCCGCCTGCTGGGGCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCAACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGa                                                                                                               <  1:831982/150‑1 (MQ=255)

                             cAGTTGCCCTGGTTTCGCCGCCTGGCCGCCCGCCTGCTGGGGCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGt                                                                                                           <  1:660975/150‑1 (MQ=255)

                              aGTTGCCCTGGTTTCGCCGCCTGGCCGCCCGCCTGCTGGGGCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTg                                                                                                          <  2:1532584/150‑1 (MQ=255)

                                   cccTGGTTTCGCCGCCTGGCCGCCCGCCTGCTGGGGCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCgg                                                                                                     >  1:1063371/1‑150 (MQ=255)

                                   cccTGGTTTCGCCGCCTGGCCGCCCGCCTGCTGGGGCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCgg                                                                                                     >  2:1617510/1‑150 (MQ=255)

                                   cccTGGTTTCGCCGCCTGGCCGCCCGCCTGCTGGGGCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCgg                                                                                                     >  2:1873547/1‑150 (MQ=255)

                                    ccTGGTTTCGCCGCCTGGCCGCCCGCCTGCTGGGGCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCggg                                                                                                    <  2:2247280/150‑1 (MQ=255)

                                    ccTGGTTTCGCCGCCTGGCCGCCCGCCTGCTGGGGCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCggg                                                                                                    >  2:141725/1‑150 (MQ=255)

                                     cTGGTTTCGCCGCCTGGCCGCCCGCCTGCTGGGGCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGGGTGAAGGGTTCGAGCGTGGTCGGGt                                                                                                   >  2:885414/1‑150 (MQ=255)

                                     cTGGTTTCGCCGCCTGGCCGCCCGCCTGCTGGGGCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGt                                                                                                   >  2:655562/1‑150 (MQ=255)

                                     cTGGTTTCGCCGCCTGGCCGCCCGCCTGCTGGGGCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGt                                                                                                   <  1:1873547/150‑1 (MQ=255)

                                     cTGGTTTCGCCGCCTGGCCGCCCGCCTGCTGGGGCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGt                                                                                                   >  2:2218950/1‑150 (MQ=255)

                                         tttCGCCGCCTGGCCGCCCGCCTGCTGGGGCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAg                                                                                               <  1:1634077/150‑1 (MQ=255)

                                         tttCGCCGCCTGGCCGCCCGCCTGCTGGGGCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAg                                                                                               <  1:1823393/150‑1 (MQ=255)

                                          ttCGCCGCCTGGCCGCCCGCCTGCTGGGGCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAgg                                                                                              >  2:1554997/1‑150 (MQ=255)

                                             gccgccTGGCCGCCCGCCTGCTGGGGCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGcc                                                                                           >  2:2112543/1‑149 (MQ=255)

                                                    ggCCGCCCGCCTGCTGGGGCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGc                                                                                    <  1:984335/150‑1 (MQ=255)

                                                       cgcccgccTGCTGGGGCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCggg                                                                                 >  1:1677974/1‑150 (MQ=255)

                                                       cgcccgccTGCTGGGGCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCAGGCCGATGGCGTGCGTGAAGGGTTAGAGCGTGGTCGAGTGGAGGGCTTTGAAGCggg                                                                                 <  1:655562/150‑1 (MQ=255)

                                                         cccgccTGCTGGGGCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGgc                                                                               >  1:1093026/1‑150 (MQ=255)

                                                           cgccTGCTGGGGCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGgcgc                                                                             >  1:2008423/1‑150 (MQ=255)

                                                           cgccTGCTGGGGCGTGGCCTCACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGgcgc                                                                             >  2:1402089/1‑150 (MQ=255)

                                                            gccTGCTGGGGCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGGCGcc                                                                            >  1:1167560/1‑150 (MQ=255)

                                                             ccTGCTGGGGCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCgg                                                                                  >  1:91692/1‑143 (MQ=255)

                                                             ccTGCTGGGGCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCgg                                                                                  <  2:91692/143‑1 (MQ=255)

                                                             ccTGCTGGGGCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGGCGCCa                                                                           <  1:2441018/150‑1 (MQ=255)

                                                             ccTGCTGGGGCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGGCGCCa                                                                           >  2:848005/1‑150 (MQ=255)

                                                                 ctgGGGCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGGCGCCAGGTg                                                                       >  1:1376450/1‑150 (MQ=255)

                                                                     ggCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGGCGCCAGGTGCTgg                                                                   <  1:1326815/150‑1 (MQ=255)

                                                                     ggCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGGCGCCAGGTGCTgg                                                                   <  1:2292038/150‑1 (MQ=255)

                                                                          ggCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGGCGCCAGGTGCTGGTGATc                                                              <  1:944697/150‑1 (MQ=255)

                                                                           gCCTGACCCGCCTGCAAGCTGAACACCGCGACTCATGGTTCCTCGGCCACGCCACCGGCAACCGCACTGGCAATGCCGATGGCGTGCTTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGGCGCCAGGTGCTGGTGATcc                                                             <  1:1945481/150‑1 (MQ=255)

                                                                             cTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCggg                                                                                 >  1:589952/1‑128 (MQ=255)

                                                                             cTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCggg                                                                                 <  2:589952/128‑1 (MQ=255)

                                                                              tGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGGCGCCAGGTGCTGGTGATCcgcg                                                          >  2:1525187/1‑150 (MQ=255)

                                                                               gACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGGCGCCAGGTGCTGGTGATCCGCGa                                                         <  1:2218950/150‑1 (MQ=255)

                                                                                  ccGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCg                                                                                   <  1:2325229/121‑1 (MQ=255)

                                                                                  ccGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCg                                                                                   >  2:2325229/1‑121 (MQ=255)

                                                                                  ccGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGGCGCCAGGTGCTGGTGATCCGCGATAc                                                      >  1:2012525/1‑150 (MQ=255)

                                                                                  ccGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGGCGCCAGGTGCTGGTGATCCGCGATAc                                                      >  1:2279137/1‑150 (MQ=255)

                                                                                  ccGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGGCGCCAGGTGCTGGTGATCCGCGATAc                                                      >  1:1489315/1‑150 (MQ=255)

                                                                                     ccTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGc                                                                                                                                       <  1:1951564/66‑1 (MQ=255)

                                                                                     ccTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGc                                                                                                                                       >  2:1951564/1‑66 (MQ=255)

                                                                                      cggCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCCAAGGCCGGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGGCGCCAGGTGCTGGTGATCCGCGATACgcgc                                                  <  2:790380/148‑1 (MQ=255)

                                                                                      cTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGGCGCCAGGTGCTGGTGATCCGCGATACgcgc                                                  <  2:2364834/150‑1 (MQ=255)

                                                                                       tGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGGCGCCAGGTGCTGGTGATCCGCGATACGCGcc                                                 >  2:1672964/1‑150 (MQ=255)

                                                                                        gCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGGCGCCAGGTGCTGGTGATCCGCGATACGCGccc                                                >  2:2416581/1‑150 (MQ=255)

                                                                                        gCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAACACGACGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGGCACAAGGTGCTGGTGATCCGCGATACGCGccc                                                >  1:148731/1‑150 (MQ=255)

                                                                                         cAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGGCGCCAGGTGCTGGTGATCCGCGATACGCGcccc                                               <  1:133085/150‑1 (MQ=255)

                                                                                         cAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGGCGCCAGGTGCTGGTGATCCGCGATACGCGcccc                                               <  2:1776372/150‑1 (MQ=255)

                                                                                          aaGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGGCGCCAGGTGCTGGTGATCCGCGATACGCGCCCCg                                              >  1:159454/1‑150 (MQ=255)

                                                                                            gCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGGCGCCAGGTGCTGGTGATCCGCGATACGCGCCCCGAc                                            >  1:939068/1‑150 (MQ=255)

                                                                                             cccAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGGCGCCAGGTGCTGGTGATCCGCGATACGCGCCCCGACa                                           <  2:939068/150‑1 (MQ=255)

                                                                                              ccAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGCAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGGCGCCAGGTGCTGGTGATCCGCGATAGGGGCCCCGAccg                                          >  2:974630/1‑148 (MQ=255)

                                                                                               cAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGa                                                                                                               <  2:2114833/80‑1 (MQ=255)

                                                                                               cAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGa                                                                                                               >  1:2114833/1‑80 (MQ=255)

                                                                                               cAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCgg                                                                                  <  1:968990/109‑1 (MQ=255)

                                                                                               cAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCgg                                                                                  >  2:968990/1‑109 (MQ=255)

                                                                                                   aCCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGGCGCCAGGTGCTGGTGATCCGCGATACGCGCCCCGACAGcgccg                                     >  1:758251/1‑150 (MQ=255)

                                                                                                          cTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGGCGCCAGGTGCTGGTGATCCGCGATACGCGCCCCGACAGCGCCGCTGTgcc                              >  2:645166/1‑150 (MQ=255)

                                                                                                            cATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCTGGTGGAGGGCTTTGAAGCGGGGCGCCAGGTGCTGGTGATCCGCGATACGCGCCCCGACAGCGCCGCTGTgccag                            <  1:586571/150‑1 (MQ=255)

                                                                                                            cATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGGCGCCAGGTGCTGGTGATCCGCGATACGCGCCCCGACAGCGCCGCTGTgccag                            <  2:1815032/150‑1 (MQ=255)

                                                                                                            cATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGGCGCCAGGTGCTGGTGATCCGCGATACGCGCCCCGACAGCGCCGCTGTgccag                            <  2:428904/150‑1 (MQ=255)

                                                                                                             aTGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGGCGCCAGGTGCTGGTGATCCGCGATACGCGCCCCGACAGCGCCGCTGTgccagg                           <  1:1525187/150‑1 (MQ=255)

                                                                                                                 ttCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGGCGCCAGGTGCTGGTGATCCGCGATACGCGCCCCGACAGCGCCGCTGTgccaggccag                       <  2:159454/150‑1 (MQ=255)

                                                                                                                   ccTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGGCGCCAGGTGCTGGTGATCCGCGATACGCGCCCCGACAGCGCCGCTGTGCCAGGCCAGGa                     <  2:758251/150‑1 (MQ=255)

                                                                                                                         ccacgccacAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGGCGCCAGGTGCTGGTGATCCGCGATACGCGcccc                                               >  1:14892/1‑118 (MQ=255)

                                                                                                                         ccacgccacAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGGCGCCAGGTGCTGGTGATCCGCGATACGCGcccc                                               <  2:14892/118‑1 (MQ=255)

                                                                                                                          cacgccacAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGGCGCCAGGTGCTGGTGATCCGCGATACGCGCCCCGACAGCGCCGCTGTGCCAGGCCAGGATGACAAc              <  1:1442034/150‑1 (MQ=255)

                                                                                                                           acgccacAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGGCGCCAGGTGCTGGTGATCCGCGATACGCGCCCCGACAGCGCCGCTGTGCCAGGCCAGGATGACAACt             >  2:1650745/1‑150 (MQ=255)

                                                                                                                             gccacAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGGCGCCAGGTGCTGGTGATCCGCGATACGCGCCCCGACAGCGCCGCTGTGCCAGGCCAGGATGACAACTTg           <  1:90300/150‑1 (MQ=255)

                                                                                                                               cacaGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGGCGGGTGGAGGGCTTTGAAGCGGGGCGCCAGGTGCTGGTGATCCGCGATACGCGCCCCGACAGCGCCGCTGTGCCAGGCCAGGATGACAACGTGtt         >  1:278550/1‑150 (MQ=255)

                                                                                                                                 atggCCAGTGCGCTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGGCGCCAGGTGCTGGTGATCCGCGATACGCGCCCCGACAGCGCCGCTGTGCCAGGCCAGGATGACAACTTGTTcg       >  2:1122235/3‑150 (MQ=255)

                                                                                                                                   ggCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGGCGCCAGGTGCTGGTGATCCGCGATACGCGCCCCGACAGCGCCGCTGTGCCAGGCCAGGATGACAACCTGTTcgac     >  2:1629524/1‑150 (MQ=255)

                                                                                                                                      cAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGGCGCCAGGTGCTGGTGATCCGCGa                                                         >  2:1303734/1‑95 (MQ=255)

                                                                                                                                      cAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGGCGCCAGGTGCTGGTGATCCGCGa                                                         <  1:1303734/95‑1 (MQ=255)

                                                                                                                                      cAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGGCGCCAGGTGCTGGTGATCCGCGATACGCGCCCCGACAGCGCCGCTGTGCCAGGCCAGGATGACAAc              <  2:825504/138‑1 (MQ=255)

                                                                                                                                      cAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGGCGCCAGGTGCTGGTGATCCGCGATACGCGCCCCGACAGCGCCGCTGTGCCAGGCCAGGATGACAAc              >  1:825504/1‑138 (MQ=255)

                                                                                                                                      cAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGGCGCCAGGTGCTGGTGATCCGCGATACGCGCCCCGACAGCGCCGCTGTGCCAGGCCAGGATGACAACTTGTTCgacgac  <  1:1402089/150‑1 (MQ=255)

                                                                                                                                        gcgcACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGGCGCCAGGt                                                                        <  2:457994/78‑1 (MQ=255)

                                                                                                                                        gcgcACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGGCGCCAGGt                                                                        >  1:457994/1‑78 (MQ=255)

TGGCCCCCGCCGACCTTACCCCTCCTGCCCAGTTGCCCTGGTTTCGCCGCCTGGCCGCCCGCCTGCTGGGGCGTGGCCTGACCCGCCTGCAAGCCCAGCACCGCGACTCATGGTTCCTCGGCCACGCCACAGGCCAGCGCACTGGCCATGCCGATGGCGTGCGTGAAGGGTTCGAGCGTGGTCGGGTGGAGGGCTTTGAAGCGGGGCGCCAGGTGCTGGTGATCCGCGATACGCGCCCCGACAGCGCCGCTGTGCCAGGCCAGGATGACAACTTGTTCGACGAC  >  NC_002947/2926796‑2927079

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: `‑`