BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,754,924	G→C	8.6%	G305A (GGT→GCT) ‡	glnA →	glutamine synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,754,924	0	G	C	8.6%	486.1 / 28.4	151	G305A (GGT→GCT) ‡	glnA	glutamine synthetase
Reads supporting (aligned to +/- strand): ref base G (71/67); new base C (11/2); total (82/69)
Fisher's exact test for biased strand distribution p-value = 3.79e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01

Only 100 of 154 total aligned reads displayed.

GCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGAC  >  NC_002947/5754778‑5755073

gCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTaa                                                                                                                                                    >  2:5260865/1‑150 (MQ=255)

 ccGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAg                                                                                                                                                   >  1:3039362/1‑150 (MQ=255)

 ccGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAg                                                                                                                                                   <  1:3776042/150‑1 (MQ=255)

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

      gCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCt                                                                                                                                              <  2:2145284/150‑1 (MQ=255)

      gCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCt                                                                                                                                              <  2:4801304/150‑1 (MQ=255)

        cACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGa                                                                                                                                            <  1:4208836/150‑1 (MQ=255)

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

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

          cTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAAc                                                                                                                                          >  2:4676378/1‑150 (MQ=255)

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

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

             tACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGc                                                                                                                                       >  2:1806377/1‑150 (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)

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

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

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

                            tCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCg                                                                                                                        >  2:3953111/1‑150 (MQ=255)

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

                              gggCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGAc                                                                                                                      >  1:5124812/1‑150 (MQ=255)

                                    gCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGTTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTc                                                                                                                <  2:1786810/150‑1 (MQ=255)

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

                                      aCGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGt                                                                                                              <  2:3586627/150‑1 (MQ=255)

                                        gTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATTGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTAc                                                                                                            >  2:2668411/1‑150 (MQ=255)

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

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

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

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

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

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

                                                 tCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTg                                                                                                   >  1:2025948/1‑150 (MQ=255)

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

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

                                                   gATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGt                                                                                                 <  2:1353469/150‑1 (MQ=255)

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

                                                    aTCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTATGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTc                                                                                                >  1:4803874/1‑150 (MQ=255)

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

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

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

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

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

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

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

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

                                                            aGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCtt                                                                                        <  1:760135/150‑1 (MQ=255)

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

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

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

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

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

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

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

                                                                            ttCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTaggc                                                                        >  2:457558/1‑147 (MQ=255)

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

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

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

                                                                                  ggTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGcc                                                                  <  1:4676378/150‑1 (MQ=255)

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

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

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

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

                                                                                        ggCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCg                                                            >  1:3658575/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)

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

                                                                                           tATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCt                                                         >  2:504822/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)

                                                                                              gCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACGCGGCTCGc                                                      >  2:4065577/1‑150 (MQ=255)

                                                                                               ccGGCCTGTCGGAAAACTGCCTGTTCTTCATCGGCGGCGTCATCAAGCACGCCAAGGCTGTGAACGCATTCACCAACCCGTCGACCAACTCCTATAAGCGTCTGGTCCCAGGTTATGAAGCTCCGGTTCTGCTGGCCTATTCGGCTCgta                                                     >  1:890907/1‑148 (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:908393/1‑150 (MQ=255)

                                                                                                           acacCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTc                                                             <  1:2227319/130‑1 (MQ=255)

                                                                                                           acacCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTc                                                             >  2:2227319/1‑130 (MQ=255)

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

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

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

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

                                                                                                                            aTCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTAc                        >  2:1563053/1‑150 (MQ=255)

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

                                                                                                                              cggcggTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGt                      >  1:626464/1‑150 (MQ=255)

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

                                                                                                                                       catcaAGCACGGTAAGGCCCTGAACGTCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGccc             <  2:2851692/150‑1 (MQ=255)

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

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

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

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

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

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

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

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

GCCGAAGCCACTGTACGGCGACAACGGCTCGGGCATGCACGTGCACATGTCGATCTGGAAAGACGGCAAGAACACCTTCTCGGGTGAAGGCTATGCCGGCCTGTCCGACACCGCCCTGTACTTCATCGGCGGTATCATCAAGCACGGTAAGGCCCTGAACGGCTTCACCAACCCGTCGACCAACTCGTACAAGCGTCTGGTCCCAGGCTTCGAAGCCCCGGTAATGCTGGCCTACTCGGCTCGCAACCGTTCCGCCTCGATCCGTATTCCTTACGTTGGCAGCCCGAAAGCCCGAC  >  NC_002947/5754778‑5755073

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