BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,629,794	A→C	9.8%	L229R (CTT→CGT)	PP_3195 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,629,794	0	A	C	9.8%	841.3 / 12.0	332	L229R (CTT→CGT)	PP_3195	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (158/128); new base C (2/40); total (160/172)
Fisher's exact test for biased strand distribution p-value = 5.59e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 337 total aligned reads displayed.

CCACTGCGCGCGGACAGTAGTAGCCGCCCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCGGGCCACTCACGCC  >  NC_002947/3629645‑3629942

ccACTGCGCGCGGACAGTAGTAGCCGCCCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGaa                                                                                                                                                      <  2:2204147/150‑1 (MQ=255)

 cACTGCGCGCGGACAGTAGTAGCCGCCCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAg                                                                                                                                                     >  1:2353079/1‑150 (MQ=255)

       gcgcgGACAGTAGTAGCCGCCCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTc                                                                                                                                               <  1:5169661/150‑1 (MQ=255)

         gcgGACAGTAGTAGCCGCCCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGCAGCGCGTCGc                                                                                                                                             >  1:1647243/1‑150 (MQ=255)

         gcgGACAGTAGTAGCCGCCCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGc                                                                                                                                             <  2:6573785/150‑1 (MQ=255)

           ggACAGTAGTAGCCGCCCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCg                                                                                                                                           <  1:925925/150‑1 (MQ=255)

               agtagtAGCCGCCCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCTCGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCTTGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGcctccac                                                                                                                                       <  1:5968852/150‑6 (MQ=255)

               agtagtAGCCGCCCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGc                                                                                                                                       >  2:1865547/1‑150 (MQ=255)

               agtagtAGCCGCCCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGc                                                                                                                                       >  2:6754297/1‑150 (MQ=255)

                 tagtagCCGCCCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGCAGCGCGTCGCCGGTGCCg                                                                                                                                     >  1:7273930/1‑150 (MQ=255)

                      gccgccCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTg                                                                                                                                >  1:5406452/1‑150 (MQ=255)

                      gccgccCAGGGAAACGCATTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCAGTCATGCCGTGGAGGCGAAGAGCTTCTCCGGTGCCGGGCTg                                                                                                                                >  1:7239964/1‑150 (MQ=255)

                       ccgccCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTgg                                                                                                                               >  2:6289347/1‑150 (MQ=255)

                        cgccCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGt                                                                                                                              <  2:7091644/150‑1 (MQ=255)

                         gccCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTc                                                                                                                             <  2:5821951/150‑1 (MQ=255)

                          cccAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCg                                                                                                                            >  1:4282862/1‑150 (MQ=255)

                          cccAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCGTCGCCGGTGCCGGGCGGGTCg                                                                                                                            >  1:5055158/1‑150 (MQ=255)

                            cAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCgtc                                                                                                                          >  1:562497/1‑148 (MQ=255)

                             aGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACa                                                                                                                         >  2:1085935/1‑150 (MQ=255)

                               ggAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAAt                                                                                                                       >  2:2129382/1‑150 (MQ=255)

                               ggAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTGCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCGTCGCCGGTGCCGGGCTGGTCGACAAt                                                                                                                       >  1:2867201/1‑150 (MQ=255)

                                gAAACGCCTTCCAGGCCGTTGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTAGGCGGCTGGCCCTGTGTTCGCTGTCTTGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATc                                                                                                                      <  1:3680340/150‑1 (MQ=255)

                                 aaaCGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCa                                                                                                                     >  2:5248569/1‑150 (MQ=255)

                                   aCGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCAGGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGTAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGc                                                                                                                   <  2:3447468/150‑1 (MQ=255)

                                      ccTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGCAGCGCGTCGCCGGGGCCGGGCGGGGCGACCAGCAGcgac                                                                                                                >  1:6906157/1‑150 (MQ=255)

                                       cTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGCAGCGCATCGCCGGGGCCGGGCGGGTCGACAATCAGccacg                                                                                                               >  1:1238973/1‑150 (MQ=255)

                                       cTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGcgacg                                                                                                               >  2:1613171/1‑150 (MQ=255)

                                        ttCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCGCAGTGGTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGTCAATCAGcgacga                                                                                                              >  2:4426202/1‑150 (MQ=255)

                                            aGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTTGTCGACAATCAGCGACGACAcg                                                                                                          >  2:1318093/1‑150 (MQ=255)

                                             ggCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACAcgc                                                                                                         >  2:784868/1‑150 (MQ=255)

                                             ggCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACAcgc                                                                                                         <  2:2357674/150‑1 (MQ=255)

                                               cggATGCGCCTGGCATCTACATCCGCAGGGGTTCCCAGCCCGCCCACTTGGCGGCTGGCCGTGTGTGCGCTGTCTTGGCTTGCCTTCTTGCGGTGGCGGCGACGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACAcgcct                                                                                                       <  1:1973853/148‑1 (MQ=255)

                                                cgttGCGCCGGGCATCGACATCAGCACGGTTTTCCAGCCAGTCCACTTTGCGGCTGGCCCAGTGTTCGCTGTCTTGGGGTGCCGTCTTGCCGTGGTGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACAcgcctc                                                                                                      <  1:610106/147‑1 (MQ=255)

                                                 gATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGGGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACAcgcctcg                                                                                                     >  2:4652325/1‑150 (MQ=255)

                                                  aTGCGCCTGGAATCGACATCCGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTGGCTGTCATGGCGTGCCGTCATGCGGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACAcgcctcgc                                                                                                    <  2:39142/150‑1 (MQ=255)

                                                    gcgcCTGGCCTCGCCATCCGCCGGGGTTGCCAGCCCGCCCAGTATGGGGGTGGCCCTGTGTTGGCTGTCTTGGCGGGCCGTCATGCGGTGGAGGGGACGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACAcgcctcgcct                                                                                                  <  1:5307427/150‑1 (MQ=255)

                                                    gcgcCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCTGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACAcgcctcgcct                                                                                                  <  2:124270/150‑1 (MQ=255)

                                                     cgcCGGGCATCGACATCAGCTGGGTTTTCCAGCCAGACCACTTTGCGGCTGGCCCAGTGTTCGCTGTCTTGGGGTGCCGTCATGCGTTGGCGGCGAAGCGCTCCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTg                                                                                                 <  1:3017291/150‑1 (MQ=255)

                                                          ggCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGgccag                                                                                            <  2:546537/150‑1 (MQ=255)

                                                             aTCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGGGGGTCGACAATAAGCGACGACACGCCTCGCCTGgccagcca                                                                                         >  1:409099/1‑150 (MQ=255)

                                                              tCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCGGGCCCAGTGTTCGCTGTCAGGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGccagccag                                                                                        <  2:2926167/150‑1 (MQ=255)

                                                              ccGACATCCGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCATGGTTCGCTTTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGccagccag                                                                                        <  2:57922/149‑1 (MQ=255)

                                                                gACATCCGCACGGGTTACCAGCCAGTCCACTTTGCGGCTGGCCCAGTGTTCGCTGTCTTGGCGTGCCGTCATGCCGTGGCGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCcaggc                                                                                      <  1:6202236/150‑1 (MQ=255)

                                                                gACAACAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGCAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCcaggc                                                                                      <  2:5862896/150‑1 (MQ=255)

                                                                  cATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCcaggcag                                                                                    >  2:843228/1‑150 (MQ=255)

                                                                   aTCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCcaggcagg                                                                                   >  2:2303871/1‑150 (MQ=255)

                                                                     cAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTTGCCTGGCTGGCcaggcgaggc                                                                                 <  2:3745103/150‑6 (MQ=255)

                                                                     cAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCcaggcaggca                                                                                 <  2:1561879/150‑1 (MQ=255)

                                                                      aGCACGGTTTTCCACCCACTCCACTATGCGGCTGGCCAAGTTTTCGCTTTAATGGCGTGCCTTAATGCCGTGGAGGCTAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCcaggcaggcag                                                                                <  1:5550182/150‑1 (MQ=255)

                                                                      aGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGGGTTCGCTGTCATTGCGTGCCTTCATGCGGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCT‑GCCTGGCTGGCcaggcgaggcgt                                                                               <  2:7384966/150‑8 (MQ=255)

                                                                        ctcGGGTTACCAGCCAGTCCACTATGCGGGTGGCCCAGTGTTGGCTGTCTTGGCGTGCCGTCATGCCGGGGAGGCGCAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCaggcaggcaggc                                                                              <  1:2712956/148‑1 (MQ=255)

                                                                        cACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCaggcaggcaggc                                                                              <  2:1712887/150‑1 (MQ=255)

                                                                        cACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCGTAGCCCGGGCCGGGCGGGGCGACCAACAGCGACCACACGCCGCCCCAGGCCACCCacgcacgcaggc                                                                              >  1:2047122/1‑150 (MQ=255)

                                                                          cGGGTTTCCAGCCAGTCCCCTAGGCGGCTGGCCCAGGGTTCGCTGTCAGGGCGGGCCGTCATGCCGGGGCGGCGCCGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCg                                                                            <  2:5090105/150‑1 (MQ=255)

                                                                           gggTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCACGCCTGGCCCGCCACGCAGGCAGGCCGa                                                                           >  1:365203/1‑150 (MQ=255)

                                                                             gTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGAcg                                                                         >  2:7236170/1‑150 (MQ=255)

                                                                             gTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGAcg                                                                         >  2:4151904/1‑150 (MQ=255)

                                                                               ttCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTACGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGAcgcg                                                                       <  1:1021859/150‑1 (MQ=255)

                                                                                   agcccgCCCACTCGGGGGGTGGCCCCGTGTTGGGTGTCTTGGGGTGCCGTCAGGCGGTGGCGGCGACGCGCTGCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTat                                                                   <  1:3897444/150‑1 (MQ=255)

                                                                                   agccagTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCGGGTCGACAAACACCGACGACACCCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGGat                                                                   >  1:4692468/1‑150 (MQ=255)

                                                                                         tCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCCACAATCAGCGACGACACGCCTCGCCAGGCCAGCCAGGCAGGCAGGCCGACGCCGTatttcatt                                                             >  1:1758486/1‑150 (MQ=255)

                                                                                         tCCACTATGCGGCTGGCCCAGTGGTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTatttcatt                                                             <  2:5999829/150‑1 (MQ=255)

                                                                                              tATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCtt                                                        <  2:2954731/150‑1 (MQ=255)

                                                                                               tggcggcgGGCCGTGTGTTCGCTGTCTGGGGGTGCCGTCTGGCGGGGGCGGCGATGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTg                                                       <  1:5925576/148‑1 (MQ=255)

                                                                                               aTGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTg                                                       <  2:3972318/150‑1 (MQ=255)

                                                                                                tGCGGCTGGCCCCGTGGTCGCGGTCAGGGCGGGCCGTCATGCCGTGGAGGCGACGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTgg                                                      <  2:3163785/150‑1 (MQ=255)

                                                                                                tGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTgg                                                      >  2:2181718/1‑150 (MQ=255)

                                                                                                 ggggCGGGCCCTGTGCTGGGTTTCTTGGCGCGCCGTCATGCGGGGGCGGCCCCGCGCCCCCCCGGCGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGt                                                     <  1:145365/148‑1 (MQ=255)

                                                                                                        gCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATcc                                              <  2:4241734/150‑1 (MQ=255)

                                                                                                         cccAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGTCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCa                                             >  2:6886666/1‑150 (MQ=255)

                                                                                                          ccAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCAGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAg                                            >  1:2255056/1‑150 (MQ=255)

                                                                                                          ccAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCGTAGCCGGTGCCGGGCGGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAg                                            >  1:790618/1‑150 (MQ=255)

                                                                                                           ctgtgtTCGTTGTCTTGGGGTGCGGTCAGGCGGTGGCGGCGACGCGCTCCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGc                                           <  1:5797142/148‑1 (MQ=255)

                                                                                                           cAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGc                                           >  2:5520981/1‑150 (MQ=255)

                                                                                                               gtgcGCTGTCTTGGGGTGCCGTCATGCCGTGGCGGCTACGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGt                                       <  1:6396797/147‑1 (MQ=255)

                                                                                                                ttCGCTGTCATGGCGTGCCGTCCTGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCggt                                      >  2:5174825/1‑148 (MQ=255)

                                                                                                                   gCTGTCAGGGCGTGCCGTAATGCCGTGGCGGCCAAGCCCTCCCCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGAc                                   <  1:6239249/150‑1 (MQ=255)

                                                                                                                    cTGTCATGGCGTGCCGTCAGGCCGTGGAGGCGAAGCGCTGCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCTGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACt                                  <  2:1678591/150‑1 (MQ=255)

                                                                                                                     tGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCCGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCAGGCCAGCCAGGCAGGCAGGCCGACGCGGGATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACtt                                 >  1:333143/1‑150 (MQ=255)

                                                                                                                       tCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAAACAGCGACGACACGCCTCGCCAGGCCAGCCAGGCAGGCAGGCCGACGCGGTCTTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGc                               >  1:484401/1‑150 (MQ=255)

                                                                                                                         aTGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTAGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACCCCTCGCCTGGCCAGCCAGGCCGGCAGGCCGACGCGGTATTTCAATTACGTGGTTGAATCCAGCCCGTTTACTTGCAc                             >  1:1079787/1‑150 (MQ=255)

                                                                                                                         aTGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTACGCCGGTGCCCGGCGGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCAc                             >  1:603808/1‑150 (MQ=255)

                                                                                                                          tGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAAAAAGCGACGACACGCCTCGCCAGGCCAGCCAGGCAGGCAGGCCGACGCGGGATTTCATTTACATGGTTGAATCCAGCCCCTTGACTTTCAcc                            >  1:4637801/1‑150 (MQ=255)

                                                                                                                           ggCGTGCCGTCATGCCGTGGCGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACca                           <  1:4196960/150‑1 (MQ=255)

                                                                                                                             cGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACcagc                         <  2:2502797/150‑1 (MQ=255)

                                                                                                                                 ccGTCAGGCCGTGGAGGCGACGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCg                     <  2:353429/150‑1 (MQ=255)

                                                                                                                                   gTCCTGCGGGGGCGGCGACGCGCTTCCCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGt                   <  1:6999589/150‑1 (MQ=255)

                                                                                                                                     cATGCGGTGGAGGCGAAGCGCTTCCCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGc                 <  1:2600257/150‑1 (MQ=255)

                                                                                                                                      aTGCCGTGGAGGCGAAGCGCATCGCCGGTGCCGGGCAGGTCGACAATAAGCGACGACACGCCACGCCAGGCCATCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCATTGTTTAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGcc                >  1:5239759/1‑150 (MQ=255)

                                                                                                                                          cGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCGGGc            <  2:5353338/150‑1 (MQ=255)

                                                                                                                                          cGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCGGGc            >  2:4803931/1‑150 (MQ=255)

                                                                                                                                           gTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCGGGcc           >  2:3404346/1‑150 (MQ=255)

                                                                                                                                           gTGGAGGCGAAGCGCTACGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCGGGcc           >  1:3259607/1‑150 (MQ=255)

                                                                                                                                            tGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCGGGCCa          >  1:3144650/1‑150 (MQ=255)

                                                                                                                                             ggAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCGGGCCAc         <  1:4875731/150‑1 (MQ=255)

                                                                                                                                               cggCGACGCGCTCCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCGGGCCACTc       <  1:2125831/149‑1 (MQ=255)

                                                                                                                                                ggCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGAAAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCGGGCCACTCa      >  1:6698735/1‑150 (MQ=255)

                                                                                                                                                 gCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCGGGCCACTCAc     >  1:3453118/1‑150 (MQ=255)

                                                                                                                                                  cGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCGGGCCACTCACg    >  2:4937533/1‑150 (MQ=255)

                                                                                                                                                   gAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCGGGCCACTCACGc   >  1:7438172/1‑150 (MQ=255)

CCACTGCGCGCGGACAGTAGTAGCCGCCCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCGGGCCACTCACGCC  >  NC_002947/3629645‑3629942

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