BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,803,335	C→A	7.0%	P849P (CCG→CCT)	PP_3361 ←	aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,803,335	0	C	A	7.0%	611.7 / 25.1	185	P849P (CCG→CCT)	PP_3361	aminotransferase
Reads supporting (aligned to +/- strand): ref base C (90/82); new base A (2/11); total (92/93)
Fisher's exact test for biased strand distribution p-value = 1.81e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 188 total aligned reads displayed.

CCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGCAAAAACACCTGCCC  >  NC_002947/3803187‑3803484

ccAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCg                                                                                                                                                      >  1:1650618/1‑150 (MQ=255)

 cAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCgg                                                                                                                                                     <  1:608532/150‑1 (MQ=255)

 cAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCgg                                                                                                                                                     <  2:1067629/150‑1 (MQ=255)

 cAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCgg                                                                                                                                                     <  1:2480007/150‑1 (MQ=255)

 cAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCgg                                                                                                                                                     >  1:4580501/1‑150 (MQ=255)

    gCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGAc                                                                                                                                                  <  1:29577/150‑1 (MQ=255)

      cTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCa                                                                                                                                                >  1:3234732/1‑150 (MQ=255)

        aCCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACg                                                                                                                                              >  1:2518073/1‑150 (MQ=255)

             gcATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTg                                                                                                                                         <  1:952983/150‑1 (MQ=255)

              cATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCCTGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTgg                                                                                                                                        <  1:4088603/150‑1 (MQ=255)

              cATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTgg                                                                                                                                        <  1:4924214/150‑1 (MQ=255)

              cATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTgg                                                                                                                                        <  1:2677806/150‑1 (MQ=255)

                 gACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGtc                                                                                                                                     >  1:3530070/1‑150 (MQ=255)

                    acCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGTAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCTTCCGGGACCACGCCCTGGGTCTCg                                                                                                                                  >  2:3836236/1‑150 (MQ=255)

                    acCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCg                                                                                                                                  >  1:1043398/1‑150 (MQ=255)

                    acCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCg                                                                                                                                  >  2:3234878/1‑150 (MQ=255)

                    acCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCg                                                                                                                                  >  1:2212589/1‑150 (MQ=255)

                          ttgccaccGAAGGTGTTGAAATATTCCATGCCATTGTT‑AAAGCTGTCTGCCACTTCCCGCGTTGTGACAACCGCAGACATTGGATGACCGTTACCAATGGGTTTGCCCATAGTGACGATGTCAGGGACAACATCCTGCAACTCGAAAGcc                                                                                                                            <  1:3021213/143‑1 (MQ=255)

                          ttACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGcc                                                                                                                            <  1:4694406/150‑1 (MQ=255)

                           tACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGccc                                                                                                                           >  1:1457912/1‑150 (MQ=255)

                           tACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGccc                                                                                                                           <  1:2597865/150‑1 (MQ=255)

                             ccGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCccac                                                                                                                         <  2:3898604/150‑1 (MQ=255)

                             ccGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCccac                                                                                                                         <  1:2823616/150‑1 (MQ=255)

                               ggCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCccacca                                                                                                                       >  2:2445079/1‑150 (MQ=255)

                               ggCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCccacca                                                                                                                       <  1:1512929/150‑1 (MQ=255)

                               accGAAGGTGTTGAAATATTCCATGCCATTGTT‑AAAGCTGTCTGCCACTTCCCGCGTTGTGACAACCGCAGACATTGGATGACCGTTACCAATGGGTTTGCCCATAGTGACGATGTCAGGGACAACATCCTGCAACTCGAAAGCccacca                                                                                                                       >  2:796280/3‑150 (MQ=255)

                                 cGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGt                                                                                                                     <  2:162416/150‑1 (MQ=255)

                                  gAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTg                                                                                                                    >  2:4592070/1‑150 (MQ=255)

                                  gAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTg                                                                                                                    >  2:1728454/1‑150 (MQ=255)

                                   aaGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTgg                                                                                                                   <  1:82950/150‑1 (MQ=255)

                                     ggTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCt                                                                                                                 <  1:1700108/150‑1 (MQ=255)

                                        aTTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGcc                                                                                                              >  1:4455298/1‑150 (MQ=255)

                                        aTTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGcc                                                                                                              >  1:717244/1‑150 (MQ=255)

                                             aaTATTCCATGCCATTGTT‑AAAGCTGTCTGCCACTTCCCGCGTTGTGACAACCGCAGACATTGGATGACCGTTACCAATGGGTTTGCCCATAGTGACGATGTCAGGGACAACATCCTGCAACTCGAAAGCCCACCAGTGGCTGCCTAcgc                                                                                                         <  1:5017415/150‑3 (MQ=255)

                                                    cATGCCGTTGTTGAAATC‑ATCGGCAACCTCGCGGGTGTTCACTACCGCCTACATCTGGTGGGCTTTACCAATGTTTTTGCCCATGCTGTCGTCTTCCGGGTCCACTCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGat                                                                                                  <  1:1601852/150‑2 (MQ=255)

                                                    cATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGaa                                                                                                  <  1:3477964/150‑1 (MQ=255)

                                                    cATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGaa                                                                                                  >  1:1920909/1‑150 (MQ=255)

                                                    cATGCCATTGTT‑AAAGCTGTCTGCCACTTCCCGCGTTGTGACAACCGCAGACATTGGATGACCGTTACCAATGGGTTTGCCCATAGTGACGATGTCAGGGACAACATCCTGCAACTCGAAAGCCCACCAGTGGCTGCCTACGCGACCaaa                                                                                                  <  1:1505374/150‑1 (MQ=255)

                                                    cATGCCATTGTT‑AAAGCTGTCTGCCACTTCCCGCGTTGTGACAACCGCAGACATTGGATGACCGTTACCAATGGGTTTGCCCATAGTGACGATGTCAGGGACAACATCCTGCAACTCGAAAGCCCACCAGTGGCTGCCTACGCGACCaaa                                                                                                  <  1:378132/150‑1 (MQ=255)

                                                     aTGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAAc                                                                                                 <  1:4167095/150‑1 (MQ=255)

                                                     aTGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAAc                                                                                                 <  2:3723453/150‑1 (MQ=255)

                                                     aTGCCGTTGTTGAAAGC‑ATCGGCAACATCGCTCGTGGTCACTACCGCCGAAATCGGGTGGCCGTTACCAATGGGTTTGCCAATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGACGATGCCGACCCGGCCGAAc                                                                                                 >  2:1863307/1‑150 (MQ=255)

                                                         cGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGa                                                                                             >  2:3417204/1‑150 (MQ=255)

                                                         cGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGa                                                                                             <  2:1263474/150‑1 (MQ=255)

                                                         cGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGGTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGa                                                                                             >  2:2935253/1‑150 (MQ=255)

                                                              ttgAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGa                                                                                        >  1:1618627/1‑150 (MQ=255)

                                                               tgAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGaa                                                                                       <  1:3268792/150‑1 (MQ=255)

                                                                       aTCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATc                                                                                <  1:2102513/150‑1 (MQ=255)

                                                                          ggCAACCTCGCGCGTGGTCACTACCGCCGACATCTGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGc                                                                             >  2:681618/1‑150 (MQ=255)

                                                                            cAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCa                                                                           <  1:5087398/150‑1 (MQ=255)

                                                                             aaCCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAg                                                                          <  1:3926053/150‑1 (MQ=255)

                                                                                 tCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGcac                                                                      >  2:3470645/1‑150 (MQ=255)

                                                                                  cgcgcgTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGcaca                                                                     <  2:5007019/150‑1 (MQ=255)

                                                                                      cgTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGccc                                                                 <  1:4774573/150‑1 (MQ=255)

                                                                                      cgTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGccc                                                                 <  1:1190677/150‑1 (MQ=255)

                                                                                            cACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTc                                                           <  2:1763051/150‑1 (MQ=255)

                                                                                            cACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTc                                                           <  1:1274266/150‑1 (MQ=255)

                                                                                            cACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTc                                                           >  2:3830490/1‑150 (MQ=255)

                                                                                               tACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGc                                                        >  2:870080/1‑150 (MQ=255)

                                                                                               tACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGc                                                        >  1:2689989/1‑150 (MQ=255)

                                                                                                  cgccgACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCg                                                     >  1:2850176/1‑150 (MQ=255)

                                                                                                          tCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATgg                                             >  1:3373369/1‑150 (MQ=255)

                                                                                                            gggTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCg                                           >  1:3742253/1‑150 (MQ=255)

                                                                                                            gggTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCg                                           >  2:2683971/1‑150 (MQ=255)

                                                                                                             ggTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGt                                          <  1:807732/150‑1 (MQ=255)

                                                                                                               tGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTac                                        <  2:2576147/150‑1 (MQ=255)

                                                                                                                ggCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTaca                                       <  1:1248112/150‑1 (MQ=255)

                                                                                                                   cGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACAcct                                    >  1:604736/1‑150 (MQ=255)

                                                                                                                       aCCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCtt                                <  2:2207175/150‑1 (MQ=255)

                                                                                                                         cAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGa                              >  2:3667576/1‑150 (MQ=255)

                                                                                                                         cAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGa                              >  1:3071266/1‑150 (MQ=255)

                                                                                                                         cAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGa                              >  1:3112535/1‑150 (MQ=255)

                                                                                                                          aaTGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCTTGGCGTACACCTCCTTGAg                             <  2:5078205/150‑1 (MQ=255)

                                                                                                                          aaTGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAg                             >  1:973438/1‑150 (MQ=255)

                                                                                                                          aaTGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAg                             <  2:3398924/150‑1 (MQ=255)

                                                                                                                          aaTGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAg                             <  1:1748596/150‑1 (MQ=255)

                                                                                                                          aaTGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAg                             <  2:5078056/150‑1 (MQ=255)

                                                                                                                          aaTGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAg                             <  2:2784337/150‑1 (MQ=255)

                                                                                                                          aaTGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAg                             <  1:1748504/150‑1 (MQ=255)

                                                                                                                          aaTGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAg                             <  1:4576257/150‑1 (MQ=255)

                                                                                                                           aTGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCccacca                                                                                                                       >  1:4260704/1‑59 (MQ=255)

                                                                                                                           aTGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCccacca                                                                                                                       <  2:4260704/59‑1 (MQ=255)

                                                                                                                           aTGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAgg                            >  1:2139912/1‑150 (MQ=255)

                                                                                                                           aTGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAgg                            <  2:4580501/150‑1 (MQ=255)

                                                                                                                             gggTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTa                          >  2:51024/1‑150 (MQ=255)

                                                                                                                             gggTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTa                          <  1:1779560/150‑1 (MQ=255)

                                                                                                                             gggTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTa                          >  1:2641545/1‑150 (MQ=255)

                                                                                                                                    cccATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCg                   >  1:1082863/1‑150 (MQ=255)

                                                                                                                                    cccATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCg                   >  2:4505312/1‑150 (MQ=255)

                                                                                                                                     ccATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCgg                  <  2:3234732/150‑1 (MQ=255)

                                                                                                                                       aTGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGc                >  1:462161/1‑150 (MQ=255)

                                                                                                                                         gCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGCaa              >  2:3437363/1‑150 (MQ=255)

                                                                                                                                           tGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGCaaaa            <  2:1176475/150‑1 (MQ=255)

                                                                                                                                             aCGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGCAAAAac          >  1:4191374/1‑150 (MQ=255)

                                                                                                                                               ggCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGCAAAAacac        >  1:3438380/1‑150 (MQ=255)

                                                                                                                                                  aTCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGCAAAAACACCTg     >  1:4223546/1‑150 (MQ=255)

                                                                                                                                                  aTCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGACAAAAACACCTg     >  2:4025418/1‑150 (MQ=255)

                                                                                                                                                   tCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGCAAAAACACCTGc    >  2:473516/1‑150 (MQ=255)

                                                                                                                                                   tCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGCAAAAACACCTGc    >  2:473529/1‑150 (MQ=255)

                                                                                                                                                     cGGGACAACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGAACACCTCCTTGAGGTAATGTTCGGGCAAAAACACCTGccc  >  2:1776081/1‑150 (MQ=255)

CCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGCAAAAACACCTGCCC  >  NC_002947/3803187‑3803484

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