BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,030,798	G→A	6.9%	D366N (GAT→AAT) ‡	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,798	0	G	A	6.9%	898.6 / 37.2	291	D366N (GAT→AAT) ‡	PP_4432	Xaa‑Pro aminopeptidase
Reads supporting (aligned to +/- strand): ref base G (135/136); new base A (12/8); total (147/144)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 311 total aligned reads displayed.

ATATCCACAAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACA  >  NC_002947/5030649‑5030946

atatCCACAAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGgatggatg                                                                                                                                                      >  1:8295894/1‑150 (MQ=255)

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

   tCCACAAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGAtt                                                                                                                                                   <  2:7190787/150‑1 (MQ=255)

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

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

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

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

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

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

       caAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTgggg                                                                                                                                               >  2:1081991/1‑150 (MQ=255)

           gACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTc                                                                                                                                           >  2:7923516/1‑150 (MQ=255)

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

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

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

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

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

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

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

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

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

                              cAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATcc                                                                                                                        <  1:3287875/150‑1 (MQ=255)

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

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

                                   ggTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATc                                                                                                                   >  1:7671189/1‑150 (MQ=255)

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

                                     cctgtccTATCCGCCGGACTGGGGCGAACGCACCATGAGCCTGCGTCCCGGCGACCGCACGGAATTGCAGCCCGGCATGACCTTTCATTTCATGACAGGTCTCTGGCTGGAAACCATGGGGCTCGAGATTACCGAGAGCATCGTGATCAc                                                                                                                 >  1:7004957/7‑150 (MQ=255)

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

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

                                          gCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATa                                                                                                            >  1:3693395/1‑150 (MQ=255)

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

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

                                              tCCGCCGGACTGGGGCGAACGCACCATGAGCCTGCGTCCCGGCGACCGCACGGAATTGCAGCCCGGCATGACCTTTCATTTCATGACAGGTCTCTGGCTGGAAACCATGGGGCTCGAGATTACCGAGAGCATCGTGATCACCGAAACAgg                                                                                                        >  1:6695905/1‑150 (MQ=255)

                                              tCCGCCGGACTGGGGCGAACGCACCATGAGCCTGCGTCCCGGCGACCGCACGGAATTGCAGCCCGGCATGACCTTTCATTTCATGACAGGTCTCTGGCTGGAAACCATGGGGCTCGAGATTACCGAGAGCATCGTGATCACCGAAACAgg                                                                                                        >  1:2537756/1‑150 (MQ=255)

                                               ccaccaGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTggg                                                                                                       <  1:9126539/150‑1 (MQ=255)

                                               ccaccaGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTggg                                                                                                       <  1:6329305/150‑1 (MQ=255)

                                               ccaccaGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTggg                                                                                                       <  2:1284564/150‑1 (MQ=255)

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

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

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

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

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

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

                                                            gCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCCTGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCt                                                                                          <  2:4058509/150‑1 (MQ=255)

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

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

                                                                   cACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGAc                                                                                   <  2:7332081/150‑1 (MQ=255)

                                                                   cACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGAc                                                                                   <  2:1142458/150‑1 (MQ=255)

                                                                      cATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCGAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGActcc                                                                                >  2:4676379/1‑150 (MQ=255)

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

                                                                            ccTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCaa                                                                          <  1:7383089/150‑1 (MQ=255)

                                                                             cTGCGCCCGAGCGACAACAGCGTGCTAAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAg                                                                         >  2:6943189/1‑150 (MQ=255)

                                                                              tGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGc                                                                        >  1:2704373/1‑150 (MQ=255)

                                                                              tGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGc                                                                        <  2:2636366/150‑1 (MQ=255)

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

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

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

                                                                                   ccGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTc                                                                   <  2:1710194/150‑1 (MQ=255)

                                                                                   ccGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTc                                                                   <  2:6057332/150‑1 (MQ=255)

                                                                                   ccGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTc                                                                   <  2:5094827/150‑1 (MQ=255)

                                                                                   ccGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTc                                                                   <  2:6876233/150‑1 (MQ=255)

                                                                                   ccGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTc                                                                   <  2:2352279/150‑1 (MQ=255)

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

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

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

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

                                                                                               aGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGa                                                       >  2:5860046/1‑150 (MQ=255)

                                                                                                  gTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGa                                                    <  1:4211957/150‑1 (MQ=255)

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

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

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

                                                                                                        aaaCCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGtt                                              <  2:5004261/150‑1 (MQ=255)

                                                                                                        aaaCCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGtt                                              >  1:574583/1‑150 (MQ=255)

                                                                                                        aaaCCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGtt                                              >  1:6899836/1‑150 (MQ=255)

                                                                                                         aaCCTGGTATGACCTTCCACTTCCTGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGNCTCCTCGCAAACTTTTCGTCAGTCCTTGACGAAAAGctt                                             <  2:6366832/150‑1 (MQ=255)

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

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

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

                                                                                                           ccTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGc                                           >  2:177869/1‑150 (MQ=255)

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

                                                                                                               gTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCAATGATACTGGGGTCTAGACGATCTGCTCGACTCAACGCAAGAATTTCGTCAAGGACTGAAGAAAAGTTAGAAAGt                                       >  1:1966295/1‑150 (MQ=255)

                                                                                                                 aTGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGTTACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCg                                     >  2:6556738/1‑150 (MQ=255)

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

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

                                                                                                                      cTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGa                                >  1:6754432/1‑150 (MQ=255)

                                                                                                                      cTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGa                                >  1:6602247/1‑150 (MQ=255)

                                                                                                                         ccACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGtt                             >  2:232283/1‑150 (MQ=255)

                                                                                                                            cTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTtata                          <  1:5885918/150‑1 (MQ=255)

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

                                                                                                                               cATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTC‑‑CAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGcc                     <  1:2967253/150‑1 (MQ=255)

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

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

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

                                                                                                                                      ggCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGc                <  1:5799428/150‑1 (MQ=255)

                                                                                                                                         cTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTAt             >  2:6190136/1‑150 (MQ=255)

                                                                                                                                           atggatggATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCt           >  1:2752882/1‑150 (MQ=255)

                                                                                                                                              gatggatgATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCa        >  2:4503142/1‑150 (MQ=255)

                                                                                                                                               atggatgATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAg       <  2:8028161/150‑1 (MQ=255)

                                                                                                                                                tggatgATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGc      <  2:3069486/150‑1 (MQ=255)

                                                                                                                                                  gatgatTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGa    >  2:6735077/1‑150 (MQ=255)

                                                                                                                                                  gatgatTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGa    <  2:8410497/150‑1 (MQ=255)

ATATCCACAAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACA  >  NC_002947/5030649‑5030946

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