BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,903,349	T→G	5.5%	R90R (CGA→CGC)	PP_2555 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,903,349	0	T	G	5.5%	388.3 / 9.7	135	R90R (CGA→CGC)	PP_2555	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (72/55); new base G (3/5); total (75/60)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

Only 100 of 136 total aligned reads displayed.

GTACGCATGGCCGACAACAGTTCAATGATGTCCTGGTGCAACAGCGCGAAGTTGGCGCAGATCAGGTCATAGTCCTTGCCAACAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCAGCAGCCAGCCCTCACCGCAGCCCAGGTCGAGCACCCGTTCGGGCTGGCGGCCAAGGATGGCCAGCAGGATGGCCTG  >  NC_002947/2903203‑2903498

gTACGCATGGCCGACAACAGTTCAATGATGTCCTGGTGCAACAGCGCGAAGTTGGCGCAGATCAGGTCATAGTCCTTGCCAACAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTcgc                                                                                                                                                     >  1:2380586/1‑150 (MQ=255)

   cGCATGGCCGACAACAGTTCAATGATGTCCTGGTGCAACAGCGCGAAGTTGGCGCAGATCAGGTCATAGTCCTTGCCAACAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCt                                                                                                                                                  <  2:1898118/150‑1 (MQ=255)

     cATGGCCGACAACAGTTCAATGATGTCCTGGTGCAACAGCGCGAAGTTGGCGCAGATCAGGTCATAGTCCTTGCCAACAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGc                                                                                                                                                <  1:2447350/150‑1 (MQ=255)

        ggCCGACAACAGTTCAATGATGTCCTGGTGCAACAGCGCGAAGTTGGCGCAGATCAGGTCATAGTCCTTGCCAACAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCGCGAGCTGcggg                                                                                                                                             <  1:2128349/150‑4 (MQ=255)

          ccGACAACAGTTCAATGATGTCCTGGTGCAACAGCGCGAAGTTGGCGCAGATCAGGTCATAGTCCTTGCCAACAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATcca                                                                                                                                           <  1:2351249/150‑1 (MQ=255)

               aacaGTTCAATGATGTCCTGGTGCAACAGCGCGAAGTTGGCGCAGATCAGGTCATAGTCCTTGCCAACAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGc                                                                                                                                      <  2:506632/150‑1 (MQ=255)

                    ttCAATGATGTCCTGGTGCAACAGCGCGAAGTTGGCGCAGATCAGGTCATAGTCCTTGCCAACAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCt                                                                                                                                 >  1:979734/1‑150 (MQ=255)

                    ttCAATGATGTCCTGGTGCAACAGCGCGAAGTTGGCGCAGATCAGGTCATAGTCCTTGCCAACAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCt                                                                                                                                 >  2:209606/1‑150 (MQ=255)

                       aaTGATGTCCTGGTGCAACAGCGCGAAGTTGGCGCAGATCAGGTCATAGTCCTTGCCAACAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTc                                                                                                                              >  1:559454/1‑150 (MQ=255)

                          gatgTCCTGGTGCAACAGCGCGAAGTTGGCGCAGATCAGGTCATAGTCCTTGCCAACAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGcc                                                                                                                           >  1:1432892/1‑150 (MQ=255)

                              tCCTGGTGCAACAGCGCGAAGTTGGCGCAGATCAGGTCATAGTCCTTGCCAACAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCg                                                                                                                       >  1:652971/1‑150 (MQ=255)

                               ccTGGTGCAACAGCGCGAAGTTGGCGCAGATCAGGTCATAGTCCTTGCCAACAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGa                                                                                                                      >  1:1564143/1‑150 (MQ=255)

                                cTGGTGCAACAGCGCGAAGTTGGCGCAGATCAGGTCATAGTCCTTGCCAACAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGAc                                                                                                                     >  1:1266560/1‑150 (MQ=255)

                                cTGGTGCAACAGCGCGAAGTTGGCGCAGATCAGGTCATAGTCCTTGCCAACAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGAc                                                                                                                     >  2:719899/1‑150 (MQ=255)

                                   gTGCAACAGCGCGAAGTTGGCGCAGATCAGGTCATAGTCCTTGCCAACAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGa                                                                                                                      >  1:2247839/1‑146 (MQ=255)

                                   gTGCAACAGCGCGAAGTTGGCGCAGATCAGGTCATAGTCCTTGCCAACAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGa                                                                                                                      <  2:2247839/146‑1 (MQ=255)

                                       aaCAGCGCGAAGTTGGCGCAGATCAGGTCATAGTCCTTGCCAACAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACg                                                                                                              >  2:211074/1‑150 (MQ=255)

                                        aCAGCGCGAAGTTGGCGCAGATCAGGTCATAGTCCTTGCCAACAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACgg                                                                                                             <  2:168604/150‑1 (MQ=255)

                                        aCAGCGCGAAGTTGGCGCAGATCAGGTCATAGTCAGTGCCAACAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACgg                                                                                                             <  1:12040/150‑1 (MQ=255)

                                                 aGTTGGCGCAGATCAGGTCATAGTCCTTGCCAACAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATAc                                                                                                    <  2:1361984/150‑1 (MQ=255)

                                                   ttGGCGCAGATCAGGTCATAGTCCTTGCCAACAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGcc                                                                                                                  >  1:2445535/1‑134 (MQ=255)

                                                   ttGGCGCAGATCAGGTCATAGTCCTTGCCAACAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGcc                                                                                                                  <  2:2445535/134‑1 (MQ=255)

                                                       cgcAGATCAGGTCATAGTCCTTGCCAACAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGt                                                                                              <  2:936722/150‑1 (MQ=255)

                                                        gcAGATCAGGTCATAGTCCTTGCCAACAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGATCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTc                                                                                             >  2:110745/1‑150 (MQ=255)

                                                           gATCAGGTCATAGTCCTTGCCAACAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCGCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGc                                                                                          <  2:559454/150‑1 (MQ=255)

                                                           aaTCAGGTCATAGTCCTTGCCAACAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTATGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGc                                                                                          <  2:1050989/149‑1 (MQ=255)

                                                             tCAGGTCATAGTCCTTGCCAACAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCa                                                                                        >  1:1142480/1‑150 (MQ=255)

                                                               aGGTCATAGTCCTTGCCAACAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAgc                                                                                      <  1:1667253/150‑1 (MQ=255)

                                                                ggTCATAGTCCTTGCCAACAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATa                                                                                                     >  2:360669/1‑134 (MQ=255)

                                                                 gTCATAGTCCTTGCCAACAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCg                                                                                            <  2:1315459/142‑1 (MQ=255)

                                                                 gTCATAGTCCTTGCCAACAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCg                                                                                            >  1:1315459/1‑142 (MQ=255)

                                                                 gTCATAGTCCTTGCCAACAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAgcgc                                                                                    <  2:1984278/150‑1 (MQ=255)

                                                                    aTAGTCCTTGCCAACAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAgcgc‑gcg                                                                                >  2:397239/1‑150 (MQ=255)

                                                                    aTAGTCCTTGCCAACAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGACGAGTCCCACGGCCTCGATACCCCGGTCGGCCAgcgc‑gcg                                                                                >  1:500282/1‑150 (MQ=255)

                                                                     tAGTCCTTGCCAACAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAgcgc‑gcgc                                                                               >  1:19396/1‑150 (MQ=255)

                                                                        tCCTTGCCAACAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCgcagc                                                                            >  1:1248708/1‑150 (MQ=255)

                                                                        tCCTTGCCAACAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCgcagc                                                                            >  2:668943/1‑150 (MQ=255)

                                                                          cTTGCCAACAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTc                                                                                                         >  1:2094862/1‑120 (MQ=255)

                                                                          cTTGCCAACAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGcagcag                                                                          >  2:1414214/1‑150 (MQ=255)

                                                                          cTTGCCAACAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGcagcag                                                                          >  1:730045/1‑150 (MQ=255)

                                                                          cTTGCCAACAGGCACCTTCGCCTCAACCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGcagcag                                                                          >  1:1519415/1‑150 (MQ=255)

                                                                           ttGCCAACAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGCTACCCCGGTCGGCCAGCGC‑GCGCagcagc                                                                         >  1:2426072/1‑150 (MQ=255)

                                                                           ttGCCAACAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCagcagc                                                                         >  2:1753911/1‑150 (MQ=255)

                                                                           ttGCCAACAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCagcagc                                                                         >  2:891313/1‑150 (MQ=255)

                                                                               cAACAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCAGcagccagc                                                                     >  1:389157/1‑150 (MQ=255)

                                                                               cAACAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCAGcagccagc                                                                     <  1:255255/150‑1 (MQ=255)

                                                                                aaCAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCAGCagccagcc                                                                    <  2:771601/150‑1 (MQ=255)

                                                                                   aGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCAGCAGCCAGCCCTc                                                                 >  2:1527245/1‑150 (MQ=255)

                                                                                     gCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCt                                                                                                          >  1:1461161/1‑108 (MQ=255)

                                                                                     gCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCAGCAGCCAGCCCTCAc                                                               >  1:356020/1‑150 (MQ=255)

                                                                                       aCCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCagcagc                                                                         >  1:726464/1‑138 (MQ=255)

                                                                                       aCCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCAGCCGCCAGACCTccccc                                                             >  2:374489/1‑149 (MQ=255)

                                                                                             gCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCAGCAGCCAGCCCTCACCGCa                                                           <  2:370866/146‑1 (MQ=255)

                                                                                             gCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCAGCAGCCAGCCCTCACCGCa                                                           >  1:370866/1‑146 (MQ=255)

                                                                                             gCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCAGCTGCGCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCAGCAGCCAGCCCTCACCGCAGccc                                                       >  1:436186/1‑150 (MQ=255)

                                                                                              ccTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCAGCAGCCAGCCCTCACCGCAGCCCa                                                      <  2:1768856/150‑1 (MQ=255)

                                                                                                    ccAACTGCGCATAGCTGGCCAGGTGCATCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTATTGCGCGGTCGGCGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCAGCAGCCAGCCCTCACCGCAGCCCAGGTCGa                                                <  2:566487/150‑1 (MQ=255)

                                                                                                    ccAACTGCGCATAGCTGGCCAGGTGCACCTCTCCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCAGCAGCCAGCCCTCACCGCAGCCCAGGTCGa                                                <  2:1248708/150‑1 (MQ=255)

                                                                                                    ccAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCAGCAGCCAGCCCTCACCGCAGCCCAGGTCGa                                                <  1:1414214/150‑1 (MQ=255)

                                                                                                     cAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGATCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCAGCAGCCAGCCCTCACCGCAGCCCAGGTCGAg                                               <  1:110745/150‑1 (MQ=255)

                                                                                                       aCTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCAGCAGCCAGCCCTCACCGCAGCCCAGGTCGAGCa                                             >  1:1719140/1‑150 (MQ=255)

                                                                                                        cTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCAGcagcca                                                                       <  1:2176640/123‑1 (MQ=255)

                                                                                                        cTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCAGcagcca                                                                       >  2:2176640/1‑123 (MQ=255)

                                                                                                        cTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCAGCAGCCAGCCCTCACCGCAGCCCAGGTCGAGCAc                                            <  2:19396/150‑1 (MQ=255)

                                                                                                         tGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATAccc                                                                                                  >  1:2276876/1‑96 (MQ=255)

                                                                                                         tGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCAGCAGCCAGCCCTCACCGCAGCCCAGGTCGAGCAcc                                           >  2:1993237/1‑150 (MQ=255)

                                                                                                         ggcgcATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATAccc                                                                                                  <  2:2276876/95‑1 (MQ=255)

                                                                                                              aTAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCAGCAGCCAGCCCTCACCGCAGCCCAGGTCGAGCACCCGTTc                                      <  1:668943/150‑1 (MQ=255)

                                                                                                               tAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCAGCAGCCAGCCCTCACCGCAGCCCAGGTCGAGCACCCGTTCg                                     <  1:30589/150‑1 (MQ=255)

                                                                                                               tAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTATGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCAGCAGCCAGCCCTCACCGCAGCCCAGGTCGAGCACCCGTTCg                                     <  2:866602/150‑1 (MQ=255)

                                                                                                                aGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCg                                                                                                >  1:578443/1‑91 (MQ=255)

                                                                                                                   tGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGc                                                                                          >  1:1962424/1‑94 (MQ=255)

                                                                                                                    ggCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCCCGATCCCCCGGTCGGCCAGCGCGGCGGTGAAG‑CACCTCGCACCGCAGCCCAGGTCGAGCACCCGTTCGGGCTg                                <  1:416134/150‑1 (MQ=255)

                                                                                                                      ccAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGtt                                                                                                                                   <  2:154055/50‑1 (MQ=255)

                                                                                                                      ccAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGtt                                                                                                                                   >  1:154055/1‑50 (MQ=255)

                                                                                                                       cAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCAGCAGCCAGCCCTCACCGCAGCCCAGGTCGAGCACCCGTTCGGGCTggcg                             <  1:395064/150‑1 (MQ=255)

                                                                                                                        aGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCAGCAGCCAGCCCTCACCGCAGCCCAGGTCGAGCACCCGTTCGGGCTgncgg                            <  2:1519415/150‑1 (MQ=255)

                                                                                                                          gTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCAGCAGCCAGCCCTCACCGCAGCCCAGGTCGAGCACCCGTTCGGGCTGGCGGcc                          >  2:405454/1‑150 (MQ=255)

                                                                                                                              aCCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCAGCAGCCAGCCCTCACCGCAGCCCAGGTCGAGCACCCGTTCGGGCt                                 >  1:2303289/1‑139 (MQ=255)

                                                                                                                              aCCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCAGCAGCCAGCCCTCACCGCAGCCCAGGTCGAGCACCCGTTCGGGCt                                 <  2:2303289/139‑1 (MQ=255)

                                                                                                                              aCCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCAGCAGCCAGCCCTCACCGCAGCCCAGGTCGAGCACCCGTTCGGGCTGGCGGCCAAgg                      >  1:1417715/1‑150 (MQ=255)

                                                                                                                               ccTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCAGCAGCCAGCCCTCACCGCAGCCCAGGTCGAGCACCCGTTCGGGCTGGCGGCCAAGGa                     >  1:828947/1‑150 (MQ=255)

                                                                                                                                cTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCAGCAGCCAGCCCTCACCGCAGCCCAGGTCGAGCACCCGTTCGGGCTGGCGGCCAAGGAt                    >  1:251449/1‑150 (MQ=255)

                                                                                                                                 tCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCAGCGGCCAGCCCTCACCGCAGCCCAGGTCGAGCACCCGTTCGGGCTGGCGGCCAAGGATg                   >  2:2394068/1‑150 (MQ=255)

                                                                                                                                 tCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCAGCAGCCAGCCCTCACCGCAGCCCAGGTCGAGCACCCGTTCGGGCTGGCGGCCAAGGATg                   >  1:1563675/1‑150 (MQ=255)

                                                                                                                                 tCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCAGCAGCCAGCCCTCACCGCAGCCCAGGTCGAGCACCCGTTCGGGCTGGCGGCCAAGGATg                   <  2:652971/150‑1 (MQ=255)

                                                                                                                                  cACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCAGCAGCCAGCCCTCACCGCAGCCCAGGTCGAGCACCCGTTCGGGCTGGCGGCCAAGGATgg                  <  2:2449557/150‑1 (MQ=255)

                                                                                                                                   aCCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGAc                                                                                                                     >  1:970459/1‑51 (MQ=255)

                                                                                                                                   aCCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCAGCAGCCAGCCCTCACCGCAGCCCAGGTCGAGCACCCGTTCGGGCTGGCGGCCAAGGATg                   >  1:1069384/1‑148 (MQ=255)

                                                                                                                                   aCCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCAGCAGCCAGCCCTCACCGCAGCCCAGGTCGAGCACCCGTTCGGGCTGGCGGCCAAGGATg                   <  2:1069384/148‑1 (MQ=255)

                                                                                                                                     ccGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCAGCAGCCAGCCCTCACCGCAGCCCAGGTCGAGCACCCGTTCGGGCTGGCGGCCAAGGATGGCca               >  1:1042012/1‑150 (MQ=255)

                                                                                                                                       gcgcCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAgcgc‑gc                                                                                 <  2:1104805/82‑1 (MQ=255)

                                                                                                                                          gatgCAGCGCGAGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCAGCAGCCAGCCCTCACCGCAGCCCAGGTCGAg                                               <  2:1678182/110‑1 (MQ=255)

                                                                                                                                          gatgCAGCGCGAGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCAGCAGCCAGCCCTCACCGCAGCCCAGGTCGAg                                               >  1:1678182/4‑113 (MQ=255)

                                                                                                                                            cGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCAGCAGCCAGCCCTCACCGCAGCCCAGGTCGAGCACCCGTTCGGGCTGGCGGCCAAGGATGGCCAGCAGGAt        >  2:859852/1‑150 (MQ=255)

                                                                                                                                             gCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCAGCGGCCAGCCCTCACCGCAGCCCAGGTCGAGCACCCGTTCGGGCTGGCGGCCAAGGATGGCCAGCAGGATg       <  1:2394068/150‑1 (MQ=255)

                                                                                                                                             gCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAACGC‑GCGCAGCAGCCAGCCCTCACCGCAGCCCAGGTCGAGCACCCGTTCGGGCTGGCGGCCAAGGATGGCCAGCAGGATg       >  1:26030/1‑150 (MQ=255)

                                                                                                                                             gCAGCGCGAGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCAGCAGCCAGCCCTCAc                                                               <  1:1062309/94‑1 (MQ=255)

                                                                                                                                               aGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCAGCAGCCAGCCCTCACCGCAGCCCAGGTCGAGCACCCGTTCGGGCTGGCGGCCAAGGATGGCCAGCAGGATGGc     <  2:1432892/150‑1 (MQ=255)

                                                                                                                                                 cTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCAGCAGCCAGCCCTCACCGCAGCCCAGGTCGAGCACCCGTTCGGGCTGGCGGCCAAGGATGGCCAGCAGGATGGCCt   >  1:596680/1‑150 (MQ=255)

GTACGCATGGCCGACAACAGTTCAATGATGTCCTGGTGCAACAGCGCGAAGTTGGCGCAGATCAGGTCATAGTCCTTGCCAACAGGCACCTTCGCCTCAGCCAACTGCGCATAGCTGGCCAGGTGCACCTCACCCGCGCCCGCAGCTCGCGCTGCATCCACCAGCGTTCTGTCGCCGTCGACGCCCACGGCCTCGATACCCCGGTCGGCCAGCGC‑GCGCAGCAGCCAGCCCTCACCGCAGCCCAGGTCGAGCACCCGTTCGGGCTGGCGGCCAAGGATGGCCAGCAGGATGGCCTG  >  NC_002947/2903203‑2903498

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