BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,626,767	A→C	12.9%	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,626,767	0	A	C	12.9%	381.8 / 5.5	163	L229R (CTT→CGT)	PP_3195	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (78/59); new base C (2/23); total (80/83)
Fisher's exact test for biased strand distribution p-value = 4.34e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 169 total aligned reads displayed.

CCACTGCGCGCGGACAGTAGTAGCCGCCCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCGGGCCACTCACGCC  >  NC_002947/3626618‑3626915

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

ccACTGCGCGCGGACAGTAGTAGCCGCCCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGaa                                                                                                                                                      >  1:1887962/1‑150 (MQ=255)

ccACTGCGCGCGGACAGTAGTAGCCGCCCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGaa                                                                                                                                                      >  2:2792307/1‑150 (MQ=255)

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

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

  aCTGCGCGCGGACAGTAGTAGCCGCCCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAgc                                                                                                                                                    >  1:1049445/1‑150 (MQ=255)

  aCTGCGCGCGGACAGTAGTAGCCGCCCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAgc                                                                                                                                                    >  1:4758251/1‑150 (MQ=255)

   cTGCGCGCGGACAGTAGTAGCCGCCCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAgcg                                                                                                                                                   >  2:2276245/1‑150 (MQ=255)

   cTGCGCGCGGACAGTAGTAGCCGCCCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAgcg                                                                                                                                                   >  2:4846981/1‑150 (MQ=255)

     gcgcgcGGACAGTAGTAGCCGCCCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCt                                                                                                                                                 <  2:2419186/150‑1 (MQ=255)

            gACAGTAGTAGCCGCCCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGcctc                                                                                                                                          <  1:4468101/150‑3 (MQ=255)

                 tagtagCCGCCCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCg                                                                                                                                     >  2:2167789/1‑150 (MQ=255)

                 tagtagCCGCCCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCg                                                                                                                                     <  2:788272/150‑1 (MQ=255)

                      gccgccCTGGGAAACGCCTTCCAGGCGGATGCGCCTGGCTTCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCTTGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTg                                                                                                                                <  1:4161517/150‑1 (MQ=255)

                      gccgccCAGGGAAACGCCTTCCAGGCCGATGCGCGTGGCATCGACATCAGCACGGGTTTCCAGCCTGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGTTGGCGAAGCGCTTCGCCGGTGCCGGGCTg                                                                                                                                <  2:1299801/150‑1 (MQ=255)

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

                       ccgccCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTgg                                                                                                                               <  2:2221946/150‑1 (MQ=255)

                         gccCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAATCGCTTCGCCGGTGCCGGGCTGGtg                                                                                                                             >  1:221321/1‑149 (MQ=255)

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

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

                         gccCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCGGGTc                                                                                                                             >  1:4849064/1‑150 (MQ=255)

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

                                 aaaCGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCa                                                                                                                     <  2:2676296/150‑1 (MQ=255)

                                 aaaCGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCa                                                                                                                     <  2:2676619/150‑1 (MQ=255)

                                  aaCGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAg                                                                                                                    <  2:3830049/150‑1 (MQ=255)

                                   aCGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGc                                                                                                                   <  2:2013106/150‑1 (MQ=255)

                                     gCCTGCCCGGCCGATGCGCCTGGCATCGACAGCAGCCCGGGTTTCCTGCCAGTCCACTCGGCGGCGGGCCCTGTGTTCGCTGTCATGGCGTGCCGTCTTGCCGTGGAGGCGACGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGcga                                                                                                                 <  2:3842013/150‑1 (MQ=255)

                                      ccTTCCAGGCCGATGCGCCTGGCATCGACAACAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGcgac                                                                                                                >  2:3244295/1‑150 (MQ=255)

                                        ttCCAGGCCGATGCGCCTGGCAACGACATCAGCACGGGTTTCCAGCCAGTCCACTATGAGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGACGGGGCCGGGCTGGGCTACAAGAAGcgacga                                                                                                              >  1:4422180/1‑150 (MQ=255)

                                            aGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACAcg                                                                                                          >  2:1355550/1‑150 (MQ=255)

                                            aGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGGGCCGGGGTGGTCGACACTCAGCGACGACAcg                                                                                                          >  1:2582909/1‑150 (MQ=255)

                                            aGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGGGCCGGGCGGGTCGACAATCAGCGACGacacc                                                                                                          >  1:4883570/1‑149 (MQ=255)

                                                cGATGCGCGTGGCATCGACATCAGCTGGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCTGTGTTCGCTGTCTTGGGGTGCCGTCTTGCGGTGGTGGCGACGCGCTGCGCCGGCGCCGGGCTGGTCGACAATCAGCGACGACAcgcctc                                                                                                      <  2:3168848/150‑1 (MQ=255)

                                                     cgcCTGGCATCGACATCCGCCCGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGACGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTg                                                                                                 <  2:2992037/150‑1 (MQ=255)

                                                     cgcCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTg                                                                                                 <  2:523820/150‑1 (MQ=255)

                                                      gcCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTgg                                                                                                >  2:4535220/1‑150 (MQ=255)

                                                      gcCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTgg                                                                                                >  2:59322/1‑150 (MQ=255)

                                                         tGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAAACAGCGACGACACGCCTCGCCAGgcca                                                                                             >  1:4657735/1‑150 (MQ=255)

                                                             aTCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGCAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACCCCACGCCTGgccagcca                                                                                         >  1:3527968/1‑150 (MQ=255)

                                                                gACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCcaggc                                                                                      >  1:4179676/1‑150 (MQ=255)

                                                                 aCATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCcaggca                                                                                     <  2:1493291/150‑1 (MQ=255)

                                                                   aTCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCT‑GCCTGGCTGGCcaggcgagg                                                                                  <  2:2520128/150‑5 (MQ=255)

                                                                     cagccggggTTTCCAGCCAGTCCACTCGGCGGCTGGCCCTGTGTGCGCTGTCATGGCGTGCCGACTTGCCGTGGTGGCGACGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCcaggcaggca                                                                                 <  1:810286/144‑1 (MQ=255)

                                                                       gCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGACGCGCTTCGCCGGTGCCGGGCGGGTCGACCATCAGCGACGACACGCCGCGCCTGGCCCGCcaggcaggccgg                                                                               >  1:4204288/1‑150 (MQ=255)

                                                                       gCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCcaggcaggcagg                                                                               <  2:5012169/150‑1 (MQ=255)

                                                                        cACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCCGCCcggcaggcaggc                                                                              >  2:827721/1‑150 (MQ=255)

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

                                                                          cGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCg                                                                            >  2:1441480/1‑150 (MQ=255)

                                                                          cGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAAAAAGCCACGACACCCCGCGCCAGGCCCCCCAGGCAGgccggccg                                                                            >  1:431282/1‑150 (MQ=255)

                                                                           gggTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGa                                                                           >  1:1207643/1‑150 (MQ=255)

                                                                           gggTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGa                                                                           >  1:5004670/1‑150 (MQ=255)

                                                                           gggTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTAGCCGGTGCCGGGCGGGTCGACAAAAAGCGACGACACGCCACGCCAGGCCAGCCAGGCAGGCCGGCCGa                                                                           >  1:1747378/1‑150 (MQ=255)

                                                                           gggTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCGTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCACGCCTGGCCAGCCAGGCAGGCAGGCCGa                                                                           >  1:1747037/1‑150 (MQ=255)

                                                                            ggTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAAACAGCGACGACACCCCACGCCTGGCCACCCAGGCAGGCAGGCCGAc                                                                          >  1:4335293/1‑150 (MQ=255)

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

                                                                              tttCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCGGGTCGACAATCCGCGACGACACGCCACGCCTGGCCAGCCAGGCCGGCAGGCCGAcgc                                                                        >  1:332373/1‑150 (MQ=255)

                                                                                 cccgcccgTCCACTATGCGGGGGGCCCCGTGTTCGCTTTCATGGCGTGCCGTCAGGCCGTGGAGGCGACGCGCGTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGt                                                                     <  1:3291381/147‑1 (MQ=255)

                                                                                  cacccagCCCACTATGCGGCTCGCCCTGTGTTCGCTGTCGGGGCGGGCCGTCATGCCGTGGCGGCCAAGCGCTGCCCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTa                                                                    <  1:538523/147‑1 (MQ=255)

                                                                                   agccagTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCGTCGCCGGTGCCGGGCGGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTat                                                                   >  1:3138130/1‑150 (MQ=255)

                                                                                    gcctgTCCACTTGGGGGGGGGCCCCGTGTTCGCTGTGTTGGGTTGCCTTCATGCGGTGGCGGCGACGCGCTCCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTatt                                                                  <  1:1474941/150‑1 (MQ=255)

                                                                                     cccgcccACTATGCGGCTGCCCCTGTGTTCGCTGTCTTGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCCCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTattt                                                                 <  1:389976/145‑1 (MQ=255)

                                                                                     ccagTCCCCTGTGCGGCGGGCCCTGGTTGCGCTTTCTTGGCGTGCGTTCAGGCGGTGGCGGGGCGGGGCGTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTattt                                                                 <  1:3026564/150‑1 (MQ=255)

                                                                                           cACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGc                                                                                                                                  >  1:5112705/1‑79 (MQ=255)

                                                                                             ctcggcggcgGGCCCTGTGTTGGGTGTCTGGGCGGGCCGTCTGGCGGTGGCGCCGACGCGCTGCCCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCt                                                         <  1:2124290/146‑1 (MQ=255)

                                                                                              tATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCtt                                                        >  2:2760299/1‑150 (MQ=255)

                                                                                              tATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCGGGTCGACAATAAGCGACGACACGCCACGCCTGGCCAGCCCGGCAGGCAGGCCGACGCGGTATTTCAATTGCtt                                                        >  1:3910281/1‑150 (MQ=255)

                                                                                              tATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCGTAGCCGGGGCCGGGCGGGGCGACAATCAGCGACCACACGCCACGCCGGGCCACCCAGGCAGGCCGGCCCACGCGGTATTTAATTTCCtt                                                        >  1:454025/1‑150 (MQ=255)

                                                                                              naTGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCGTCGCCCGTGCCCGGCGGGGCGGCAATCAGCGACCACACGCCACGCCAGGCCAGCCCCGCAGGCAGGCCGCCGCCGTATGTAATTTccgg                                                        >  1:856153/2‑148 (MQ=255)

                                                                                               tggcggcgGGCCCAGTGTTCGGGTTACTGGCGTGCCGTCATGCGGTGGAGGCGACGCGCTCCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTg                                                       <  1:2940604/148‑1 (MQ=255)

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

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

                                                                                                ggcggcgGGCCCTGTGTTCGCTGTCATGGCGTGCCGTCAGGCGGTGGAGGCGAAGCGCTGCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTgg                                                      <  1:3758405/149‑1 (MQ=255)

                                                                                                 ggggCTGGCCCTGTGTTGTGTGTCCTGGCGTGCCGTCTTGCGGGGGCGGCTACGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGt                                                     <  1:972772/148‑1 (MQ=255)

                                                                                                  cggcggGCCCCGGGGGCGCTGTCCGGGCCTGCGGGCTGGCCGCGGCGGCGACGCGCGCCGCCGGTGCCGGGCTGGTCGCCACTCAGCGACGACCCGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCACTTCCCTGGtt                                                    <  1:2703126/150‑1 (MQ=255)

                                                                                                  cGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCCCTTCGCCGGTGCCGGGCTGGTCGACAAACAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGtt                                                    <  2:3346907/150‑1 (MQ=255)

                                                                                                  cGGCTGCCCCAGTGTTCGCTGTCTTGGCGTGCCGTCATGCCGTGGCGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGtt                                                    <  1:4629155/150‑1 (MQ=255)

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

                                                                                                        gCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATAAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTAATTTCCTTTGTTGAATcc                                              >  1:2460209/1‑150 (MQ=255)

                                                                                                           cAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCGTCGCCCGTGCCGGGCGGGTCGACAATCAGCGACGACACGCCTCGCCAGGCCAGCCACGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGc                                           >  1:2072168/1‑150 (MQ=255)

                                                                                                             gttttCGCTGTGTTGGCGTGCCGTCCTGCCGTGGCGGCGAAGGGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGGGGTATTTCATTTCCTTGGTTGAATCCAGccc                                         <  2:19357/147‑1 (MQ=255)

                                                                                                               gtTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGt                                       <  2:5118432/150‑1 (MQ=255)

                                                                                                                ttCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCggt                                      >  2:4402574/1‑148 (MQ=255)

                                                                                                                   gCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGCATCGCGTAGCCGGTGCCGGGCTGGTCGACAAAAAGCGACGACACCCCACCCCTGGCCACCCACGCAGGCAGGCCGACGCGGTATTTAATTTACTTGGTTGAATACAACCCCTTGAc                                   >  1:4373296/1‑150 (MQ=255)

                                                                                                                     tGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACtt                                 <  2:924650/150‑1 (MQ=255)

                                                                                                                     tGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACtt                                 >  2:1117811/1‑150 (MQ=255)

                                                                                                                     tGTCATGGCGTGCCGTCAGGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACtt                                 <  2:176044/150‑1 (MQ=255)

                                                                                                                       tgttGGCGTGCCGTATTGCGGTGGTGGCGACGCGCTTCGCCGGTGCCGGGCTGGTCGACACTCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGc                               <  1:3775242/147‑1 (MQ=255)

                                                                                                                       tCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCCGTGCCGGGCTGGTCGACAAACACCGACGACACGCCTCGCCAGGCCAGCCAGGCAGGCAGGCCGACGCCGTATTTCATTTACTTGGTTGAATCCAGCCCCTTGACTTGc                               >  1:3924585/1‑150 (MQ=255)

                                                                                                                           ggCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACca                           <  2:5004670/150‑1 (MQ=255)

                                                                                                                              gTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTTTTTCATTTACTTGGTTGTATCCAGCCCGTTGACTTGCACcagca                        >  1:3651416/1‑150 (MQ=255)

                                                                                                                              gTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACcagca                        >  2:4656587/1‑150 (MQ=255)

                                                                                                                                gCCGTCATGCCGTGGAGGCGACGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCagcagc                      <  2:617861/150‑1 (MQ=255)

                                                                                                                                       tGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCg               <  2:2582909/150‑1 (MQ=255)

                                                                                                                                        gCCGTGGAGGCGAAGCGCGTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACGTGCACCAGCAGCGGTGCCgg              >  1:4534342/1‑150 (MQ=255)

                                                                                                                                            tGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCGGGCCa          <  2:1302530/150‑1 (MQ=255)

                                                                                                                                                ggCGACGCGCTCCCCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCGGGCCACTCa      <  1:4783952/150‑1 (MQ=255)

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

                                                                                                                                                 gCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCGGGCCACTCAc     <  2:1911526/150‑1 (MQ=255)

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

                                                                                                                                                   gtaGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCGGGCCACTCACGc   <  2:2294069/148‑1 (MQ=255)

CCACTGCGCGCGGACAGTAGTAGCCGCCCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCGGGCCACTCACGCC  >  NC_002947/3626618‑3626915

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