BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_008463	5,386,198	G→A	A289V (GCG→GTG)	proB ←	glutamate 5‑kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_008463	5,386,198	0	G	A	99.1%	2340.6 / 2.2	567	A289V (GCG→GTG)	proB	glutamate 5‑kinase
Reads supporting (aligned to +/- strand): ref base G (0/5); new base A (314/247); total (314/253)
Fisher's exact test for biased strand distribution p-value = 1.71e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 579 total aligned reads displayed.

CAGCGCACTGTAGTTGACCAGGCCGCGCGCCACTTCGCGCCCGCCCTGGTCCACACAGACCACCATCTCGCCGCGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCGCCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGATCAGGGGTCAACAGGGTGCCGAGGCGTTCGCCGGCGCGCAGGCGGTCGAGGACACGCTCGATACGCCCGC  >  NC_008463/5386050‑5386348

cAGCGCACTGTAGTTGACCAGGCCGCGCGCCACTTCGCGCCCGCCCTGGTCCACACAGACCACCATCTCGCCGCGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCaccacc                                                                                                                                                      <  2:10001535/151‑1 (MQ=255)

 aGCGCACTGTAGTTGACCAGGCCGCGCGCCACTTCGCGCCCGCCCTGGTCCACACAGACCACCATCTCGCCGCGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCAccc                                                                                                                                                     <  1:9110917/151‑1 (MQ=255)

   cgcACTGTAGTTGACCCGGCCGCGCGCCACTTCGCGCCCGCCCTGGTCCACACAGACCACCATCTCGCCGCGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCa                                                                                                                                                   <  1:2339539/151‑1 (MQ=255)

    gcACTGTAGTTGACCAGGCCGCGCGCCACTTCGCGCCCGCCCTGGTCCACACAGACCACCATCTCGCCGCGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCag                                                                                                                                                  >  2:7537495/1‑150 (MQ=255)

     cACTGTAGTTGACCAGGCCGCGCGCCACTTCGCGCCCGCCCTGGTCCACACAGACCACCATCTCGCCGCGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATc                                                                                                                                                 <  2:4158856/151‑1 (MQ=255)

     cACTGTAGTTGACCAGGCCGCGCGCCACTTCGCGCCCGCCCTGGTCCACACAGACCACCATCTCGCCGCGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATc                                                                                                                                                 <  1:10034251/151‑1 (MQ=255)

         gTAGTTGACCAGGCCGCGCGCCACTTCGCGCCCGCCCTGGTCCACACAGACCACCATCTCGCCGCGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTcc                                                                                                                                             <  1:6985166/151‑1 (MQ=255)

          tAGTTGACCAGGCCGCGCGCCACTTCGCGCCCGCCCTGGTCCACACAGACCACCATCTCGCCGCGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCa                                                                                                                                            <  1:10739528/151‑1 (MQ=255)

          tAGGTGCCCAGGCCGCGCGCCACTTCGCGCCCGCCCTGGTCCACCCAGACCACCATCTCGCCGCGCCGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCa                                                                                                                                            <  2:1728448/151‑1 (MQ=255)

           aGTTGACCAGGCCGCGCGCCACTTCGCGCCCGCCCTGGTCCACACAGACCACCATCTCGCCGCGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAg                                                                                                                                           >  2:3825768/1‑151 (MQ=255)

               gACCAGGCCGCGCGCCACTTCGCGCCCGCCCTGGTCCACACAGACCACCATCTCGCCGCGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCAcc                                                                                                                                       >  2:2614517/1‑151 (MQ=255)

               gACCAGGCCGCGCGCCACTTCGCGCCCGCCCTGGTCCACACAGACCACCATCTCGCCGCGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCAcc                                                                                                                                       >  1:9500428/1‑151 (MQ=255)

               gACCAGGCCGCGCGCCACTTCGCGCCCGCCCTGGTCCACACAGACCACCATCTCGCCGCGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCAcc                                                                                                                                       <  2:5585045/151‑1 (MQ=255)

               gACCAGGCCGCGCGCCACTTCGCGCCCGCCCTGGTCCACACAGACCACCATCTCGCCGCGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTGGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCAcc                                                                                                                                       >  2:2793606/1‑151 (MQ=255)

                aCCAGGCCGCGCGCCACGTCGCGCCCGCCCTGGTCCACACAGACCACCATCTCGCCGCGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCGTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCAccc                                                                                                                                      >  1:4631242/1‑150 (MQ=255)

                   aGGCCGCGCGCCACTTCGCGCCCGCCCTGGTCCACACAGACCACCATCTCGCCGCGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAggg                                                                                                                                   >  1:5182514/1‑151 (MQ=255)

                   aGGCCGCGCGCCACTTCGCGCCCGCCCTGGTCCACACAGACCACCATCTCGCCGCGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAggg                                                                                                                                   >  2:7244057/1‑151 (MQ=255)

                    ggCCGCGCGCCACTTCGCGCCCGCCCTGGTCCACACAGACCACCATCTCGCCGCGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGt                                                                                                                                  >  1:1836838/1‑151 (MQ=255)

                         cgcgcCACTTCGCGCCCGCCCTGGTCCACACAGACCACCATCTCGCCGCGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCcgc                                                                                                                             >  1:1589973/1‑151 (MQ=255)

                         cgcgcCACTTCGCGCCCGCCCTGGTCCACACAGACCACCATCTCGCCGCGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCcgc                                                                                                                             <  1:419061/151‑1 (MQ=255)

                           cgcCACTTCGCGCCCGCCCTGGTCCACACAGACCACCATCTCGCCGCGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCgcgc                                                                                                                           >  1:8063328/1‑151 (MQ=255)

                            gcCACTTCGCGCCCGCCCTGGTCCACACAGACCACCATCTCGCCGCGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCgcgcc                                                                                                                          >  2:7794772/1‑150 (MQ=255)

                            gcCACTTCGCGCCCGCCCTGGTCCACACAGACCACCATCTCGCCGCGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCa                                                                                                                          >  2:35719/1‑151 (MQ=255)

                            gcCACTTCGCGCCCGCCCTGGTCCACACAGACCACCATCTCGCCGCGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCa                                                                                                                          >  2:3012402/1‑151 (MQ=255)

                             ccACTTCGCGCCCGCCCTGGTCCACACAGACCACCATCTCGCCGCGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCAt                                                                                                                         <  2:6638866/151‑1 (MQ=255)

                               aCTTCGCGCCCGCCCTGGTCCACACAGACCACCATCTCGCCGCGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCt                                                                                                                       >  1:9252497/1‑151 (MQ=255)

                               aCTTCGCGCCCGCCCTGGTCCACACAGACCACCATCTCGCCGCGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCt                                                                                                                       >  2:5099890/1‑151 (MQ=255)

                                  tCGCGCCCGCCCTGGTCCACACAGACCACCATCTCGCCGCGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCa                                                                                                                    <  2:396012/151‑1 (MQ=255)

                                    gcgcCCGCCCTGGTCCACACAGACCACCATCTCGCCGCGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAg                                                                                                                  >  2:4567695/1‑151 (MQ=255)

                                       cccgcccTGGTCCACACAGACCACCATCTCGCCGCGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTgg                                                                                                               <  2:4881215/151‑1 (MQ=255)

                                          gcccTGGTCCACACAGACCACCATCTCGCCGCGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTggccg                                                                                                            <  1:1049288/151‑1 (MQ=255)

                                            ccTGGTCCACACAGACCACCATCTCGCCGCGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTggccggc                                                                                                          >  1:8716007/1‑151 (MQ=255)

                                            ccTGGTCCACACAGACCACCATCTCGCCGCGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTggccggc                                                                                                          >  2:8152719/1‑151 (MQ=255)

                                             cTGGTCCACACAGACCACCATCTCGCCGCGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTggccggcc                                                                                                         >  1:10250266/1‑151 (MQ=255)

                                             cTGGTCCACACAGACCACCATCTCGCCGCGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTggccggcc                                                                                                         >  1:4808553/1‑151 (MQ=255)

                                               ggTCCACACAGACCACCATCTCGCCGCGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGgccag                                                                                                       >  2:9230772/1‑151 (MQ=255)

                                                      acaGACCACCATCTCGCCGCGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGc                                                                                                <  2:5962158/151‑1 (MQ=255)

                                                            caccaTCTCGCCGCGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCACGTGGCCGGCCAGCCACTGCTTgcgc                                                                                          >  2:1851407/1‑151 (MQ=255)

                                                             accaTCTCGCCGCGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTgcgcg                                                                                         <  2:6989592/151‑1 (MQ=255)

                                                               caTCTCGCCGCGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGcgcc                                                                                       >  1:6671028/1‑151 (MQ=255)

                                                                 tctcGCCGCGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGcgccgc                                                                                     >  1:4426957/1‑151 (MQ=255)

                                                                  ctcGCCGCGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCgccgcc                                                                                    >  2:8365494/1‑151 (MQ=255)

                                                                  ctcGCCGCGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCgccgcc                                                                                    >  1:7140334/1‑151 (MQ=255)

                                                                   tcGCCGCGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCt                                                                                   <  1:9620208/151‑1 (MQ=255)

                                                                    cgccgcGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCtt                                                                                  <  2:10797707/151‑1 (MQ=255)

                                                                    cgccgcGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCtt                                                                                  >  2:9839081/1‑151 (MQ=255)

                                                                      ccgcGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCTCCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGc                                                                                <  2:4217335/151‑1 (MQ=255)

                                                                       cgcgACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCg                                                                               >  1:857880/1‑151 (MQ=255)

                                                                       cgcgACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCg                                                                               >  2:3243453/1‑151 (MQ=255)

                                                                        gcgACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGa                                                                              >  1:2431957/1‑151 (MQ=255)

                                                                         cgACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGAc                                                                             >  2:6890518/1‑151 (MQ=255)

                                                                         cgACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGAc                                                                             >  2:1489296/1‑151 (MQ=255)

                                                                                aGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGATc                                                                      >  1:543663/1‑151 (MQ=255)

                                                                                aGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGATc                                                                      >  1:32680/1‑151 (MQ=255)

                                                                                 gCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGATCa                                                                     <  2:10073318/151‑1 (MQ=255)

                                                                                    gCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGATCAggg                                                                  >  2:3379989/1‑151 (MQ=255)

                                                                                      ccTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGATCAGGGGt                                                                >  1:2817735/1‑151 (MQ=255)

                                                                                        tGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGATCAGGGGTCa                                                              >  1:3600813/1‑151 (MQ=255)

                                                                                        tGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGATCAGGGGTCa                                                              >  2:7730566/1‑151 (MQ=255)

                                                                                         ggACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGATCAGGGGTCaa                                                             >  1:6077636/1‑151 (MQ=255)

                                                                                          gACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGATCAGGGGTCAAc                                                            >  2:7977190/1‑151 (MQ=255)

                                                                                          gACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGACCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGACAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGTTAAGGGGTCAAc                                                            >  2:12211/1‑151 (MQ=255)

                                                                                             cGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGATCAGGGGTCAACAgg                                                         >  2:2103480/1‑151 (MQ=255)

                                                                                              gCCTTGTCCCCAACCCGCAACCGGCTCTTGAGAGCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGATCAGGGGTCAACAggg                                                        <  2:6257803/151‑1 (MQ=255)

                                                                                               ccTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCTAGTGGCCGGCCTGCCACTGCTTGCGCGCCGCCTTGCGACTGCGAGCAGGGGTCAACAGGGt                                                       >  1:4000797/1‑151 (MQ=255)

                                                                                                cTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGATCAGGGGTCAACAGGGTg                                                      <  2:7400371/151‑1 (MQ=255)

                                                                                                  tGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGATCAGGGGTCAACAGGGTGcc                                                    <  2:9446764/151‑1 (MQ=255)

                                                                                                  tGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGATCAGGGGTCAACAGGGTGcc                                                    <  1:7857915/151‑1 (MQ=255)

                                                                                                    aCCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGATCAGGGGTCAACAGGGTGCCGa                                                  >  2:5404550/1‑151 (MQ=255)

                                                                                                     ccccAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGATCAGGGGTCAACAGGGTGCCGAg                                                 >  1:1883618/1‑151 (MQ=255)

                                                                                                     ccccAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGATCAGGGGTCAACAGGGTGCCGAg                                                 <  2:6400569/151‑1 (MQ=255)

                                                                                                       ccAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGATCAGGGGTCAACAGGGTGCCGAGGc                                               >  1:7424226/1‑151 (MQ=255)

                                                                                                          aCCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGATCAGGGGTCAACAGGGTGCCGAGGCGtt                                            >  2:4520758/1‑151 (MQ=255)

                                                                                                          aCCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGATCAGGGGTCAACAGGGTGCCGAGGCGtt                                            >  1:1660671/1‑151 (MQ=255)

                                                                                                          aCCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGATCAGGGGTCAACAGGGTGCCGAGGCGtg                                            >  2:9131620/1‑150 (MQ=255)

                                                                                                              gCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGATCAGGGGTCAACAGGGTGCCGAGGCGTTCGcc                                        >  1:10198649/1‑151 (MQ=255)

                                                                                                               cAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGATCAGGGGTCAACAGGGTGCCGAGGCGTTCGCCg                                       >  1:6536334/1‑151 (MQ=255)

                                                                                                               cAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGATCAGGGGTCAACAGGGTGCCGAGGCGTTCGCCg                                       <  1:120304/151‑1 (MQ=255)

                                                                                                                   aGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGATCAGGGGTCAACAGGGTGCCGAGGCGTTCGCCGgcgc                                   >  1:8252868/1‑151 (MQ=255)

                                                                                                                   aGGCTCTTGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGATCAGGGGTCAACAGGGTGCCGAGGCGTTCGCCGgcgc                                   >  2:2350399/1‑151 (MQ=255)

                                                                                                                         ttGTGATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGATCAGGGGTCAACAGGGTGCCGAGGCGTTCGCCGGCGCGCAGGc                             <  1:2755491/151‑1 (MQ=255)

                                                                                                                           gtgATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGATCAGGGGTCAACAGGGTGCCGAGGCGTTCGCCGGCGCGCAGGCgg                           <  2:5713884/151‑1 (MQ=255)

                                                                                                                            tgATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGATCAGGGGTCAACAGGGTGCCGAGGCGTTCGCCGGCGCGCAGGCGGt                          >  2:8789883/1‑151 (MQ=255)

                                                                                                                            tgATCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGATCAGGGGTCAACAGGGTGCCGAGGCGTTCGCCGGCGCGCAGGCGGt                          <  2:6717066/151‑1 (MQ=255)

                                                                                                                              aTCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGATCAGGGGTCAACAGGGTGCCGAGGCGTTCGCCGGCGCGCAGGCGGTCg                        <  1:3756959/151‑1 (MQ=255)

                                                                                                                               tCCTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGATCAGGGGTCAACAGGGTGCCGAGGCGTTCGCCGGCGCGCAGGCGGTCGa                       <  1:8468454/151‑1 (MQ=255)

                                                                                                                                ccTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGATCAGGGGTCAACAGGGTGCCGAGGCGTTCGCCGGCGCGCAGGCGGTCGAg                      <  1:6780974/151‑1 (MQ=255)

                                                                                                                                 cTGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGATCAGGGGTCAACAGGGTGCCGAGGCGTTCGCCGGCGCGCAGGCGGTCGAgg                     >  1:9813407/1‑151 (MQ=255)

                                                                                                                                  tGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGATCCGGGGTCAACAGGGTGCCGAGGCGTTCGCCGGCGCGCAGGCGGTCGAGGa                    >  1:7538335/1‑151 (MQ=255)

                                                                                                                                  tGGCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGATCAGGGGTCAACAGGGTGCCGAGGCGTTCGCCGGCGCGCAGGCGGTCGAGGa                    <  1:2208097/151‑1 (MQ=255)

                                                                                                                                   ggCTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGATCAGGGGTCAACAGGGTGCCGAGGCGTTCGCCGGCGCGCAGGCGGTCGAGGac                   <  1:2610129/151‑1 (MQ=255)

                                                                                                                                     cTCACCGCCTTCACCACCCCCTCGTCCAGCACCCGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGATCAGGGGTCAACAGGGTGCCGAGGCGTTCGCCGGCGCGCAGGCGGTCGAGGacac                 <  2:1047156/151‑1 (MQ=255)

                                                                                                                                     cTCACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGATCAGGGGTCAACAGGGTGCCGAGGCGTTCGCCGGCGCGCAGGCGGTCGAGGacac                 <  2:8886970/151‑1 (MQ=255)

                                                                                                                                       cACCGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGATCAGGGGTCAACAGGGTGCCGAGGCGTTCGCCGGCGCGCAGGCGGTCGAGGACACGc               <  2:10627627/151‑1 (MQ=255)

                                                                                                                                         ccGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGATCAGGGGTCAACAGGGTGCCGAGGCGTTCGCCGGCGCGCAGGCGGTCGAGGACACGCTc             <  2:8626699/151‑1 (MQ=255)

                                                                                                                                         ccGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACAGCGATCAGGGGTCAACAGGGTGCCGAGGCGTTCGCCGGCGCGCAGGCGGTCGAGGACACGCTc             <  1:5659072/151‑1 (MQ=255)

                                                                                                                                          cGCCTTCACCACCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGATCAGGGGTCAACAGGGTGCCGAGGCGTTCGCCGGCGCGCAGGCGGTCGAGGACACGCTCg            >  1:5838434/1‑151 (MQ=255)

                                                                                                                                                caccacCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGATCAGGGGTCAACAGGGTGCCGAGGCGTTCGCCGGCGCGCAGGCGGTCGAGGACACGCTCGATAcgc      >  1:6453394/1‑151 (MQ=255)

                                                                                                                                                caccacCCCATCGTCCAGCACCAGGGGGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGATCAGGGGTCAACAGGGTGCCGAGGCGTTCGCCGGCGCGCAGGCGGTCGAGGACACGCTCGATAcgc      <  2:10270676/151‑1 (MQ=255)

                                                                                                                                                 accaccCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGATCAGGGGTCAACAGGGTGCCGAGGCGTTCGCCGGCGCGCAGGCGGTCGAGGACACGCTCGATAcgcc     <  1:1329512/151‑1 (MQ=255)

CAGCGCACTGTAGTTGACCAGGCCGCGCGCCACTTCGCGCCCGCCCTGGTCCACACAGACCACCATCTCGCCGCGACGGAAGCTGCCCTGGACCGCCTTGACCCCAACCGGCAACAGGCTCTTGTGATCCTGGCTCACCGCCTTCACCGCCCCATCGTCCAGCACCAGGGTGCCGCGCATCTGCAAGTGGCCGGCCAGCCACTGCTTGCGCGCCGCCTTGCGACTGCGATCAGGGGTCAACAGGGTGCCGAGGCGTTCGCCGGCGCGCAGGCGGTCGAGGACACGCTCGATACGCCCGC  >  NC_008463/5386050‑5386348

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