BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,754,967	C→T	14.9%	Y319Y (TAC→TAT)	glnA →	glutamine synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,754,967	0	C	T	14.9%	363.8 / 47.2	147	Y319Y (TAC→TAT)	glnA	glutamine synthetase
Reads supporting (aligned to +/- strand): ref base C (74/51); new base T (11/11); total (85/62)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 149 total aligned reads displayed.

GCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCATCGGCCAACCCGTACCT  >  NC_002947/5754820‑5755116

tcacttctCGATCTGGAAAGAAGGCAAGCCAACGTTTGCTGGTAACGAATATGCCGGCCTGTCGGAAAACTGCCTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATaa                                                                                                                                                     <  2:177099/143‑1 (MQ=255)

tcacttctCGATCTGGAAAGAAGGCAAGCCAACGTTTGCTGGTAACGAATATGCCGGCCTGTCGGAAAACTGCCTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATaa                                                                                                                                                     <  2:1281151/143‑1 (MQ=255)

  acaTGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACaatc                                                                                                                                                   >  2:416701/1‑148 (MQ=255)

  acaTGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGc                                                                                                                                                   >  2:696580/1‑150 (MQ=255)

  acaTGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGc                                                                                                                                                   <  1:137359/150‑1 (MQ=255)

  acaTGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGc                                                                                                                                                   >  2:448607/1‑150 (MQ=255)

  acaTGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGc                                                                                                                                                   >  2:3241057/1‑150 (MQ=255)

    ttctCGATCTGGAAAGAAGGCAAGCCAACGTTTGCTGGTAACGAATATGCCGGCCTGTCGGAAAACTGCCTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATAAGCGt                                                                                                                                                 >  1:748965/4‑150 (MQ=255)

    aTGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGt                                                                                                                                                 >  1:2055028/1‑150 (MQ=255)

     tGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTc                                                                                                                                                <  1:2094032/150‑1 (MQ=255)

     tGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTc                                                                                                                                                >  1:1887770/1‑150 (MQ=255)

     tGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGATTCACCAACCCGTCGACCAACTCGTACAAGCGTc                                                                                                                                                <  1:1250506/150‑1 (MQ=255)

      gTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCt                                                                                                                                               >  2:1140819/1‑150 (MQ=255)

      ctCGATCTGGAAAGAAGGCAAGCCAACGTTTGCTGGTAACGAATATGCCGGCCTGTCGGAAAACTGCCTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATAAGCGTCt                                                                                                                                               >  2:4582686/2‑150 (MQ=255)

         gATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGt                                                                                                                                            >  1:4843381/1‑150 (MQ=255)

         gATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGt                                                                                                                                            >  1:482366/1‑150 (MQ=255)

          aTCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTc                                                                                                                                           >  2:4402304/1‑150 (MQ=255)

          aTCTGGAAAGAAGGCAAGCCAACGTTTGCTGGTAACGAATATGCCGGCCTGTCGGAAAACTGCCTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATAAGCGTCTGGTc                                                                                                                                           >  2:1049903/1‑150 (MQ=255)

            cTGGAAAGAAGGCAAGCCAACGTTTGCTGGTAACGAATATGCCGGCCTGTCGGAAAACTGCCTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATAAGCGTCTGGTccc                                                                                                                                         <  1:4160340/150‑1 (MQ=255)

                  aGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCtt                                                                                                                                   >  1:4139359/1‑150 (MQ=255)

                      ggCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGTGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGaa                                                                                                                               <  2:2306260/150‑1 (MQ=255)

                      ggCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGaa                                                                                                                               <  1:1633489/150‑1 (MQ=255)

                      ggCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGaa                                                                                                                               <  1:902172/150‑1 (MQ=255)

                       gCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAg                                                                                                                              >  2:2472470/1‑150 (MQ=255)

                         aaGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGcc                                                                                                                            >  1:1999939/1‑150 (MQ=255)

                           gAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGcccc                                                                                                                          >  2:4290315/1‑150 (MQ=255)

                               acCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTa                                                                                                                      >  1:4928409/1‑150 (MQ=255)

                                   tctcGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGc                                                                                                                  >  2:2964333/1‑150 (MQ=255)

                                   tctcGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGc                                                                                                                  >  2:3810017/1‑150 (MQ=255)

                                     tcGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTg                                                                                                                <  1:2845382/150‑1 (MQ=255)

                                     tcGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTg                                                                                                                <  1:3692842/150‑1 (MQ=255)

                                          taacgaatATGCCGGCCTGTCGGAAAACTGCCTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATAAGCGTCTGGTCCCAGGTTATGAAGCTCCGGTTCTGCTGGCCTa                                                                                                           >  1:3090569/8‑150 (MQ=255)

                                             aGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTc                                                                                                        >  2:3151830/1‑150 (MQ=255)

                                              gaatATGCCGGCCTGTCGGAAAACTGCCTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATAAGCGTCTGGTCCCAGGTTATGAAGCTCCGGTTCTGCTGGCCTATTCg                                                                                                       >  1:2706268/4‑150 (MQ=255)

                                               gCTATGCCGTCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCgg                                                                                                      >  1:4682483/1‑150 (MQ=255)

                                               gCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCgg                                                                                                      >  2:3142784/1‑150 (MQ=255)

                                               gCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCgg                                                                                                      >  2:586122/1‑150 (MQ=255)

                                                 tATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCt                                                                                                    >  1:4118091/1‑150 (MQ=255)

                                                  aTGCCGGCCTGTCGGAAAACTGCCTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATAAGCGTCTGGTCCCAGGTTATGAAGCTCCGGTTCTGCTGGCCTATTCGGCTc                                                                                                   <  2:4368945/150‑1 (MQ=255)

                                                  aTGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTc                                                                                                   <  1:5096372/150‑1 (MQ=255)

                                                  aTGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTc                                                                                                   <  2:3212508/150‑1 (MQ=255)

                                                  aTGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTc                                                                                                   <  1:111699/150‑1 (MQ=255)

                                                  aTGCCGCCCTGTCCGACACCGGCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTc                                                                                                   <  1:223332/150‑1 (MQ=255)

                                                      cGGCCTGTCGGAAAACTGCCTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATAAGCGTCTGGTCCCAGGTTATGAAGCTCCGGTTCTGCTGGCCTATTCGGCTCGTaa                                                                                               <  2:2719284/150‑1 (MQ=255)

                                                      cGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCaa                                                                                               >  1:4893180/1‑150 (MQ=255)

                                                            gTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTc                                                                                         <  1:3086335/150‑1 (MQ=255)

                                                            gTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGATGGCCTACTCGGCTCGCAACCGTTc                                                                                         <  1:2328743/150‑1 (MQ=255)

                                                             tCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTcc                                                                                        >  2:419579/1‑150 (MQ=255)

                                                             tCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTcc                                                                                        <  1:1459705/150‑1 (MQ=255)

                                                               cGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGc                                                                                      <  1:1615779/150‑1 (MQ=255)

                                                                gACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGcc                                                                                     >  1:1939491/1‑150 (MQ=255)

                                                                     cGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGAt                                                                                >  1:2303692/1‑150 (MQ=255)

                                                                        ccTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATAAGCGTCTGGTCCCAGGTTATGAAGCTCCGGTTCTGCTGGCCTATTCGGCTCGTAACCGTTCCGCTTcctgccg                                                                             <  1:3597516/150‑7 (MQ=255)

                                                                        ccTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATAAGCGTCTGGTCCCAGGTTATGAAGCTCCGGTTCTGCTGGCCTATTCGGCTCGTAACCGTTCCGCTTcctgccg                                                                             <  1:874061/150‑7 (MQ=255)

                                                                         cTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATAAGCGTCTGGTCCCAGGTTATGAAGCTCCGGTTCTGCTGGCCTATTCGGCTCGTAACCGTTCCGCTTcctgccgt                                                                            >  2:1799262/1‑143 (MQ=255)

                                                                           gTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTAt                                                                          <  2:2922998/150‑1 (MQ=255)

                                                                             tcTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATAAGCGTCTGGTCCCAGGTTATGAAGCTCCGGTTCTGCTGGCCTATTCGGCTCGTAACCGTTCCGCTTCCTGCCGTATTc                                                                        <  1:1479215/149‑1 (MQ=255)

                                                                              cTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTcc                                                                       <  1:627922/150‑1 (MQ=255)

                                                                               ttCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCt                                                                      >  1:2382176/1‑150 (MQ=255)

                                                                                   tCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACg                                                                  >  2:1246615/1‑150 (MQ=255)

                                                                                   tCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACg                                                                  >  1:3728937/1‑150 (MQ=255)

                                                                                     ggcggTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGtt                                                                <  2:2610940/150‑1 (MQ=255)

                                                                                     ggcggTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGtt                                                                <  2:2610682/150‑1 (MQ=255)

                                                                                               tcaAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGa                                                      >  1:2964169/1‑150 (MQ=255)

                                                                                               tcaAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGa                                                      >  1:2584009/1‑150 (MQ=255)

                                                                                                 aaGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGaaa                                                    <  1:4331578/150‑1 (MQ=255)

                                                                                                 aaGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGaaa                                                    <  2:140013/150‑1 (MQ=255)

                                                                                                 aaGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGGTCCGCCGCGATCCGGATTCCTTACGTTGGCAGCCCGaaa                                                    <  2:1098345/150‑1 (MQ=255)

                                                                                                   gCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTACTCCTTACGTTGGCAGCCCGAAAGc                                                  >  1:3360346/1‑150 (MQ=255)

                                                                                                     aCGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGccc                                                <  1:4958914/150‑1 (MQ=255)

                                                                                                       ggTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGa                                              <  2:1220317/150‑1 (MQ=255)

                                                                                                       ggTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGa                                              >  2:33574/1‑150 (MQ=255)

                                                                                                        gTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGAc                                             <  2:2870125/150‑1 (MQ=255)

                                                                                                         tAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACg                                            >  1:1838797/1‑150 (MQ=255)

                                                                                                         tAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACg                                            >  1:1838456/1‑150 (MQ=255)

                                                                                                            ggCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTAt                                         <  1:2004332/150‑1 (MQ=255)

                                                                                                                cTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGaa                                     <  1:696580/150‑1 (MQ=255)

                                                                                                                cTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGaa                                     <  1:1738925/150‑1 (MQ=255)

                                                                                                                 tGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAg                                    <  2:1331829/150‑1 (MQ=255)

                                                                                                                    aCGGCTTCACCAACCCGTCGACCAACTCGGACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCAc                                 <  2:1876748/150‑1 (MQ=255)

                                                                                                                       gCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCt                              >  1:4405296/1‑150 (MQ=255)

                                                                                                                        cTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCtt                             >  1:3260861/1‑150 (MQ=255)

                                                                                                                         ttCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTc                            >  1:4563676/1‑150 (MQ=255)

                                                                                                                         ttCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCATTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTc                            >  1:194592/1‑150 (MQ=255)

                                                                                                                           cACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTccc                          <  2:2570695/150‑1 (MQ=255)

                                                                                                                            aCCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCg                         >  2:123767/1‑150 (MQ=255)

                                                                                                                            aCCAACCCGTCGACCAACGCGTACCAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCg                         >  2:2891922/1‑150 (MQ=255)

                                                                                                                             ccAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCgg                        >  2:4905795/1‑150 (MQ=255)

                                                                                                                                aCCCGTCGACCCACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGAcc                     >  1:3746977/1‑150 (MQ=255)

                                                                                                                                  ccGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGAAGCCCGAAAGCCCGAAGTATCGAAGCACGCTTCCCGGACCCa                   >  1:209191/1‑150 (MQ=255)

                                                                                                                                          cAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGtt                                                                <  2:1838712/97‑1 (MQ=255)

                                                                                                                                          cAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCATCGGCCaa           <  2:4139359/150‑1 (MQ=255)

                                                                                                                                          cAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCATCGGCCaa           >  2:3372754/1‑150 (MQ=255)

                                                                                                                                          cAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCATCGGCCaa           >  2:1645727/1‑150 (MQ=255)

                                                                                                                                           aaCTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCATCGGCCAAc          <  1:4402304/150‑1 (MQ=255)

                                                                                                                                           aaCTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCATCGGCCAAc          >  2:2649486/1‑150 (MQ=255)

                                                                                                                                           aaCTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCATCGGCCAAc          <  2:4928409/150‑1 (MQ=255)

                                                                                                                                                 tACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCATCGGCCAACCCGTAc    >  1:2723721/1‑150 (MQ=255)

                                                                                                                                                  aCAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCATCGGCCAACCCGTAcc   >  2:2210810/1‑150 (MQ=255)

                                                                                                                                                  aCAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCATCGGCCAACCCGTAcc   >  2:3638152/1‑150 (MQ=255)

GCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGAAGCACGCTTCCCGGACCCATCGGCCAACCCGTACCT  >  NC_002947/5754820‑5755116

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