BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,163,456	C→G	7.7%	A158G (GCG→GGG)	PP_3661 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,163,456	0	C	G	7.7%	515.3 / 3.1	189	A158G (GCG→GGG)	PP_3661	membrane protein
Reads supporting (aligned to +/- strand): ref base C (92/80); new base G (10/7); total (102/87)
Fisher's exact test for biased strand distribution p-value = 8.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 202 total aligned reads displayed.

ATCGTGTCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACT  >  NC_002947/4163308‑4163604

aTCGTGTCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGgcg                                                                                                                                                     <  1:771407/150‑1 (MQ=255)

aTCGTGTCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGgcg                                                                                                                                                     >  1:583102/1‑150 (MQ=255)

 tCGTGTCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGgggc                                                                                                                                                    >  1:4378387/1‑150 (MQ=255)

 tCGTGTCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGgcgc                                                                                                                                                    <  1:1766498/150‑1 (MQ=255)

 tCGTGTCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGAATTCGAATCGgcgc                                                                                                                                                    <  2:5240670/150‑1 (MQ=255)

       cGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCTgcgcc                                                                                                                                              >  1:2635624/1‑149 (MQ=255)

       cGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCa                                                                                                                                              <  1:3891259/150‑1 (MQ=255)

         gTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAgc                                                                                                                                            >  2:2769853/1‑150 (MQ=255)

          tGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAgcg                                                                                                                                           >  2:2656451/1‑150 (MQ=255)

             tGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGccc                                                                                                                                        <  1:251850/150‑1 (MQ=255)

             tGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGccc                                                                                                                                        >  2:1314374/1‑150 (MQ=255)

              ggTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCGcccc                                                                                                                                       >  1:4516961/1‑150 (MQ=255)

              ggTAGTCACCAGTACCCTCATGATCGGCCTGTGGTGTGGCATGAAATTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGAGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGcccc                                                                                                                                       >  2:2989709/1‑150 (MQ=255)

               gttgTCACCTGTACCCTGCTGTTGGGCCTGGGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATCCCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCGGTGCTGGCATTCGATTCGGCGCTGCGCAGCGccccc                                                                                                                                      <  2:1584388/147‑1 (MQ=255)

                tAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGcccccc                                                                                                                                     >  2:4273992/1‑150 (MQ=255)

                 aGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCTTTCGAATCGGCGCTGCGCAGCGCCCCCCa                                                                                                                                    <  2:4143925/150‑1 (MQ=255)

                   tCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCGCGCCCCCCaaa                                                                                                                                  >  2:348193/1‑148 (MQ=255)

                   tCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGcccccccca                                                                                                                                  >  1:4733682/1‑147 (MQ=255)

                      cctgTACCCTGCTGATCGGCCTGTGGTGTGGCTTGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGGTGACCGCCGCCGGCTGCGCATTCTGCGGCGCCGCGGCCTTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATaaaa                                                                                                                               <  2:70394/147‑1 (MQ=255)

                      ccAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATaaaa                                                                                                                               <  2:2094784/150‑1 (MQ=255)

                       cAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGTCGCAGCGCCCCCCagcaaac                                                                                                                              >  1:4482228/1‑144 (MQ=255)

                        aGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGGGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAgc                                                                                                                             <  2:2001075/150‑1 (MQ=255)

                        aGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCGCCCCCCAAAAAAgc                                                                                                                             >  1:3710544/1‑150 (MQ=255)

                             ccTGCTGTTGGGCGTGTGTTGTGGCATTAATGTGTTCATGCTCGACGGCGATACCGCCTTGCTGACCGCGGCGGGCTGCGCCTTCTGCGGCGCCGCCGCCGTGCTGGCTTTCGAATGGGCGCTGCGATGCGCCCCCCATAAAAGCGCCAt                                                                                                                        <  2:5340179/150‑1 (MQ=255)

                             ccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCGCCCCCCAAAAAAGCGCCAt                                                                                                                        >  2:2514484/1‑150 (MQ=255)

                             ccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCCCCCCCCCAaaaaacaccca                                                                                                                        >  1:201095/1‑143 (MQ=255)

                              ctgcttATCGGCCTGTGGTGTGGCATTAATGTTTTCATGCTCGACCGCTTTACCGCGTTGCTGTCGGCCGCGGGCTGCGCCATCTGGGGCGCCGCCGCCGTGCTGGCATTCGAATGGGGGCTGCGCAGCGCCCCCCATAAAAGCGCCATg                                                                                                                       <  2:4564159/150‑1 (MQ=255)

                               ttgtgATCGGCCTGTGGTGTGGCATTAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGAAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGAATTCGAATGGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATgg                                                                                                                      <  2:5497350/147‑1 (MQ=255)

                                 cctaTCTGCCTGTGGTGTGGCATGAATGTGTTAAAGCTCGACCGCGTTACGGCGTTGCTGACCGCCGCGGGCTGCGCAATGTGGGGCGCGGCCGCGGTGGTGGCATTCTATTCGGGGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCa                                                                                                                    <  2:1136256/147‑1 (MQ=255)

                                   gATCGGCCTGTGGTGTGGCTTGAATGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCTTATGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGt                                                                                                                  <  2:4205296/150‑1 (MQ=255)

                                     tCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGCGCAGCGCCCCCCCCAAAAGCGCCCAGGCAGGCg                                                                                                                >  1:5196808/1‑150 (MQ=255)

                                        gCCTGTGGTGTGGCATGATAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGGCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGGTGGCATTCTAATCGGCGGTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGaa                                                                                                             <  2:4615610/150‑1 (MQ=255)

                                        gCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCAGCCCCGCGCCCCCCAAAAAAAACCCCAGGCAGTCGGaa                                                                                                             >  1:3189551/1‑150 (MQ=255)

                                        gCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCATTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCTTTCGATTCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGaa                                                                                                             <  2:4076258/150‑1 (MQ=255)

                                             tggtgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCgga                                                                                                        >  2:3545183/1‑148 (MQ=255)

                                              ggtgtggCATGAATGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCGGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAg                                                                                                       <  2:632627/150‑1 (MQ=255)

                                                tgtggAATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCGGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATGGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTg                                                                                                     <  2:166555/150‑1 (MQ=255)

                                                  tggCATGAATGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCTGCGCCATCTGCGGCGCCGCCGCCTTGTTGGCATTCGAATCGGCGCTGCGAAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCt                                                                                                   <  2:439563/150‑1 (MQ=255)

                                                   ggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCCCCCCCCAAAAAACCGCCAAGGCAGTCGGAAGCCTAGTGCTg                                                                                                  >  1:3361671/1‑150 (MQ=255)

                                                      aTGAACGTGTTCAAGCTCGACCGCGATACCGCCTGGCTGCCCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTc                                                                                               <  2:2850422/150‑1 (MQ=255)

                                                      aTGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAAAAAAGCGCCAAGGCAGTCGGAAGCGTAGTGCGGTTc                                                                                               >  2:2377334/1‑150 (MQ=255)

                                                       tGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCg                                                                                              <  2:4547777/150‑1 (MQ=255)

                                                        gAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAGAAAAGCGCCAAGGCACTCGGAAGCGTAGTGCTGTTAgg                                                                                             >  1:5499842/1‑150 (MQ=255)

                                                        gAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCGCCCCCCATAAAAGCGCCCAGGCAGTCGGACGCGTAGTGGTGTTCgg                                                                                             >  1:4559664/1‑150 (MQ=255)

                                                          aaGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGCGCACCGCCCCCCAAAAAAACGCCAAGGCACTCGGAACCGGAGTGCTGTTAGGTa                                                                                           >  1:4486714/1‑150 (MQ=255)

                                                              gtTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCCCCCCCCAAAAAAGCCCCCAGGCACTCGGAACCGGAGTGGTGGTCGGTAccgt                                                                                       >  1:2886428/1‑148 (MQ=255)

                                                               ttCAAGCTCGACTGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCCTAAAAGCGCCATGGCCGGCGGGAGCGTGGTGGTGGTCGGTACCtgg                                                                                      >  2:5485320/1‑148 (MQ=255)

                                                               ttCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCGCGCCCCCCCAAAAACCGCCCAGGCAGGCGGTAGCCTAGGGCGGGTCGGGGCCTTg                                                                                      >  1:5399519/1‑150 (MQ=255)

                                                                tCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCAAAAAAGCGCCATGGCAGTCGGAAGCGGAGTGCTGTTCGGTACCATGt                                                                                     >  2:3206678/1‑150 (MQ=255)

                                                                 cAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGGTCGGTACCTTGTc                                                                                    >  2:1883817/1‑150 (MQ=255)

                                                                 cAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCGCCCCCCCAAAAAGCGCCCAGGCCGTCGGAAGCCTAGTGCGGTTAGGTAccaagtg                                                                                    >  1:2467864/1‑145 (MQ=255)

                                                                 cAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCGCCCCCCCAAAAAGCGCCATGGCAGTCGGAAGCGGAGTGGTGGTCGGTACCgtgtg                                                                                    >  1:5014195/1‑149 (MQ=255)

                                                                  aaGCTCGACCGCTTTACCGCCTTGCTGCCCGCCGCGGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAAGCGGCGCTGCGCACCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCg                                                                                   <  1:1258546/150‑1 (MQ=255)

                                                                  aaGCACGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTAGGCAGCGTAGTGCGGTTCGGTACCTTGTCg                                                                                   >  2:2395696/1‑150 (MQ=255)

                                                                    gCTCGCCCGCTTTCCCGCCTTGCTGCCCGCGGCCGGCAGCGCCCTGGGGGGCGCCGCCGCCGGGGTGTCTTTCGACTCGGGGCGGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGAt                                                                                 <  1:781766/150‑1 (MQ=255)

                                                                      tCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCGCCCCCCCGAAAATCGCCCGGGCCGTCGGAAGCGTAGTGGTGTTCGGTACCGTGTCGgtgt                                                                               >  1:2800317/1‑150 (MQ=255)

                                                                       cGACCGCGATACCGCCTTGCTGACCGCCGCGGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAAGCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGtt                                                                              <  2:3133225/150‑1 (MQ=255)

                                                                           cgcgATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCAAGGCAGTCGGAAGCGTAGTGCAGTTCGGTACCTTGGCGATGTTccat                                                                          >  2:5176060/1‑148 (MQ=255)

                                                                           cgcgATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAGAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCtt                                                                          >  2:4610358/1‑150 (MQ=255)

                                                                            gcgATACCGCCTTGTTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCAAGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCttt                                                                         >  2:2223784/1‑150 (MQ=255)

                                                                             cgATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCGGCGGCGCCGCCGCCGTGCTGGCTTTCGACGGGGCGCTGCCCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTa                                                                        <  1:1287356/150‑1 (MQ=255)

                                                                             cgACACCGCCTGGGTGCCCGCCGCCGCCAGCGCCAGCTGCGGCGCCGCCGCGGTGGTGGATTTCGACGGGGCGCGGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTa                                                                        <  1:1539502/150‑1 (MQ=255)

                                                                                 cccGCTTTGCGGCCCGCCGCCGGCAGCGCCATGGGGGGCGCCGCCGCGGTGCTGGCATTCGATGCGGCGCTGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCa                                                                    <  1:315020/149‑1 (MQ=255)

                                                                                   cGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGGGGCGCAGCGCCCCCCAGAAAAGCGCCAAGGCAGTCGGAAGCGTAGTGGTGGTGGGTACCTTGTCGATGTTGCGTTAACCAca                                                                  >  1:360662/1‑149 (MQ=255)

                                                                                   cGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCACCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACt                                                                  >  1:30279/1‑150 (MQ=255)

                                                                                    gCCTTGCTGCCCGCCGCGGGCCGCGCCATGTGGGGCGCCGCCGCGGTGCTGGAATTCTAAGGGGCGGTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTg                                                                 <  1:2109979/150‑1 (MQ=255)

                                                                                    gCCTTGCTGACCGCCGCCGGCCGCGCCATCTGCGGCGCCGCCGCGGTGGTGGAATTCGAATGGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTg                                                                 <  2:1403218/150‑1 (MQ=255)

                                                                                     ccTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCAGAAAAGCGCCAAGGCAGTCGGAAGCGTAGTGGTGTTCGGTACCATGTCGAGGTTCCGTTAACCACAgg                                                                >  1:1365106/1‑150 (MQ=255)

                                                                                     ccTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCAGTCGAATCGGCGCTGCGCCGCGCCCCCCCTACAAGCGCCATGGCAGGCGGAAGCGGAGTGCGGGTCGGTCCCGGGGCGCTGTGCCGTTACCCACAgg                                                                >  2:3705439/1‑150 (MQ=255)

                                                                                     cattGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAGAAAAGCGCCAAGGCAGACGGTAGCGAAGTGCTGTTCGGTAACTAGTCGATGTTCAGTTAACCACAgg                                                                >  2:2918931/3‑150 (MQ=255)

                                                                                      gtggCTGCCCGCCGCGGGCAGCGCATTCGGGGGCGCCGCCGCCGTGCTGTAATTCCAATGGGCGCGGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGc                                                               <  1:1895931/147‑1 (MQ=255)

                                                                                         gCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGAt                                                            <  2:3139849/150‑1 (MQ=255)

                                                                                             cccgccgccgGCAGCGCCCTCGGCGGCGCCGCCGCCGTGCTGGCATTCGAAGCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAAt                                                        <  1:4205643/149‑1 (MQ=255)

                                                                                             aCCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAAt                                                        >  2:2113545/1‑150 (MQ=255)

                                                                                               cgccgccgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCCTAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACAGGCGATCAATGc                                                      >  2:2441790/1‑150 (MQ=255)

                                                                                                 ccgccgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGGCGTGCTGTTCGGTACCTTGTCGAGGGTCCTTTACCCACTGGCGCTCAATGCCg                                                    >  2:4786625/1‑150 (MQ=255)

                                                                                                 ccgccgGCAGCGCCATCGGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCg                                                    <  2:3861487/150‑1 (MQ=255)

                                                                                                        cAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAAAAAAGCGCCAAGGCAGTCGGAAGCGTAGTGGTGTTCGGTACCTTGTCGATGTTACTTTACCCACAGGCGATAAAATCCggctggct                                             >  1:2783085/1‑150 (MQ=255)

                                                                                                        cAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCggctggct                                             <  2:4496017/150‑1 (MQ=255)

                                                                                                         aGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACAGGCGATCAATGCCGgctggctg                                            >  2:3972645/1‑150 (MQ=255)

                                                                                                         aGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCGGTTCGGTACCGTGTCGAAGTTCCTTTACCCACAGGCGATCAATGCCGgctggctg                                            >  2:2297956/1‑150 (MQ=255)

                                                                                                              cttCTGCGGCGCCGCCGCCGTGCTGGCATTCTAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCt                                       <  2:1158475/148‑1 (MQ=255)

                                                                                                                 cTGCGGCGCCGCCGCCGTGGTGGATTTCGAGTCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGa                                    <  1:3677117/150‑1 (MQ=255)

                                                                                                                   gcggcgCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTGGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATa                                  >  2:1767234/1‑150 (MQ=255)

                                                                                                                    cggcgCCGCCGCGGTGGGGGCATGCCACGCGGCGCGGCGCAGCCCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATAc                                 <  1:3856750/150‑1 (MQ=255)

                                                                                                                     ggcgCCGCCGCCGTGCTGGCATTGGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATAcc                                <  2:1616227/150‑1 (MQ=255)

                                                                                                                     ggcgCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAGAAAAGCGCCCTGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCGTTACCCACGGGCGATCAAAGCCGGCTGGCTGCATCTGGATAcc                                >  2:4967241/1‑150 (MQ=255)

                                                                                                                                cgTGCTGGCATTCGAAGCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCgg                     <  2:2601258/150‑1 (MQ=255)

                                                                                                                                 gTGCTGGCATTCTAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCggg                    <  2:1861174/150‑1 (MQ=255)

                                                                                                                                   gCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACGGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgct                  <  2:1484716/150‑1 (MQ=255)

                                                                                                                                          ttcgccgcggcgcggcgCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGctgctgctg           <  1:348193/143‑1 (MQ=255)

                                                                                                                                          ttCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGctgctgctg           >  2:2262632/1‑150 (MQ=255)

                                                                                                                                           tCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAACCGTAGTGCTGTTCGGTACCATGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCggg          >  2:3218374/1‑150 (MQ=255)

                                                                                                                                             gattgggcgcggcgcTGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggc        <  2:702102/145‑1 (MQ=255)

                                                                                                                                             gAATCGGCGCTGCGCAGCGCCCCCCATAAAATCGCCAAGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTAACCACTGGCGATCAAAGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggc        >  2:4038809/1‑150 (MQ=255)

                                                                                                                                             gAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGGTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggc        >  2:4646209/1‑150 (MQ=255)

                                                                                                                                             gAATCGGCGCTGCGCAGCGCCCCCCAAAAAATCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATTTTCAATTAACCACAGGCGATCAATGCCCGCTGGCTGCATCTGGATAACAATGGCGCCGGGCTGCTGCTGggc        >  1:827577/1‑150 (MQ=255)

                                                                                                                                              aaTCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTAGGTACCTTGTCGATGTTCCTTTACCCACAGGCGATCAATTCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggcg       >  1:382451/1‑150 (MQ=255)

                                                                                                                                                  ggCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCAc   >  1:4620335/1‑150 (MQ=255)

                                                                                                                                                   gcgcTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACt  >  2:1509614/1‑150 (MQ=255)

ATCGTGTCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACT  >  NC_002947/4163308‑4163604

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