BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_008463	6,070,098	G→C	pseudogene (567/1086 nt)	PA14_RS27735 →	oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_008463	6,070,098	0	G	C	99.8%	2231.7 / NA	574	pseudogene (567/1086 nt)	PA14_RS27735	oxidoreductase
Reads supporting (aligned to +/- strand): ref base G (0/0); major base C (278/295); minor base A (1/0); total (279/295)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 574 total aligned reads displayed.

AGGCCGCCCTGCAGGCCGGCAGCATCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTTCCGGGCTGTGGTTTGACCTTGGCCCGCAACCTGGTTGACCAAGGCACTGTGCCTGTTTGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTGGCGCTGCGCGCCAGCATGCT  >  NC_008463/6069953‑6070244

aGGCCGCCCTGCAGGCCGGCAGCATCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGAcc                                                                                                                                              <  1:5196336/150‑1 (MQ=255)

aGGCCGCCCTGCAGGCCGGCAGCATCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGAcc                                                                                                                                              <  1:6469401/150‑1 (MQ=255)

 ggCCGCCCTGCAGGCCGGCAGCATCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACCa                                                                                                                                             >  2:961846/1‑150 (MQ=255)

  gccgccCTGCAGGCCGGCAGCATCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACCag                                                                                                                                            <  1:5012169/150‑2 (MQ=255)

      cccTGCAGGCCGGCAGCATCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCAc                                                                                                                                       <  1:2843411/150‑1 (MQ=255)

      cccTGCAGGCCGGCAGCATCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCAc                                                                                                                                       <  2:3889067/150‑1 (MQ=255)

      cccTGCAGGCCGGCAGCATCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCAc                                                                                                                                       <  1:8200967/150‑1 (MQ=255)

       ccTGCAGGCCGGCAGCATCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACt                                                                                                                                      <  2:1351320/150‑1 (MQ=255)

        cTGCAGGCCGGCAGCATCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACtg                                                                                                                                     <  1:2198936/150‑1 (MQ=255)

         tGCAGGCCGGCAGCATCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCAGCAGGTCCGCGACCGCTGGCGCGATCACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACtgt                                                                                                                                    <  2:6426736/150‑1 (MQ=255)

          gCAGGCCGGCAGCATCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACtgtg                                                                                                                                   >  2:7370666/1‑150 (MQ=255)

           cAGGCCGGCAGCATCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGc                                                                                                                                  <  1:2096200/150‑1 (MQ=255)

           cAGGCCGGCAGCATCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGc                                                                                                                                  <  2:4033387/150‑1 (MQ=255)

            aGGCCGGCAGCATCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGcc                                                                                                                                 >  2:2153926/1‑150 (MQ=255)

            aGGCCGGCAGCATCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGcc                                                                                                                                 <  1:10602660/150‑1 (MQ=255)

             ggCCGGCAGCATCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCt                                                                                                                                >  2:8040638/1‑150 (MQ=255)

             ggCCGGCAGCATCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCt                                                                                                                                >  2:5263349/1‑150 (MQ=255)

             ggCCGGCAGCATCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCt                                                                                                                                <  2:426009/150‑1 (MQ=255)

             ggCCGGCAGCAACGGCCCGGTGCGCTACTTAGAGTCGGCCATCGAGCGCTTCCGTCCGATGGTCCGCGACCGCTGGCGCGACCACGACCTGCCCGGC‑TCCGGGATGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCt                                                                                                                                >  2:9304773/1‑150 (MQ=255)

              gCCGGCAGCATCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTg                                                                                                                               <  1:10052798/150‑1 (MQ=255)

               ccGGCAGCAACGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGt                                                                                                                              >  2:5944597/1‑150 (MQ=255)

                 ggCAGCATCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACCCGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGttc                                                                                                                            <  2:2657596/150‑2 (MQ=255)

                  gcagcaTCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGttcg                                                                                                                           >  2:768285/1‑148 (MQ=255)

                    agcaTCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCGCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGc                                                                                                                         <  2:8906950/150‑1 (MQ=255)

                     gcaTCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCa                                                                                                                        <  1:6984566/150‑1 (MQ=255)

                       aTCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATa                                                                                                                      <  1:8761888/150‑1 (MQ=255)

                              cGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGc                                                                                                               >  2:4473458/1‑150 (MQ=255)

                              cGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGc                                                                                                               >  1:10190964/1‑150 (MQ=255)

                              cGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGc                                                                                                               >  1:1502994/1‑150 (MQ=255)

                               ggTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCg                                                                                                              <  1:11026507/150‑1 (MQ=255)

                                 tGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTa                                                                                                            <  1:1756682/150‑1 (MQ=255)

                                  gcgcTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTAt                                                                                                           >  1:5846660/1‑150 (MQ=255)

                                       aCTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACg                                                                                                      <  1:7564363/150‑1 (MQ=255)

                                         ttCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCGTACCGCAGCGTATGCACGGa                                                                                                    >  2:10637433/1‑150 (MQ=255)

                                          tCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGac                                                                                                   <  2:10412401/150‑1 (MQ=255)

                                             aGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACAcc                                                                                                <  1:4934369/150‑1 (MQ=255)

                                                cGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTgc                                                                                             <  2:8749756/150‑1 (MQ=255)

                                                  gCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTgcgc                                                                                           >  1:10544065/1‑150 (MQ=255)

                                                      tCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGc                                                                                       <  1:1217932/150‑1 (MQ=255)

                                                      tCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGc                                                                                       <  1:4223400/150‑1 (MQ=255)

                                                       cGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCt                                                                                      >  2:1771336/1‑150 (MQ=255)

                                                        gAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTg                                                                                     >  1:7689937/1‑150 (MQ=255)

                                                         aGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTgc                                                                                    <  2:10838044/150‑1 (MQ=255)

                                                                 cgtccgCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGgc                                                                            >  1:9773780/1‑150 (MQ=255)

                                                                   tccgCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGgcgc                                                                          >  2:3171008/1‑150 (MQ=255)

                                                                   tccgCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGgcgc                                                                          >  1:4000925/1‑150 (MQ=255)

                                                                     cgCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGgcgcgc                                                                        >  2:4886804/1‑150 (MQ=255)

                                                                     cgCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGgcgcgc                                                                        <  2:5527668/150‑1 (MQ=255)

                                                                     cgCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGgcgcgc                                                                        >  2:8657834/1‑150 (MQ=255)

                                                                       cAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTg                                                                      >  1:388026/1‑150 (MQ=255)

                                                                       cAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTg                                                                      >  2:1129713/1‑150 (MQ=255)

                                                                         ggTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGAc                                                                    <  1:1846424/150‑1 (MQ=255)

                                                                          gTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGAcc                                                                   <  2:9858547/150‑1 (MQ=255)

                                                                          gTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGAcc                                                                   <  1:3514530/150‑1 (MQ=255)

                                                                            ccGCGACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACcga                                                                 <  1:10865516/150‑1 (MQ=255)

                                                                              gcgAGCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTTTGCACGGACACCTGCGCCGGCTGCGCGTGGGCGCGCTGAGcgtcg                                                               <  2:10081443/150‑1 (MQ=255)

                                                                               cgACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACcgacga                                                              >  1:3189595/1‑150 (MQ=255)

                                                                                gACCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCgacgac                                                             <  2:5007977/150‑1 (MQ=255)

                                                                                 aCCGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACt                                                            <  1:10412523/150‑1 (MQ=255)

                                                                                   cGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTgg                                                          <  1:5520788/150‑1 (MQ=255)

                                                                                   cGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTgg                                                          <  2:2344821/150‑1 (MQ=255)

                                                                                   cGCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTgg                                                          <  2:6665219/150‑1 (MQ=255)

                                                                                    gCTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGt                                                         <  1:7212898/150‑1 (MQ=255)

                                                                                     cTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGtt                                                        >  1:6764054/1‑150 (MQ=255)

                                                                                     cTGGCGCGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGtt                                                        >  1:3358726/1‑150 (MQ=255)

                                                                                        gcgcGAACACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGa                                                     <  2:570151/150‑1 (MQ=255)

                                                                                             aaCACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTcc                                                <  2:8325104/150‑1 (MQ=255)

                                                                                             aaCACGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTcc                                                <  1:4714353/150‑1 (MQ=255)

                                                                                              acacGACCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCt                                               <  1:4875479/150‑1 (MQ=255)

                                                                                                   aCCTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCg                                          >  2:7236808/1‑150 (MQ=255)

                                                                                                    ccTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCTCTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGa                                         <  1:2409400/150‑1 (MQ=255)

                                                                                                    ccTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGa                                         <  1:6059324/150‑1 (MQ=255)

                                                                                                    ccTGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGa                                         <  2:10244345/150‑1 (MQ=255)

                                                                                                      tGCCCGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACt                                       <  1:6632512/150‑1 (MQ=255)

                                                                                                        cccGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTAt                                     <  1:9243493/150‑1 (MQ=255)

                                                                                                         ccGGC‑TCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATc                                    <  1:7417977/150‑1 (MQ=255)

                                                                                                             gctCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGc                                 >  1:3739005/3‑150 (MQ=255)

                                                                                                             gctCCGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTAGCCGc                                 >  2:9453579/3‑150 (MQ=255)

                                                                                                                ccGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACt                              >  1:5393624/1‑150 (MQ=255)

                                                                                                                ccGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACt                              <  1:8474903/150‑1 (MQ=255)

                                                                                                                 cGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTg                             <  2:8463149/150‑1 (MQ=255)

                                                                                                                 cGGGCTGTGGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTg                             <  2:1421159/150‑1 (MQ=255)

                                                                                                                       gtgGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGt                       >  2:680762/1‑150 (MQ=255)

                                                                                                                        tgGTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTg                      >  2:3319460/1‑150 (MQ=255)

                                                                                                                         ggTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTgg                     >  2:7862950/1‑150 (MQ=255)

                                                                                                                         ggTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTgg                     >  1:9962065/1‑150 (MQ=255)

                                                                                                                          gTTCGACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTGgc                    <  1:7960086/150‑1 (MQ=255)

                                                                                                                             cgACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTGGCGCt                 >  1:6290245/2‑150 (MQ=255)

                                                                                                                             agACCTCGGCCCGC‑ACCTGCTCGACC‑AGGCCCTGTGCCTGTTCGGCATACCGCAGCGGATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACGGGTGCGAGGTCCTGCGGGACGATCCGCCCTGCGAGGTGGGGCt                 >  2:6228679/2‑150 (MQ=255)

                                                                                                                               aCCTCGGCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTGGCGCTgc               <  2:3516691/150‑1 (MQ=255)

                                                                                                                                 ctcggCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTGGCGCTgcgc             >  2:2522801/4‑150 (MQ=255)

                                                                                                                                     gCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTGGCGCTGCGCGCCa         <  2:753120/150‑1 (MQ=255)

                                                                                                                                     gCCCGC‑ACCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTGGCGCTGCGCGCCa         <  1:8942400/150‑1 (MQ=255)

                                                                                                                                        ccgcaCCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTGGCGCTGCGCGCCAGc       <  2:7149950/146‑1 (MQ=255)

                                                                                                                                         cgcaCCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACTGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTGGCGCTGCGCGCCAGCa      <  2:2855277/147‑1 (MQ=255)

                                                                                                                                         cgcaCCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTGGCGCTGCGCGCCAGCa      <  1:1984469/147‑1 (MQ=255)

                                                                                                                                          gcaCCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTGGCGCTGCGCGCCAGCAt     >  2:3716288/3‑150 (MQ=255)

                                                                                                                                            aCCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTGGCGCTGCGCGCCAGCATGc   <  1:5680338/150‑1 (MQ=255)

                                                                                                                                            aCCTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTGGCGCTGCGCGCCAGCATGc   >  2:1223027/1‑150 (MQ=255)

                                                                                                                                             ccTGCTCGACC‑AGGCACTGTGCCTGTTCGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTGGCGCTGCGCGCCAGCATGCt  >  2:5324366/1‑150 (MQ=255)

AGGCCGCCCTGCAGGCCGGCAGCATCGGCCCGGTGCGCTACTTCGAGTCGGCCATCGAGCGCTTCCGTCCGCAGGTCCGCGACCGCTGGCGCGAACACGACCTGCCCGGCTTCCGGGCTGTGGTTTGACCTTGGCCCGCAACCTGGTTGACCAAGGCACTGTGCCTGTTTGGCATACCGCAGCGTATGCACGGACACCTGCGCCGGCTGCGCGAGGGCGCGCTGACCGACGACTGGTTCGAGGTCCTGCTCGACTATCCGCACTGCCAGGTGGCGCTGCGCGCCAGCATGCT  >  NC_008463/6069953‑6070244

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