BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,317,256	C→G	7.0%	P608P (CCG→CCC)	treSA ←	trehalose synthase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,317,256	0	C	G	7.0%	310.2 / 12.8	100	P608P (CCG→CCC)	treSA	trehalose synthase A
Reads supporting (aligned to +/- strand): ref base C (50/43); new base G (3/4); total (53/47)
Fisher's exact test for biased strand distribution p-value = 7.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.44e-01

Only 100 of 107 total aligned reads displayed.

GACCGTCTCGTTGATCATGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCGTATAGCGCCCGGGCACGCG  >  NC_002947/3317109‑3317404

gACCGTCTCGTTGATCATGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCgg                                                                                                                                                    <  1:2394930/150‑1 (MQ=255)

gACCGTCTCGTTGATCATGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCgg                                                                                                                                                    <  2:371073/150‑1 (MQ=255)

     tctcGTTGATCATGTCTACCATCGGCCCAGGTGTGAACCCGGTCAACAGCAGATCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATGTGGATACCGCGCCCGGTTGGCAACTCATGCACCATAACCAGCCACCCCGGgctgc                                                                                                                                               >  2:556654/1‑150 (MQ=255)

     tctcGTTGATCATGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCGGGgctgc                                                                                                                                               <  1:1787335/150‑1 (MQ=255)

     tctcGTTGATCATGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGgctgc                                                                                                                                               >  1:618822/1‑150 (MQ=255)

      ctcGTTGATCATGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGgctgct                                                                                                                                              <  2:1407425/150‑1 (MQ=255)

      ctcGTTGATCATGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGgctgct                                                                                                                                              <  2:764869/150‑1 (MQ=255)

       tcGTTGATCATGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCGGGgctgctg                                                                                                                                             <  2:1580694/150‑2 (MQ=255)

                aTGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGt                                                                                                                                    >  1:814996/1‑150 (MQ=255)

                 tGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTa                                                                                                                                   <  2:539794/150‑1 (MQ=255)

                 tGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTa                                                                                                                                   >  1:2031327/1‑150 (MQ=255)

                        caccGGCCCAGGTGTGAACCCGGTTAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTAGCAGCGCATGTACCATCACCAGCAGCCCCGAGCTGCTCACCTCAGGTACAAGCAc                                                                                                                            >  2:477504/1‑150 (MQ=255)

                        caccGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGc                                                                                                                              <  1:1358178/148‑1 (MQ=255)

                        caccGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGc                                                                                                                              >  2:1358178/1‑148 (MQ=255)

                         accGGCCCAGGTGTGAACCCGGTCAACAGCAGTAACTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCAcc                                                                                                                           >  1:190176/1‑150 (MQ=255)

                          ccGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCt                                                                                                                          <  1:1228090/150‑1 (MQ=255)

                               ccAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCGGGgctgctgg                                                                                                                                            >  2:1943062/1‑125 (MQ=255)

                               ccAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCGGGgctgctgg                                                                                                                                            <  1:1943062/127‑3 (MQ=255)

                               ccAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACgg                                                                                                                     >  1:1432188/1‑150 (MQ=255)

                               ccAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACgg                                                                                                                     >  1:2290117/1‑150 (MQ=255)

                                cAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGc                                                                                                                    <  2:1528128/150‑1 (MQ=255)

                                cAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGc                                                                                                                    <  1:660564/150‑1 (MQ=255)

                                   gtgtGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCa                                                                                                                                       <  2:1716979/128‑1 (MQ=255)

                                   gtgtGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCa                                                                                                                                       >  1:1716979/1‑128 (MQ=255)

                                   gtgtGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCtgg                                                                                                                 >  2:2113524/1‑150 (MQ=255)

                                       gAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGc                                                                                                             >  1:389345/1‑150 (MQ=255)

                                       gAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGc                                                                                                             >  1:389360/1‑150 (MQ=255)

                                        aaCCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCg                                                                                                            >  1:184968/1‑150 (MQ=255)

                                        aaCCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCg                                                                                                            >  2:1782120/1‑150 (MQ=255)

                                            cGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGc                                                                                                        >  2:274840/1‑150 (MQ=255)

                                             ggTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGcc                                                                                                       <  2:618822/150‑1 (MQ=255)

                                              gTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGaggtg                                                                                                                    >  1:2060365/1‑132 (MQ=255)

                                              gTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGaggtg                                                                                                                    <  2:2060365/136‑5 (MQ=255)

                                                 aaCAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAg                                                                                                   <  2:389345/150‑1 (MQ=255)

                                                 aaCAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAg                                                                                                   <  2:389360/150‑1 (MQ=255)

                                                  aCAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAgg                                                                                                  >  1:1952604/1‑150 (MQ=255)

                                                    agcagTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGcc                                                                                                >  2:867217/1‑150 (MQ=255)

                                                    agcagTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATACACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGcc                                                                                                >  2:2284838/1‑150 (MQ=255)

                                                        gTTCCTCGGCAATCGCGTCCTGCACGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGgc                                                                                            <  1:1819697/150‑1 (MQ=255)

                                                           ccTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGca                                                                                         >  1:1179680/1‑150 (MQ=255)

                                                             tCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGcaca                                                                                       >  2:1988720/1‑150 (MQ=255)

                                                             tCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGcaca                                                                                       >  2:2374824/1‑150 (MQ=255)

                                                             tCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGcaca                                                                                       >  2:1973529/1‑150 (MQ=255)

                                                             tCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGcaca                                                                                       >  1:2078268/1‑150 (MQ=255)

                                                               ggCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCAcagc                                                                                     <  2:184968/150‑1 (MQ=255)

                                                               ggCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCAcagc                                                                                     <  2:2407316/150‑1 (MQ=255)

                                                                          ccTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGa                                                                                                           <  1:691537/117‑1 (MQ=255)

                                                                          ccTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGa                                                                                                           >  2:691537/1‑117 (MQ=255)

                                                                           cTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTttctt                                                                         <  1:846711/150‑1 (MQ=255)

                                                                            tGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTttcttc                                                                        >  2:147448/1‑150 (MQ=255)

                                                                            tGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTttcttc                                                                        >  2:139654/1‑150 (MQ=255)

                                                                            tGGCCGAAGTTCAGCCCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTttcttc                                                                        <  2:1179680/150‑1 (MQ=255)

                                                                              gCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCAc                                                                      <  1:2374824/150‑1 (MQ=255)

                                                                                cGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCt                                                                    >  2:953262/1‑150 (MQ=255)

                                                                                 gAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCtt                                                                   <  1:867217/150‑1 (MQ=255)

                                                                                    gTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATACACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCa                                                                <  1:2284838/150‑1 (MQ=255)

                                                                                      tCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAgg                                                              <  1:1832424/150‑1 (MQ=255)

                                                                                           gcAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGAGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAAt                                                         <  2:814996/150‑1 (MQ=255)

                                                                                               tGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAAt                                                     >  2:1609688/1‑150 (MQ=255)

                                                                                               tGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAAt                                                     <  1:1988720/150‑1 (MQ=255)

                                                                                                 aTCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCa                                                   >  1:1760393/1‑150 (MQ=255)

                                                                                                 aTCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCa                                                   <  2:2419747/150‑1 (MQ=255)

                                                                                                 aTCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCa                                                   >  2:2097260/1‑150 (MQ=255)

                                                                                                  tCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCAc                                                  <  1:1782120/150‑1 (MQ=255)

                                                                                                   cTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCAct                                                 <  2:1432188/150‑1 (MQ=255)

                                                                                                   cTGGATACCGCGCCCCGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCAct                                                 <  1:1074863/150‑1 (MQ=255)

                                                                                                    tGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAgg                                                                                                  >  2:726278/1‑100 (MQ=255)

                                                                                                    tGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAgg                                                                                                  <  1:726278/100‑1 (MQ=255)

                                                                                                      gATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACtctc                                              <  1:1465696/150‑1 (MQ=255)

                                                                                                      gATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATCACATAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACtctc                                              >  2:917978/1‑150 (MQ=255)

                                                                                                             cgcCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGc                                       <  1:274840/150‑1 (MQ=255)

                                                                                                                    cTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGt                                <  2:2031327/150‑1 (MQ=255)

                                                                                                                        cAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCGGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAg                            <  2:47025/150‑1 (MQ=255)

                                                                                                                          gCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGc                          >  2:42339/1‑150 (MQ=255)

                                                                                                                          gCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGc                          >  1:1462713/1‑150 (MQ=255)

                                                                                                                            tCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTg                        >  1:497435/1‑150 (MQ=255)

                                                                                                                                gCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCGt                    >  2:2136064/1‑150 (MQ=255)

                                                                                                                                 cACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTttc                                                                           >  1:509440/1‑94 (MQ=255)

                                                                                                                                  aCCNTCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCGtat                  >  2:1508193/1‑150 (MQ=255)

                                                                                                                                  aCCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCGtat                  >  1:421747/1‑150 (MQ=255)

                                                                                                                                    cATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCGTATAg                >  2:1391138/1‑150 (MQ=255)

                                                                                                                                       cACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAg                                                                                                   <  1:1073/64‑1 (MQ=255)

                                                                                                                                       cACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAg                                                                                                   >  2:1073/1‑64 (MQ=255)

                                                                                                                                        aCCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCGTATAGCGcc            >  2:1272725/1‑150 (MQ=255)

                                                                                                                                         ccAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCGTATAGCGccc           <  2:348775/150‑1 (MQ=255)

                                                                                                                                         ccAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGCCCAGGCTGCCGTATAGCGccc           >  2:775943/1‑150 (MQ=255)

                                                                                                                                          gagcagcCCGGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCGtata                 <  2:452963/142‑1 (MQ=255)

                                                                                                                                          gagcagcCCGGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCGtata                 >  1:452963/2‑143 (MQ=255)

                                                                                                                                          cagcagCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCGTATAGCGCCCg          <  2:497435/150‑1 (MQ=255)

                                                                                                                                          cagcagCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCGTATAGCGCCCg          >  2:622308/1‑150 (MQ=255)

                                                                                                                                           agcagcCCGGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCGTATAGCGCCCgg         >  2:1289813/1‑150 (MQ=255)

                                                                                                                                           agcagcCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCGTATAGCGCCCgg         <  1:1559019/150‑1 (MQ=255)

                                                                                                                                            gcagcCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCGTATAGCGCCCggg        <  1:1272725/150‑1 (MQ=255)

                                                                                                                                                ccccGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCGTATAGCGCCCGGGCAcg    >  2:212868/1‑150 (MQ=255)

                                                                                                                                                 cccGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAAt                                                     >  2:714892/1‑100 (MQ=255)

                                                                                                                                                 cccGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAAt                                                     <  1:714892/100‑1 (MQ=255)

                                                                                                                                                 cccGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCGTATAGCGCCCGGGCAcgc   <  1:477504/150‑1 (MQ=255)

                                                                                                                                                  ccGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTttc                                                                           >  1:534294/1‑77 (MQ=255)

                                                                                                                                                  ccGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTttc                                                                           <  2:534294/77‑1 (MQ=255)

                                                                                                                                                  acGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTAGTCCAGCTGCCTGTCCAGGCTGCCGTATAGCGCCCGGGCAcgcg  <  1:1762266/149‑1 (MQ=255)

GACCGTCTCGTTGATCATGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCGTATAGCGCCCGGGCACGCG  >  NC_002947/3317109‑3317404

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