BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,030,782	T→A	6.3%	P360P (CCT→CCA) ‡	PP_4432 →	Xaa‑Pro aminopeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,030,782	0	T	A	6.3%	1047.2 / 42.2	316	P360P (CCT→CCA) ‡	PP_4432	Xaa‑Pro aminopeptidase
Reads supporting (aligned to +/- strand): ref base T (139/156); new base A (12/8); total (151/165)
Fisher's exact test for biased strand distribution p-value = 3.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.14e-01

Only 100 of 321 total aligned reads displayed.

GTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGG  >  NC_002947/5030633‑5030931

   cTTGAGAAGTTCGATATCCACAAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGc                                                                                                                                                    <  2:6820415/150‑1 (MQ=255)

    ttGAGAAGTTCGATATCCACAAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGTGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGcc                                                                                                                                                   >  1:3836078/1‑150 (MQ=255)

    ttGAGAAGTTCGATATCCACAAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGcc                                                                                                                                                   >  2:7332921/1‑150 (MQ=255)

    ttGAGAAGTTCGATATCCACAAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGcc                                                                                                                                                   >  2:6252423/1‑150 (MQ=255)

      gagaAGTTCGATATCCACAAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTa                                                                                                                                                 <  2:5824596/150‑1 (MQ=255)

           gTTCGATATCCACAAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCGGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTatggat                                                                                                                                            <  2:4675190/150‑1 (MQ=255)

            ttCGATATCCACAAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTatggatg                                                                                                                                           <  2:2651084/150‑1 (MQ=255)

             tCGATATCCACAAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTatggatgg                                                                                                                                          <  1:4402538/150‑1 (MQ=255)

               gATATCCACAAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATggatggat                                                                                                                                        >  1:9094246/1‑150 (MQ=255)

                 tatCCACAAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGgatga                                                                                                                                      <  1:8704839/150‑1 (MQ=255)

                   tCCACAAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTAGGGATGGATGAtt                                                                                                                                    >  1:2562057/1‑150 (MQ=255)

                    ccACAAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTAtggatgga                                                                                                                                         >  1:5714906/1‑144 (MQ=255)

                     cacaAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTgg                                                                                                                                  <  1:2953636/150‑1 (MQ=255)

                     cacaAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTgg                                                                                                                                  >  1:8204182/1‑150 (MQ=255)

                     cacaAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTgg                                                                                                                                  >  2:4399295/1‑150 (MQ=255)

                     cacaAGGACAGCCGGTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTgg                                                                                                                                  <  2:2266061/150‑1 (MQ=255)

                      acaAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTggg                                                                                                                                 <  2:2772061/150‑1 (MQ=255)

                      acaAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTggg                                                                                                                                 >  2:6096574/1‑150 (MQ=255)

                       tgaagataGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGAGTgggg                                                                                                                                >  2:3356342/8‑150 (MQ=255)

                             cAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCga                                                                                                                          <  1:8168495/150‑1 (MQ=255)

                             cAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCga                                                                                                                          >  2:8205609/1‑150 (MQ=255)

                              aGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCgag                                                                                                                         >  1:7693526/1‑150 (MQ=255)

                                  gTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTa                                                                                                                     >  2:6712704/1‑150 (MQ=255)

                                   ttGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTAc                                                                                                                    <  1:8342972/150‑1 (MQ=255)

                                    tGCGGTTACCCAATCGGTATCAGCTATCCACTAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTAcc                                                                                                                   <  2:6184044/150‑1 (MQ=255)

                                    tGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTAcc                                                                                                                   <  2:3221762/150‑1 (MQ=255)

                                      cGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCga                                                                                                                 <  1:3418729/150‑1 (MQ=255)

                                      cGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCga                                                                                                                 <  1:8548900/150‑1 (MQ=255)

                                        gTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCgaga                                                                                                               <  1:1158012/150‑1 (MQ=255)

                                         ttACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGagag                                                                                                              >  1:8933364/1‑150 (MQ=255)

                                             ccAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATc                                                                                                          <  2:2170624/150‑1 (MQ=255)

                                             ccAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATc                                                                                                          >  2:6000289/1‑150 (MQ=255)

                                             caaaTCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATc                                                                                                          >  1:8073940/3‑150 (MQ=255)

                                               aaTCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCt                                                                                                        >  2:2967253/1‑150 (MQ=255)

                                               aaTCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCt                                                                                                        >  2:6960470/1‑150 (MQ=255)

                                               aaTCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCt                                                                                                        >  2:6472472/1‑150 (MQ=255)

                                               aaTCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTAGGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCt                                                                                                        >  2:4506539/1‑150 (MQ=255)

                                               aaTCGGTATCAGCTATCCACAAGATTGGGGCGAGCGCACAATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGAACTTCCACTTCATGCCTGGCCAATGGATGGATGATTGGGGGCTCGAGATTAAAGAGAGCATCCt                                                                                                        >  2:5059265/1‑150 (MQ=255)

                                                  cGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGAt                                                                                                     <  1:2344498/150‑1 (MQ=255)

                                                    gcctgtccTATCCGCCGGACTGGGGCGAACGCACCATGAGCCTGCGTCCCGGCGACCGCACGGAATTGCAGCCCGGCATGACCTTTCATTTCATGACAGGTCTCTGGCTGGAAACCATGGGGCTCGAGATTACCGAGAGCATCGTGATCa                                                                                                   <  2:8628963/143‑1 (MQ=255)

                                                       tCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTg                                                                                                <  1:7425607/150‑1 (MQ=255)

                                                       tCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTg                                                                                                >  1:6366832/1‑150 (MQ=255)

                                                         tccTATCCGCCGGACTGGGGCGAACGCACCATGAGCCTGCGTCCCGGCGACCGCACGGAATTGCAGCCCGGCATGACCTTTCATTTCATGACAGGTCTCTGGCTGGAAACCATGGGGCTCGAGATTACCGAGAGCATCGTGATCACCGaa                                                                                              >  1:7090556/3‑149 (MQ=255)

                                                         aGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGAt                                                                                              >  1:6417252/1‑150 (MQ=255)

                                                            tATCCGCCGGACTGGGGCGAACGCACCATGAGCCTGCGTCCCGGCGACCGCACGGAATTGCAGCCCGGCATGACCTTTCATTTCATGACAGGTCTCTGGCTGGAAACCATGGGGCTCGAGATTACCGAGAGCATCGTGATCACCGAAAca                                                                                           >  1:5107315/1‑149 (MQ=255)

                                                            tATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACt                                                                                           >  2:4930929/1‑150 (MQ=255)

                                                            tATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACt                                                                                           >  1:4724692/1‑150 (MQ=255)

                                                            tATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCGATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACt                                                                                           >  2:5628785/1‑150 (MQ=255)

                                                             aTCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTg                                                                                          <  1:3959043/150‑1 (MQ=255)

                                                             aTCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTg                                                                                          >  2:750063/1‑150 (MQ=255)

                                                                 gccggactGGGGCGAACGCACCATGAGCCTGCGTCCCGGCGACCGCACGGAATTGCAGCCCGGCATGACCTTTCATTTCATGACAGGTCTCTGGCTGGAAACCATGGGGCTCGAGATTACCGAGAGCATCGTGATCACCGAAACAGGCGt                                                                                      >  1:8011057/8‑150 (MQ=255)

                                                                 accaGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGt                                                                                      >  2:805329/1‑150 (MQ=255)

                                                                 accaGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGt                                                                                      <  2:6909231/150‑1 (MQ=255)

                                                                 accaGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGt                                                                                      >  1:6265270/1‑150 (MQ=255)

                                                                   cggactGGGGCGAACGCACCATGAGCCTGCGTCCCGGCGACCGCACGGAATTGCAGCCCGGCATGACCTTTCATTTCATGACAGGTCTCTGGCTGGAAACCATGGGGCTCGAGATTACCGAGAGCATCGTGATCACCGAAACAGGCGTCg                                                                                    >  1:1566088/6‑150 (MQ=255)

                                                                   caGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCg                                                                                    >  1:6876233/1‑150 (MQ=255)

                                                                   caGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCg                                                                                    <  2:1341508/150‑1 (MQ=255)

                                                                    ggactGGGGCGAACGCACCATGAGCCTGCGTCCCGGCGACCGCACGGAATTGCAGCCCGGCATGACCTTTCATTTCATGACAGGTCTCTGGCTGGAAACCATGGGGCTCGAGATTACCGAGAGCATCGTGATCACCGAAACAGGCGTCga                                                                                   >  1:3172236/5‑150 (MQ=255)

                                                                    aGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCga                                                                                   <  1:8436342/150‑1 (MQ=255)

                                                                        tGGGGCGAACGCACCATGAGCCTGCGTCCCGGCGACCGCACGGAATTGCAGCCCGGCATGACCTTTCATTTCATGACAGGTCTCTGGCTGGAAACCATGGGGCTCGAGATTACCGAGAGCATCGTGATCACCGAAACAGGCGTCgaatgc                                                                               <  1:1656029/150‑5 (MQ=255)

                                                                        tGGGGCGAACGCACCATGAGCCTGCGTCCCGGCGACCGCACGGAATTGCAGCCCGGCATGACCTTTCATTTCATGACAGGTCTCTGGCTGGAAACCATGGGGCTCGAGATTACCGAGAGCATCGTGATCACCGAAACAGGCGTCgaatgc                                                                               >  2:5475152/1‑146 (MQ=255)

                                                                            gCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCt                                                                           >  1:7213440/1‑150 (MQ=255)

                                                                                gcgcACCCTGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTAGGACCTTCCCCTTCATGCCGGGCCGATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTc                                                                       <  2:1155778/150‑1 (MQ=255)

                                                                                gcgcACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTc                                                                       <  2:393667/150‑1 (MQ=255)

                                                                                  gcACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGa                                                                     <  1:4318263/150‑1 (MQ=255)

                                                                                    aCCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGAct                                                                   >  2:5276929/1‑150 (MQ=255)

                                                                                      cATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGActcc                                                                 >  2:7202767/1‑150 (MQ=255)

                                                                                         gAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCg                                                              <  1:8807365/150‑1 (MQ=255)

                                                                                                cgcCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCtt                                                       <  1:519280/150‑1 (MQ=255)

                                                                                                   ccGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGa                                                                     >  1:1928268/1‑133 (MQ=255)

                                                                                                   ccGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGa                                                                     <  2:1928268/133‑1 (MQ=255)

                                                                                                   ccGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTc                                                    >  2:3782811/1‑150 (MQ=255)

                                                                                                     gAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGt                                                  >  2:2375335/1‑150 (MQ=255)

                                                                                                        cGACAACAGCGTGTTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCaa                                               <  1:6703827/150‑1 (MQ=255)

                                                                                                        cGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCaa                                               <  2:4356175/150‑1 (MQ=255)

                                                                                                        cGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCaa                                               <  1:7370602/150‑1 (MQ=255)

                                                                                                        cGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTAGGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCaa                                               <  2:9081540/150‑1 (MQ=255)

                                                                                                        cGACAACAGCGTGCTAAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCaa                                               <  1:2822205/150‑1 (MQ=255)

                                                                                                            aacaGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGAc                                           >  2:2338326/1‑150 (MQ=255)

                                                                                                              caGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTg                                         >  2:5747733/1‑150 (MQ=255)

                                                                                                                  gTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGa                                     >  2:1823751/1‑150 (MQ=255)

                                                                                                                     cTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGaaaa                                  >  2:8788480/1‑150 (MQ=255)

                                                                                                                      tCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAg                                 >  1:5412238/1‑150 (MQ=255)

                                                                                                                         aaCCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGttt                              <  2:2534896/150‑1 (MQ=255)

                                                                                                                           ccTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAACTTTTCGTCAGTCCTTGACGAAAAGCTTGc                            <  1:4624887/150‑1 (MQ=255)

                                                                                                                            cTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCa                           >  1:6967024/1‑150 (MQ=255)

                                                                                                                              ggTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTTGCAAACTTTTCGTCAGTCCTTGACGAAAAGCTTGCGAg                         <  2:1092556/150‑1 (MQ=38)

                                                                                                                                   gACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCt                    >  1:3069486/1‑150 (MQ=255)

                                                                                                                                    aCCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTg                   <  2:1719442/150‑1 (MQ=255)

                                                                                                                                    aCCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTg                   <  1:957914/150‑1 (MQ=255)

                                                                                                                                      cTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGaaa                                   >  1:6671096/1‑132 (MQ=255)

                                                                                                                                      cTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGaaa                                   <  2:6671096/132‑1 (MQ=255)

                                                                                                                                      cTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGa                 <  1:3664484/150‑1 (MQ=255)

                                                                                                                                      cTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGa                 <  2:4052870/150‑1 (MQ=255)

                                                                                                                                       ttCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAg                <  1:2713963/150‑1 (MQ=255)

                                                                                                                                           aCTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTtat            >  2:7778808/1‑150 (MQ=255)

                                                                                                                                               cATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTg        <  1:7138954/150‑1 (MQ=255)

                                                                                                                                                  gCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCa     >  1:5857565/1‑150 (MQ=255)

                                                                                                                                                   ccTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCAc    >  1:5332546/1‑150 (MQ=255)

                                                                                                                                                     tGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACgg  >  1:1761254/1‑150 (MQ=255)

GTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGG  >  NC_002947/5030633‑5030931

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