BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,787,624:1	+A	6.0%	coding (787/1260 nt)	envC →	activator of AmiB,C murein hydrolases, septal ring factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,787,624	1	.	A	6.0%	592.6 / 32.3	134	coding (787/1260 nt)	envC	activator of AmiB,C murein hydrolases, septal ring factor
Reads supporting (aligned to +/- strand): ref base . (64/62); new base A (4/4); total (68/66)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 141 total aligned reads displayed.

GGGCTGGAGTCTTCCATCCAGCAAGGTCAGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACGCTGCATCGCTATGGC  >  NC_000913/3787477‑3787770

      gAGTCTTCCACCCAGCAAGGTCAGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGcc                                                                                                                                              >  4:755234/1‑148 (MQ=255)

          cTTCCATCCAGCAAGGTCAGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCagaa                                                                                                                                          >  1:806041/1‑148 (MQ=255)

            tccatccaGCAAGGTCAGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCTGGCGG‑TCGCGAACCGcctg                                                                                                                                            >  4:827933/1‑142 (MQ=255)

            tccatccaGCAAGGTCAGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCagaaag                                                                                                                                        >  4:539736/1‑148 (MQ=255)

            tccatccaGCAAGGTCAGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGG‑TCGCG‑aa                                                                                                                                                  <  3:837324/137‑2 (MQ=255)

            tccatccaGCAAGGTCAGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATAGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCagaaag                                                                                                                                        >  1:916371/1‑148 (MQ=255)

              cacccaGCAAGGTCAGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCagaaagaa                                                                                                                                      <  3:755234/145‑1 (MQ=255)

                tccaGCAAGGTCAGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGcc                                                                                                                                              <  1:953152/138‑1 (MQ=255)

                tccaGCAAGGTCAGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGcc                                                                                                                                              >  2:953152/1‑138 (MQ=255)

                tccaGCAAGGTCAGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGc                                                                                                                                    >  1:901588/1‑148 (MQ=255)

                tccaGCAAGGTCAGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGG‑TCGCGAACCGcc                                                                                                                                              >  2:880203/1‑138 (MQ=255)

                 ccaGCAAGGTCAGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCg                                                                                                                                   <  4:546143/148‑1 (MQ=255)

                    gCAAGGTCAGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGG‑TCGCGAACCGcctg                                                                                                                                            >  4:921471/1‑134 (MQ=255)

                    gCAAGGTCAGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGG‑TCGCGAACCGcctg                                                                                                                                            <  3:921471/136‑3 (MQ=255)

                     cAAGGTCAGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCagaaagaa                                                                                                                                      >  4:1085778/1‑141 (MQ=255)

                         gTCAGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAaagaag                                                                                                                                     >  1:317494/1‑138 (MQ=255)

                         gTCAGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTATGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAaagaag                                                                                                                                     <  2:317494/138‑1 (MQ=255)

                          tCAGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAATGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCaaa                                                                                                                          <  2:916371/148‑1 (MQ=255)

                          tCAGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCaaa                                                                                                                          >  3:50584/1‑148 (MQ=255)

                            aGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAgg                                                                                                                        <  3:539736/148‑1 (MQ=255)

                             gCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAgg                                                                                                                        <  1:784439/147‑1 (MQ=255)

                              cAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCAGAAGCCGCGGCTAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGca                                                                                                                      <  2:806041/148‑1 (MQ=255)

                               aaCAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCagaaaga                                                                                                                                       <  1:942122/130‑1 (MQ=255)

                               aaCAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCagaaaga                                                                                                                                       >  2:942122/1‑130 (MQ=255)

                                aCAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCagaaaga                                                                                                                                       >  3:621580/1‑129 (MQ=255)

                                aCAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGcacc                                                                                                                    >  4:445012/1‑148 (MQ=255)

                                aCAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGcacc                                                                                                                    <  3:298574/148‑1 (MQ=255)

                                    ttGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCaa                                                                                                                           <  2:144562/137‑1 (MQ=255)

                                    ttGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCaa                                                                                                                           >  1:144562/1‑137 (MQ=255)

                                                      aaCGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAaccg                                                                                                        >  2:426055/1‑138 (MQ=255)

                                                      aaCGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGG‑TCGCGAACCGcc                                                                                                                                              <  1:993457/100‑1 (MQ=255)

                                                      aaCGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGG‑TCGCGAACCGcc                                                                                                                                              >  2:993457/1‑100 (MQ=255)

                                                       aCGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGcacca                                                                                                                   >  2:319484/1‑126 (MQ=255)

                                                             cccGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGcc                                                                                                                                              <  1:386300/93‑1 (MQ=255)

                                                             cccGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCaa                                                                                                                           <  3:696695/112‑1 (MQ=255)

                                                             cccGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTa                                                                                                               >  2:762807/1‑124 (MQ=255)

                                                              ccGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCg                                                                                      <  3:917629/148‑1 (MQ=255)

                                                               cGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGc                                                                                                                       <  4:908349/114‑1 (MQ=255)

                                                               cGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACaa                                                                                                            >  1:1012461/1‑125 (MQ=255)

                                                               cGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACaa                                                                                                            <  2:1012461/125‑1 (MQ=255)

                                                                 tCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTg                                                                                   >  2:742732/1‑148 (MQ=255)

                                                                  cTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGa                                                                                  <  1:742732/148‑1 (MQ=255)

                                                                  cTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGa                                                                                  >  2:185459/1‑148 (MQ=255)

                                                                     cGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCagaa                                                                                                                                          >  3:255717/1‑89 (MQ=255)

                                                                     cGTAACAGCATTGACCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGt                                                                               >  4:13021/1‑148 (MQ=255)

                                                                      gTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTc                                                                              <  4:50584/148‑1 (MQ=255)

                                                                      gTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTc                                                                              <  1:976985/148‑1 (MQ=255)

                                                                       tAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGa                                                                                                                                  >  4:167127/1‑95 (MQ=255)

                                                                       tAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGcacca                                                                                                                   <  2:390554/110‑1 (MQ=255)

                                                                       tAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGcacca                                                                                                                   >  1:390554/1‑110 (MQ=255)

                                                                       tAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCg                                                                                      <  2:575278/139‑1 (MQ=255)

                                                                       tAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCg                                                                                      >  1:575278/1‑139 (MQ=255)

                                                                         aCAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTcc                                                                             <  1:163191/146‑1 (MQ=255)

                                                                         aCAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTcc                                                                             >  2:163191/1‑146 (MQ=255)

                                                                         aCAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCg                                                                           >  2:239582/1‑148 (MQ=255)

                                                                          cAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCagaa                                                                                                                                          >  2:549871/1‑84 (MQ=255)

                                                                           aGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCa                                                                                                                                             <  2:133532/80‑1 (MQ=255)

                                                                           aGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCa                                                                                                                                             >  1:133532/1‑80 (MQ=255)

                                                                           aGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATc                                                                                       >  2:700088/1‑134 (MQ=255)

                                                                           aGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATc                                                                                       <  1:313966/134‑1 (MQ=255)

                                                                              aTTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCg                                                                      >  2:338265/1‑148 (MQ=255)

                                                                                      tGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTg                                                              <  1:239582/148‑1 (MQ=255)

                                                                                                 cggcgAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCagaaagaa                                                                                                                                      >  1:20433/1‑65 (MQ=255)

                                                                                                 cggcgAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCcgcg                                                   >  3:425813/1‑148 (MQ=255)

                                                                                                      aaaGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTgg                                                             <  2:495532/133‑1 (MQ=255)

                                                                                                      aaaGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTgg                                                             >  1:495532/1‑133 (MQ=255)

                                                                                                      aaaGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGt                                                  >  3:851515/1‑144 (MQ=255)

                                                                                                      aaaGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGt                                                  <  4:851515/144‑1 (MQ=255)

                                                                                                         gcgcGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAgg                                           <  4:425813/148‑1 (MQ=255)

                                                                                                           gcgTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCa                                         <  4:419428/148‑1 (MQ=255)

                                                                                                             gTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCAtt                                       <  1:338265/148‑1 (MQ=255)

                                                                                                               gCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGc                                                                                     <  1:1075645/100‑1 (MQ=255)

                                                                                                                cAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTgg                                                <  3:602455/136‑1 (MQ=255)

                                                                                                                cAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTgg                                                >  4:602455/1‑136 (MQ=255)

                                                                                                                cAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTc                                              >  1:479223/1‑138 (MQ=255)

                                                                                                                   aaCGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGCAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGcc                                 >  2:953615/1‑148 (MQ=255)

                                                                                                                   aaCGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTg                                                              <  1:288022/119‑1 (MQ=255)

                                                                                                                   aaCGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTg                                                              >  2:288022/1‑119 (MQ=255)

                                                                                                                   aaCGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCa                                             >  1:609643/1‑136 (MQ=255)

                                                                                                                   aaCGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCa                                             <  2:609643/136‑1 (MQ=255)

                                                                                                                   aaCGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGcc                                 <  3:13021/148‑1 (MQ=255)

                                                                                                                   aaCGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGcc                                 >  2:787024/1‑148 (MQ=255)

                                                                                                                    aCGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTg                                                              <  3:509665/118‑1 (MQ=255)

                                                                                                                    aCGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTg                                                 >  4:653317/1‑131 (MQ=255)

                                                                                                                              cGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGaa                                                                                                  <  3:116889/72‑1 (MQ=255)

                                                                                                                                tGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAg                                                                                                                         >  3:284114/1‑47 (MQ=255)

                                                                                                                                tGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAg                                                                                                                         <  4:284114/47‑1 (MQ=255)

                                                                                                                                tGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGt                              <  3:443076/138‑1 (MQ=255)

                                                                                                                                  aGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAAt                                                                                        <  3:1064612/78‑1 (MQ=255)

                                                                                                                                  aGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAAt                                                                                        >  4:1064612/1‑78 (MQ=255)

                                                                                                                                      ccAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACGCTg              >  4:96505/1‑148 (MQ=255)

                                                                                                                                        aGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGc                                                                                     >  3:979833/1‑75 (MQ=255)

                                                                                                                                        aGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACGCTGCa            >  4:690989/1‑148 (MQ=255)

                                                                                                                                          gCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTc                                              <  4:823606/112‑1 (MQ=255)

                                                                                                                                          gCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTc                                              >  3:823606/1‑112 (MQ=255)

                                                                                                                                              ttCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTcgc                          <  1:72429/128‑1 (MQ=255)

                                                                                                                                               tCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCggtgg                                                                  >  1:922908/1‑87 (MQ=255)

                                                                                                                                               tCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCggtgg                                                                  <  2:922908/87‑1 (MQ=255)

                                                                                                                                                  cg‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGc                                  <  1:554329/116‑1 (MQ=255)

                                                                                                                                                  cg‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACGCTGCATCGCTATGGc  >  2:1096493/1‑148 (MQ=255)

GGGCTGGAGTCTTCCATCCAGCAAGGTCAGCAACAGTTGAGCGAGCTGCGCGCCAACGAATCCCGTCTGCGTAACAGCATTGCCCGTGCGGAAGCCGCGGCGAAAGCGCGTGCAGAACGAGAAGCACGTGAGGCCCAGGCGGTTCGCG‑ACCGCCAGAAAGAAGCGACGCGCAAAGGCACCACCTACAAACCGACCGAAAGCGAAAAATCGCTGATGTCCCGTACCGGTGGTCTGGGCGCGCCGCGTGGTCAGGCATTCTGGCCGGTTCGCGGGCCGACGCTGCATCGCTATGGC  >  NC_000913/3787477‑3787770

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