BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_008463	2,001,811	A→C	Q66P (CAG→CCG)	pgl →	6‑phosphogluconolactonase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_008463	2,001,811	0	A	C	100.0%	1522.6 / NA	401	Q66P (CAG→CCG)	pgl	6‑phosphogluconolactonase
Reads supporting (aligned to +/- strand): ref base A (0/0); new base C (172/229); total (172/229)

Only 100 of 401 total aligned reads displayed.

CAGCGCCGCCGAGCAGGCCCGCGGCCTGGCCGCCGAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGATCAGGTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAGGCGGCGAAGGCGCGCTTCATCGGCCTCTACCAG  >  NC_008463/2001668‑2001953

 aGCGCCGCCGAGCAGGCCCGCGGCCTGGCCGCCGAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGAc                                                                                                                                         <  1:10529041/150‑1 (MQ=255)

 aGCGCCGCCGAGCAGGCCCGCGGCCTGGCCGCCGAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCGCGAAGCCTTGAGCGAGGAGCCGCTGGAc                                                                                                                                         >  1:8651834/1‑150 (MQ=255)

  gcgcCGCCGAGCAGGCCCGCGGCCTGGCCGCCGAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACt                                                                                                                                        >  1:6803281/1‑150 (MQ=255)

  gcgcCGCCGAGCAGGCCCGCGGCCTGGCCGCCGAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACt                                                                                                                                        >  1:3110111/1‑150 (MQ=255)

      cgccgAGCAGGCCCGCGGCCTGGCCGCCGAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTc                                                                                                                                    <  1:4892075/150‑1 (MQ=255)

      cgccgAGCAGGCCCGCGGCCTGGCCGCCGAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTc                                                                                                                                    <  2:4193278/150‑1 (MQ=255)

      cgccgAGCAGGCCCGCGGCCTGGCCGCCGAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTc                                                                                                                                    >  2:10586027/1‑150 (MQ=255)

       ccctaGCAGGCCCGCGTCCTCGCCGGCGAGGTCCCCGGCCGGTTCCGCTCTGCCATGGCCGTGCAGCGCCAGGCGCAGCTGGTGGTGTCCGGTGGGCGGAGTCCGGTGGCCTTCCTCGAACCCTTGAGCGAGGAGCCGCTGGACTGGTcg                                                                                                                                   <  2:7272522/146‑1 (MQ=255)

          gAGCAGGCCCGCGGCCTGGCCGCCGAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCgg                                                                                                                                <  2:1653729/150‑1 (MQ=255)

             cAGGCCCGCGGCCTGGCCGCCGAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATc                                                                                                                             >  2:307909/1‑150 (MQ=255)

              aGGCCCGCGGCCTGGCCGCCGAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCa                                                                                                                            >  1:5809589/1‑150 (MQ=255)

                 cccGCGGCCTGGCCGCCGAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCAccc                                                                                                                         <  1:6478214/150‑1 (MQ=255)

                  ccGCGGCCTGGCCGCCGAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCt                                                                                                                        >  2:4928989/1‑150 (MQ=255)

                   cgcgGCCTGGCCGCCGAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCGGCTGGTGGTGTCCGGTGGGCGCAGTCCGGGGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTc                                                                                                                       <  1:1874114/150‑1 (MQ=255)

                      ggCCTGGCCGCCGAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGc                                                                                                                    <  2:9585865/150‑1 (MQ=255)

                       gCCTGGCCGCCGAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGcc                                                                                                                   <  2:6902911/150‑1 (MQ=255)

                        ccTGGCCGCCGAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCt                                                                                                                  <  2:666799/150‑1 (MQ=255)

                         cTGGCCGCCGAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTg                                                                                                                 <  1:8408727/150‑1 (MQ=255)

                          tGGCCGCCGAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTgg                                                                                                                >  2:4421322/1‑150 (MQ=255)

                          tGGCCGCCGAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCGTCCTCGAAGCCTTGTGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCGCGGCCTgg                                                                                                                >  1:9294260/1‑150 (MQ=255)

                           ggCCGCCGAGGTCTCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGACTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGc                                                                                                               >  1:442595/1‑150 (MQ=255)

                                 cgAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCcgacga                                                                                                         <  1:8098484/150‑1 (MQ=255)

                                      tCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCt                                                                                                    <  2:3837566/150‑1 (MQ=255)

                                         ccggccggTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTggg                                                                                                 >  2:5078877/1‑150 (MQ=255)

                                          cggccggTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGt                                                                                                <  2:1375048/150‑1 (MQ=255)

                                           ggccggTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTg                                                                                               >  2:8312869/1‑150 (MQ=255)

                                            gccggTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCCGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGGCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGc                                                                                              >  2:616113/1‑150 (MQ=255)

                                            gccggTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGc                                                                                              <  2:5383690/150‑1 (MQ=255)

                                             ccggTGGCGCTCGGCGCAGGCCGAGCAGGGCCAGGCGCTGCCGGTGCTGGCCGGTGGGCGCAGGCCGGTGGCCTCCCTCGAACCCGTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCACAGCCTGGCCGACGAGCGCTGGGTGcc                                                                                             <  2:8893406/150‑1 (MQ=255)

                                                gTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGa                                                                                          <  1:278752/150‑1 (MQ=255)

                                                   gcgcTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTc                                                                                       >  2:6729921/1‑150 (MQ=255)

                                                    cgcTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCg                                                                                      >  1:4568428/1‑150 (MQ=255)

                                                    cgcTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCg                                                                                      <  1:1615098/150‑1 (MQ=255)

                                                     gcTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGc                                                                                     >  1:9563830/1‑150 (MQ=255)

                                                     gcTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGc                                                                                     <  1:7267881/150‑1 (MQ=255)

                                                          gcgcTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGcc                                                                                <  1:1111274/150‑1 (MQ=255)

                                                          gcgcTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGcc                                                                                >  2:10851306/1‑150 (MQ=255)

                                                            gcTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGa                                                                              >  1:10929426/1‑150 (MQ=255)

                                                            gcTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGa                                                                              <  1:10075455/150‑1 (MQ=255)

                                                            gcTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCGTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGa                                                                              <  2:555014/150‑1 (MQ=255)

                                                             cTGGCCGAGCAGGGCTAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGAt                                                                             <  2:3619943/150‑1 (MQ=255)

                                                             cTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGAt                                                                             <  2:3379826/150‑1 (MQ=255)

                                                              gggCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATa                                                                            <  1:7426595/149‑1 (MQ=255)

                                                               ggCCGAGCCGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAg                                                                           <  2:2692699/150‑1 (MQ=255)

                                                               ggCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAg                                                                           <  2:7756072/150‑1 (MQ=255)

                                                                 ccGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCa                                                                         <  1:7212771/150‑1 (MQ=255)

                                                                 ccGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCa                                                                         >  1:7081948/1‑150 (MQ=255)

                                                                 ccGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCGTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCa                                                                         <  2:7387406/150‑1 (MQ=255)

                                                                 ccGAGCAGGGCCAGGCGCTGCGGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCa                                                                         <  1:5582485/150‑1 (MQ=255)

                                                                  cGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCaa                                                                        <  1:10367392/150‑1 (MQ=255)

                                                                  cGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCaa                                                                        <  1:478076/150‑1 (MQ=255)

                                                                   gAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAAc                                                                       >  1:5114123/1‑150 (MQ=255)

                                                                    aGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACg                                                                      <  2:9752851/150‑1 (MQ=255)

                                                                    aGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACg                                                                      <  1:9606798/150‑1 (MQ=255)

                                                                     gCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCGTCCTCGAAGCCTTGAGCGAGGAGCCGCGGGACGGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGc                                                                     <  2:2185019/150‑1 (MQ=255)

                                                                       aGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCTAAGCCTTGAGCGAGTAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGACGGAGTCGCATGCCGATAGCAACGCCg                                                                   <  2:4717125/150‑1 (MQ=255)

                                                                       aGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCg                                                                   <  1:6191107/150‑1 (MQ=255)

                                                                        gggCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCgg                                                                  <  2:9172240/150‑1 (MQ=255)

                                                                        gggCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCgg                                                                  <  2:9210043/150‑1 (MQ=255)

                                                                           ccAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCt                                                               <  1:7911929/150‑1 (MQ=255)

                                                                            cAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTg                                                              >  1:8267568/1‑150 (MQ=255)

                                                                               gcgcTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGtt                                                           <  2:7085100/150‑1 (MQ=255)

                                                                               gcgcTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGtt                                                           >  2:579425/1‑150 (MQ=255)

                                                                                        gtggtgTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCAc                                                  <  1:9534821/150‑1 (MQ=255)

                                                                                        gtggtgTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCAc                                                  >  1:2313666/1‑150 (MQ=255)

                                                                                         tggtgTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCAcc                                                 <  2:6481036/150‑1 (MQ=255)

                                                                                             gtCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCt                                             <  2:1861227/150‑1 (MQ=255)

                                                                                             gtCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCt                                             <  2:402139/150‑1 (MQ=255)

                                                                                               ccGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTcc                                           <  1:265436/150‑1 (MQ=255)

                                                                                               ccGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTcc                                           >  1:9011953/1‑150 (MQ=255)

                                                                                                  gTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTg                                        >  1:4225334/1‑150 (MQ=255)

                                                                                                   tGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTgg                                       <  2:8169425/150‑1 (MQ=255)

                                                                                                    gggCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGc                                      >  2:2587169/1‑150 (MQ=255)

                                                                                                     ggCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCg                                     <  1:1684648/150‑1 (MQ=255)

                                                                                                       cgcAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAg                                   <  1:3396262/150‑1 (MQ=255)

                                                                                                       cgcAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAg                                   <  1:3508494/150‑1 (MQ=255)

                                                                                                        gcAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAgg                                  >  2:2314567/1‑150 (MQ=255)

                                                                                                          aGTCCGGTGGCCTTCCTCGAAGCCTAGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCCGCCTGGCCGACGCGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAggcg                                <  2:3140168/150‑1 (MQ=255)

                                                                                                           gTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAggcgg                               >  1:2018390/1‑150 (MQ=255)

                                                                                                              cGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAGGCGGCGa                            >  2:1125779/1‑150 (MQ=255)

                                                                                                                 tGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCGGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAGGCGGCGAAgg                         >  1:2418622/1‑150 (MQ=255)

                                                                                                                  ggCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAGGCGGCGAAGgc                        >  2:5243129/1‑150 (MQ=255)

                                                                                                                   gCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAGGCGGCGAAGgcg                       <  1:140627/150‑1 (MQ=255)

                                                                                                                        ccTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAGGCGGCGAAGGCGCGCtt                  <  1:10193193/150‑1 (MQ=255)

                                                                                                                        ccTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAGGCGGCGAAGGCGCGCtt                  >  2:1038295/1‑150 (MQ=255)

                                                                                                                        ccTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAGGCGGCGAAGGCGCGCtt                  <  1:6912420/150‑1 (MQ=255)

                                                                                                                          tCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAGGCGGCGAAGGCGCGCTTCa                <  2:8548144/150‑1 (MQ=255)

                                                                                                                           cGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAGGCGGCGAAGGCGCGCTTCAt               >  1:9795615/1‑150 (MQ=255)

                                                                                                                           cGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAGGCGGCGAAGGCGCGCTTCAt               <  1:3988820/150‑1 (MQ=255)

                                                                                                                             aaGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAGGCGGCGAAGGCGCGCTTCATCg             >  1:1081765/1‑150 (MQ=255)

                                                                                                                              aGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAGGCGGCGAAGGCGCGCTTCATCgg            >  1:8966545/1‑150 (MQ=255)

                                                                                                                                 cTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAGGCGGCGAAGGCGCGCTTCATCGGCct         >  1:5411904/1‑150 (MQ=255)

                                                                                                                                  ttGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAGGCGGCGAAGGCGCGCTTCATCGGCctc        >  2:3295817/1‑150 (MQ=255)

                                                                                                                                  ttGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAGGCGGCGAAGGCGCGCTTCATCGGCctc        <  2:10474146/150‑1 (MQ=255)

                                                                                                                                    gAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAGGCGGCGAAGGCGCGCTTCATCGGCCTCTa      <  1:1246441/150‑1 (MQ=255)

                                                                                                                                     aGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAGGCGGCGAAGGCGCGCTTCATCGGCCTCTAc     >  2:5657755/1‑150 (MQ=255)

                                                                                                                                     aGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAGGCGGCGAAGGCGCGCTTCATCGGCCTCTAc     >  2:3927082/1‑150 (MQ=255)

                                                                                                                                        gaggagCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAGGCGGCGAAGGCGCGCTTCATCGGCCTCTACCAg  >  2:10963454/1‑150 (MQ=255)

                                                                                                                                        gaggagCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAGGCGGCGAAGGCGCGCTTCATCGGCCTCTACCAg  <  1:9998460/150‑1 (MQ=255)

                                                                                                                                        gaggagCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAGGCGGCGAAGGCGCGCTTCATCGGCCTCTACCAg  >  1:8836494/1‑150 (MQ=255)

CAGCGCCGCCGAGCAGGCCCGCGGCCTGGCCGCCGAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGATCAGGTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAGGCGGCGAAGGCGCGCTTCATCGGCCTCTACCAG  >  NC_008463/2001668‑2001953

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