BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,754,931	C→T	8.6%	A307A (GCC→GCT)	glnA →	glutamine synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,754,931	0	C	T	8.6%	439.0 / 22.9	151	A307A (GCC→GCT)	glnA	glutamine synthetase
Reads supporting (aligned to +/- strand): ref base C (72/66); new base T (11/2); total (83/68)
Fisher's exact test for biased strand distribution p-value = 3.82e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 152 total aligned reads displayed.

GCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGA  >  NC_002947/5754784‑5755080

gCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCt                                                                                                                                                     >  2:3944354/1‑150 (MQ=255)

   aCTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGaa                                                                                                                                                  >  1:4203126/1‑150 (MQ=255)

    cTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAAc                                                                                                                                                 <  2:2667135/150‑1 (MQ=255)

      gTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACgg                                                                                                                                               <  1:1376877/150‑1 (MQ=255)

       tACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGc                                                                                                                                              >  2:1806377/1‑150 (MQ=255)

         cGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCt                                                                                                                                                     <  2:1752387/141‑1 (MQ=255)

          ggCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTc                                                                                                                                           <  1:3371683/150‑1 (MQ=255)

            cGACAACGGCTCGGGCATGCATGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCAc                                                                                                                                         <  1:2229803/150‑1 (MQ=255)

            cGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCAc                                                                                                                                         <  2:4659799/150‑1 (MQ=255)

            cGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCAc                                                                                                                                         <  2:3254073/150‑1 (MQ=255)

                  cGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAAccc                                                                                                                                   >  1:3586627/1‑150 (MQ=255)

                   ggCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCg                                                                                                                                  <  1:3039082/150‑1 (MQ=255)

                   ggCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCg                                                                                                                                  >  1:197947/1‑150 (MQ=255)

                   ggCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCg                                                                                                                                  >  1:1032759/1‑150 (MQ=255)

                        gggCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCtt                                                                                                                                            <  2:2739447/135‑1 (MQ=255)

                        gggCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCtt                                                                                                                                            >  1:2739447/1‑135 (MQ=255)

                              gCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTc                                                                                                                       <  1:94926/150‑1 (MQ=255)

                                aCGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGt                                                                                                                     >  2:1858787/1‑150 (MQ=255)

                                    gCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACaa                                                                                                                 <  2:1943905/150‑1 (MQ=255)

                                    gCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACaa                                                                                                                 <  2:3551557/150‑1 (MQ=255)

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

                                        ttctCGATCTGGAAAGAAGGCAAGCCAACGTTAGCTGGTAACGAATATGCCGGCCTGTCGGAAAACTGCCTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCAATAAGCGt                                                                                                             >  2:4021242/4‑150 (MQ=255)

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

                                        aTGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGt                                                                                                             <  1:4732319/150‑1 (MQ=255)

                                            cGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTgg                                                                                                         <  1:1849958/150‑1 (MQ=255)

                                            cGATCTGGAAAGAAGGCAAGCCAACGTTTGCTGGTAACGAATATGCCGGCCTGTCGGAAAACTGCCTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATAAGCGTCTgg                                                                                                         <  1:1113884/150‑1 (MQ=255)

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

                                             gATCTGGAAAGAAGGCAAGCCAACGTTTGCTGGTAACGAATATGCCGGCCTGTCGGAAAACTGCCTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATAAGCGTCTGGt                                                                                                        >  2:2501139/1‑150 (MQ=255)

                                                cTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTccc                                                                                                     <  1:2630438/150‑1 (MQ=255)

                                                 tGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCa                                                                                                    <  1:5240303/150‑1 (MQ=255)

                                                 tGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCa                                                                                                    <  1:5240338/150‑1 (MQ=255)

                                                 tGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCa                                                                                                    >  2:230523/1‑150 (MQ=255)

                                                 tGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCa                                                                                                    >  2:230665/1‑150 (MQ=255)

                                                 tGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCa                                                                                                    >  1:3932204/1‑150 (MQ=255)

                                                  ggAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGc                                                                                                               <  2:4350140/138‑1 (MQ=255)

                                                  ggAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAg                                                                                                   >  2:2011858/1‑150 (MQ=255)

                                                   gAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAgg                                                                                                  <  2:3341850/150‑1 (MQ=255)

                                                     aaGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCt                                                                                                >  2:4402968/1‑150 (MQ=255)

                                                          ggCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCt                                                                                                                                                     >  2:4600828/1‑92 (MQ=255)

                                                          ggCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCt                                                                                                                                                     <  1:4600828/92‑1 (MQ=255)

                                                          ggCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGATCGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGaa                                                                                           <  2:3429364/150‑1 (MQ=255)

                                                          ggCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGaa                                                                                           >  1:445212/1‑150 (MQ=255)

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

                                                          ggCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGaa                                                                                           >  1:1548501/1‑150 (MQ=255)

                                                           gCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAg                                                                                          <  1:4322215/150‑1 (MQ=255)

                                                           gCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAg                                                                                          <  2:3766908/150‑1 (MQ=255)

                                                               gAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGcccc                                                                                      <  1:3644694/150‑1 (MQ=255)

                                                                aaCACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCg                                                                                     <  1:1858787/150‑1 (MQ=255)

                                                                     cTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAAt                                                                                >  1:3946779/1‑150 (MQ=255)

                                                                         tcGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTc                                                                                                                       >  1:3356326/1‑107 (MQ=255)

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

                                                                         tcGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTg                                                                            >  1:2584311/1‑150 (MQ=255)

                                                                                aaGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACt                                                                     <  1:230523/150‑1 (MQ=255)

                                                                                aaGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACt                                                                     >  1:4250248/1‑150 (MQ=255)

                                                                                  ggCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCg                                                                   >  1:630356/1‑150 (MQ=255)

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

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

                                                                                     tATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGAcc                                                                                                                            <  2:5334272/90‑1 (MQ=255)

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

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

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

                                                                                       tGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCg                                                              >  1:1375584/1‑150 (MQ=255)

                                                                                        gCCGGCCTGTCGGAAAACTGCCTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATAAGCGTCTGGTCCCAGGTTATGAAGCTCCGGTTCTGCTGGCCTATTCGGCTCgt                                                             >  1:3088796/1‑149 (MQ=255)

                                                                                        gCCGGCCTGTCGGAAAACTGCCTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATAAGCGTCTGGTCCCAGGTTATGAAGCTCCGGTTCTGCTGGCCTATTCGGCTCgt                                                             >  1:3088834/1‑149 (MQ=255)

                                                                                        gCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGc                                                             >  1:1798740/1‑150 (MQ=255)

                                                                                         ccGGCCTGTCGGAAAACTGCCTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATAAGCGTCTGGTCCCAGGTTATGAAGCTCCGGTTCTGCTGGCCTATTCGGCTCgta                                                            >  1:890907/1‑148 (MQ=255)

                                                                                          cGGCCTGTCGGAAAACTGCCTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATAAGCGTCTGGTCCCAGGTTATGAAGCTCCGGTTCTGCTGGCCTATTCGGCTCGTaa                                                           >  1:586/1‑150 (MQ=255)

                                                                                           ggCCTGTCGGAAAACTGCCTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATAAGCGTCTGGTCCCAGGTTATGAAGCTCCGGTTCTGCTGGCCTATTCGGCTCGTAAc                                                          >  1:1861110/1‑150 (MQ=255)

                                                                                            gCCTGTCGGAAAACTGCCTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATAAGCGTCTGGTCCCAGGTTATGAAGCTCCGGTTCTGCTGGCCTATTCGGCTCGTAAcc                                                         >  2:2154271/1‑150 (MQ=255)

                                                                                              cTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGt                                                       >  1:1600652/1‑150 (MQ=255)

                                                                                               tGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGtt                                                      >  1:4913681/1‑150 (MQ=255)

                                                                                                 tCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGc                                                                          <  2:2150189/128‑1 (MQ=255)

                                                                                                 tCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGc                                                                          >  1:2150189/1‑128 (MQ=255)

                                                                                                        ccGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGa                                             >  1:3814812/1‑150 (MQ=255)

                                                                                                          tgccTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATAAGCGTCTGGTCCCAGGTTATGAAGCTCCGGTTCTGCTGGCCTATTCGGCTCGTAACCGTTCCGCTTcctgc                                           >  2:4940629/3‑146 (MQ=255)

                                                                                                            ccTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATAAGCGTCTGGTCCCAGGTTATGAAGCTCCGGTTCTGCTGGCCTATTCGGCTCGTAACCGTTCCGCTTcctgccg                                         >  1:1920139/1‑144 (MQ=255)

                                                                                                             cTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGt                                        >  1:4033511/1‑150 (MQ=255)

                                                                                                              tGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTa                                       <  2:523877/150‑1 (MQ=255)

                                                                                                                  cTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTcc                                   >  2:4751188/1‑150 (MQ=255)

                                                                                                                    tCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCtt                                 <  2:1574506/150‑1 (MQ=255)

                                                                                                                     cATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTc                                                                                                                                <  2:1554871/54‑1 (MQ=255)

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

                                                                                                                            ggTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGc                         <  2:3039362/150‑1 (MQ=255)

                                                                                                                             gTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCTACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCa                        <  1:727067/150‑1 (MQ=255)

                                                                                                                             gTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCa                        <  2:2844833/150‑1 (MQ=255)

                                                                                                                                tcatcaAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGcc                     <  2:3689691/150‑1 (MQ=255)

                                                                                                                                 catcaAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGccc                    <  2:3932204/150‑1 (MQ=255)

                                                                                                                                  atcaAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTg                                                                            <  1:3722822/93‑1 (MQ=255)

                                                                                                                                  atcaAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTcc                                                    >  2:4771212/1‑117 (MQ=255)

                                                                                                                                  atcaAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTcc                                                    <  1:4771212/117‑1 (MQ=255)

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

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

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

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

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

                                                                                                                                     aaGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGaaa                >  1:2693364/1‑150 (MQ=255)

                                                                                                                                      aGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGaaa                >  2:4168868/1‑149 (MQ=255)

                                                                                                                                      aGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGaaa                <  1:4168868/149‑1 (MQ=255)

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

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

GCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGACGTATCGA  >  NC_002947/5754784‑5755080

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