BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,629,794	A→C	8.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	8.8%	881.5 / 9.3	338	L229R (CTT→CGT)	PP_3195	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (185/112); new base C (2/38); total (188/150)
Fisher's exact test for biased strand distribution p-value = 7.08e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 342 total aligned reads displayed.

CCACTGCGCGCGGACAGTAGTAGCCGCCCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCGGGCCACTCACGC  >  NC_002947/3629645‑3629941

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

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

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

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

  aCTGCGCGCGGACAGTAGTAGCCGCCCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAgc                                                                                                                                                   <  1:8549054/150‑1 (MQ=255)

    tGCGCGCGGACAGTAGTAGCCGCCCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAgcgc                                                                                                                                                 >  1:6897216/1‑150 (MQ=255)

       gcgcgGACAGTAGTAGCCGCCCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCAACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTc                                                                                                                                              <  2:3464324/150‑1 (MQ=255)

              cAGTAGTAGCCGCCCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGcctcca                                                                                                                                       <  2:5378863/150‑5 (MQ=255)

               agtagtAGCCGCCCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGGGAAGCGCTTCGCCGGTgg                                                                                                                                      >  1:7466636/1‑149 (MQ=255)

               agtagtAGCCGCCCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGTCGCGCTGCGCCGGTGc                                                                                                                                      >  1:5059450/1‑150 (MQ=255)

                 ttgtagCCGCCCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCg                                                                                                                                    <  2:2938031/148‑1 (MQ=255)

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

                      gccgccCAGGGAAACGCCTTCCAGGCCTTTGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCATTCCACTATGCGGCTGGCCCAGTTTTCGCTTTCATGGCTTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTg                                                                                                                               <  2:530396/150‑1 (MQ=255)

                      cccgccCAGGCCAACGCCTCCCCGGCCGAGGCGCCTGGCATCGACACCAGCACGGGTTTCCAGCCAGTCCCCTCTGCGGCGGGCCCAGTGTGCGCTGTCAGGGCGTGCCGTCTTGCCGTGGAGGCGACGCGCTTCGCCGGTGCCGGGCTg                                                                                                                               <  2:8257979/149‑1 (MQ=255)

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

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

                       ccgccCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCGTCGCCGGTGCcgggcggg                                                                                                                              >  1:4371779/1‑150 (MQ=255)

                        cgccCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCAGGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGt                                                                                                                             <  2:1097626/150‑1 (MQ=255)

                         gccCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTc                                                                                                                            >  2:8060313/1‑150 (MQ=255)

                                gAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCGTCGCCGGGGCCGGGCGGGTCGACAATc                                                                                                                     >  1:1610576/1‑150 (MQ=255)

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

                                   aCGCCTTCCAGGCCGATGCGCATGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGACGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCTACAATCAGc                                                                                                                  >  2:7426560/1‑150 (MQ=255)

                                     gCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAATCGCTTCGCCGGTGCCGGGCTGGTCGACAATCaccca                                                                                                                >  1:4196938/1‑146 (MQ=255)

                                     gCCGTCCCGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTAGGCGGCTGGCCCAGTGTTCGCTGTCATGGCGGGCCGTCATGCCGTGGCGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGcga                                                                                                                <  1:529937/150‑1 (MQ=255)

                                      ccTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGCAGCGCGTCGCCGGGGCCGGGCGGGGCGCCCAACAGcgcc                                                                                                               >  1:5326912/1‑148 (MQ=255)

                                      ccTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGcgac                                                                                                               >  2:1496608/1‑150 (MQ=255)

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

                                       atgccAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGcgacg                                                                                                              >  2:4071871/4‑150 (MQ=255)

                                          ccAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGaca                                                                                                           >  2:6377272/1‑150 (MQ=255)

                                                     cgcCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTTTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTg                                                                                                <  2:1673631/150‑1 (MQ=255)

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

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

                                                      gcCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCGGGTCGACAATCAGCGCCGACACGCCTCGCCTgg                                                                                               >  1:2744187/1‑150 (MQ=255)

                                                       ccTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGgc                                                                                              <  2:7847218/150‑1 (MQ=255)

                                                          ggCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGTCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGgccag                                                                                           >  2:1632501/1‑150 (MQ=255)

                                                           gCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGACATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGgccagc                                                                                          <  2:6128362/150‑1 (MQ=255)

                                                             aTCTACAACAGCAGGGGTTTCCAGCCCGCCAACTTGGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGCCATGCCGTGGCGGCGACGCGCTCCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGgccagcca                                                                                        <  1:3623398/150‑1 (MQ=255)

                                                               cGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCcagg                                                                                      <  2:2694110/150‑1 (MQ=255)

                                                                gACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCGTCGCCGGTGCCGGGCGGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCccggc                                                                                     >  1:2630333/1‑150 (MQ=255)

                                                                   aTCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGCAGCGCGTCGCCGGTGCCGGGCTGGTCGACAAACAGCGACCACACGCCTCCCCAGGCCAGCccggcagg                                                                                  >  1:3936147/1‑150 (MQ=255)

                                                                    tCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCAGGCCAGCcaggcaggc                                                                                 >  1:1286087/1‑150 (MQ=255)

                                                                     ccgCAGGGGTTTCCCGCCTGTCCACTTTGCGGCGGGCCCCGTGTCGGCTTTCTTGGCGGGCCGTCATGCCGTGGTGGCGACGCGCTCCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTTGCCTGGCTGGCcaggcgaggc                                                                                <  1:6388550/148‑6 (MQ=255)

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

                                                                       gCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCTGGCcaggcgaggcgt                                                                              <  2:2142250/150‑8 (MQ=255)

                                                                       gCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCcaggcaggcagg                                                                              >  2:1761298/1‑150 (MQ=255)

                                                                       gCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCCCCTGGCCAGCcaggcaggcagg                                                                              >  1:3341706/1‑150 (MQ=255)

                                                                         tgggggTTCCAGCCCGCCCCCTTGGCGGCTGGCCCCGTGTGCGCGGTCTGGGGGTGCCTTCCGGCGGTGGCGGCGCCGCGCTGCCCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGcc                                                                            <  1:171757/148‑1 (MQ=255)

                                                                         aCGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGcc                                                                            >  2:6207750/1‑150 (MQ=255)

                                                                         aCGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCCCTTCGCCGGTGCCGGGCTGGTCGACAAACAGCGACGACACGCCGCGCCAGGCCAGCCAGGCAGGCAGGcc                                                                            >  1:3023531/1‑150 (MQ=255)

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

                                                                                tCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCGTCGCCGGTGCCGGGCGGGTCGACCAACAGCGACGACACGCCTCGCCAGGCCCGCCAGGCAGGCAGGCCGACGCgg                                                                     >  1:604467/1‑150 (MQ=255)

                                                                                 ccagccagTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGt                                                                    <  2:8583694/150‑1 (MQ=255)

                                                                                 ccagccagTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCGGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGt                                                                    >  1:618798/1‑150 (MQ=255)

                                                                                    gccagTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTatt                                                                 >  2:4688680/1‑150 (MQ=255)

                                                                                    gccagTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTatt                                                                 >  2:5772991/1‑150 (MQ=255)

                                                                                      ccgTCCACTAGGCGGCTGGCCCCGTGTGCGCTGTCTTGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTatttc                                                               <  2:359430/148‑1 (MQ=255)

                                                                                      cagTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAAACACCGACGACACCCCTCGCCTGGCCAGCCACGCAGGCAGGCCGACCCGGTatttc                                                               >  1:8268848/1‑150 (MQ=255)

                                                                                       cgTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTatttca                                                              <  2:749184/149‑1 (MQ=255)

                                                                                          ccACTAGGCGGCTGGCCCAGTGTTCGCGGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTatttcattt                                                           <  2:5083351/150‑1 (MQ=255)

                                                                                            aCTAGGCGGCTGGCCCAGTGTTCGCTGTCAGGGCGTGCCGTCATGCCGGGGAGGCGACGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTcc                                                         <  2:4920151/150‑1 (MQ=255)

                                                                                             cTAGGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGACGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCt                                                        <  2:3818512/150‑1 (MQ=255)

                                                                                             cTAGGCGGCGGGCCCTGTGTTCGCTGTCTGGGCGTGCCGTCTTGCGGTGGCGGCGACGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCt                                                        <  1:4857219/150‑1 (MQ=255)

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

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

                                                                                               aTGCGGCTGGCCCCGTGTTCGATGTCTTGGCGTGCCGTCATGCCGTGGAGGCTAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCTCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTg                                                      <  1:2705378/150‑1 (MQ=255)

                                                                                               aTGCGGCGGGCCCTGTCTCGGCGTACATGCCGTGCCCACATGCCGTGGAGCCCACCCCCTCCCCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTg                                                      <  1:8792/150‑1 (MQ=255)

                                                                                                tggggCTGGCCCTGTGTTCGCTGTCATGGGGTGCCGTCTTGCGGTGGAGGCGCCGCGCTCCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTgg                                                     <  1:3418094/147‑1 (MQ=255)

                                                                                                ggcggggggCCCCGTGTGGGCTGTCTTGGCGTGCCGTCTGGCGGGGGCGGCGACGCGCGTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTgg                                                     <  1:4571037/143‑1 (MQ=255)

                                                                                                 gCGGCTGGCCCAGTTTTCGCTGTCTTGGGGTGCCGTCATGCGGTGGCGGCGACGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGt                                                    <  1:8028540/150‑1 (MQ=255)

                                                                                                  cGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGtt                                                   <  2:6090122/150‑1 (MQ=255)

                                                                                                  cGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGtt                                                   <  2:515292/150‑1 (MQ=255)

                                                                                                    gCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGa                                                 >  2:2124503/1‑150 (MQ=255)

                                                                                                    gCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGa                                                 >  2:5947248/1‑150 (MQ=255)

                                                                                                    gCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGa                                                 <  2:4285983/150‑1 (MQ=255)

                                                                                                    gCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTAGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGa                                                 <  2:3178979/150‑1 (MQ=255)

                                                                                                     cTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGaa                                                <  2:8669610/150‑1 (MQ=255)

                                                                                                      tGGCCTTGTGTTGGGTTTCTTGGGGTGCCGTCTTGCGGTGGCGGCGACGCGCTCCCCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAAt                                               <  1:215610/150‑1 (MQ=255)

                                                                                                        gCCCTGTGTTCGCTGCGAGGGCGTGCCGTCATGCGGTGGAGGCGCCGCGCTCCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATcc                                             <  1:3293456/150‑1 (MQ=255)

                                                                                                        gCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAAACAGCGACGACACGCCTCGCCAGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATcc                                             >  1:7060250/1‑150 (MQ=255)

                                                                                                          ccAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGGGCCCGGCGGGTCGACCATCAGCGACGACACGCCACGCCTGGCCAGCCCGGCAGGCAGGCCCACGCGGTATTTCATTTCCTTGGTTGAATcccc                                           >  1:4889674/1‑148 (MQ=255)

                                                                                                           cAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCGTGGGTGAATCCAGc                                          >  2:4646452/1‑150 (MQ=255)

                                                                                                                ttCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCGGGTCGACAATAAGCGACGACACCCCTCGCCTGGCCAGCCACGCAGGCAGGCCGACGCGGGATTTCATTTCCGTGGTTGAAAACAGCCCCtt                                     >  1:5132002/1‑150 (MQ=255)

                                                                                                                ttCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCGTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACCCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGtt                                     >  1:1962487/1‑150 (MQ=255)

                                                                                                                   ggtggcttGGCGGGCCTTCCTGGGGTGGTGGCGCCGCGCGCCCCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGAc                                  <  1:1018707/143‑1 (MQ=255)

                                                                                                                     tggcATGGCGTGCCGTCATGCCGGGGAGGCGAAGCGCATCGCCGGTGCCGCGCTGGTCGACAATCAGCCACGACACGCCTCCCCTGGCCAGCCAGGCACGCACGCCGACGCGGGAGTTCATATCCGTGGTTGAAAACAGCCCAGTGAcat                                >  2:499938/4‑148 (MQ=255)

                                                                                                                     tGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGGCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGTCGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCATTGTCtt                                >  2:390510/1‑150 (MQ=255)

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

                                                                                                                         aTGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCAc                            <  2:6551935/150‑1 (MQ=255)

                                                                                                                          tGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCACCGACGACACGCCACGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCGTGGTTGAAACCAGCCCCATGACATGCAcc                           >  1:8655279/1‑150 (MQ=255)

                                                                                                                           ggCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAAACAGCGACGACACGCCACGCCAGGCCAGCCAGGCAGGCCGGCCGACGCCGTATTTCATTTCCTGGGTTTAATCCAGCCCGTTGACTTGCACca                          >  1:7917329/1‑150 (MQ=255)

                                                                                                                                        gCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCCGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCgg             >  1:6557857/1‑150 (MQ=255)

                                                                                                                                         ccGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCCGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCCGCCCGTTGACGTGCACCCGCAGCGGTGCCggg            >  2:6841554/1‑150 (MQ=255)

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

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

                                                                                                                                              gAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCGGGCCACt       >  1:7313873/1‑150 (MQ=255)

                                                                                                                                              gAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCGGGCCACt       >  1:590227/1‑150 (MQ=255)

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

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

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

                                                                                                                                                   gAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCGGGCCACTCACGc  <  2:1953178/150‑1 (MQ=255)

CCACTGCGCGCGGACAGTAGTAGCCGCCCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCGGGCCACTCACGC  >  NC_002947/3629645‑3629941

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