BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,700,783	G→C	5.6%	P625A (CCC→GCC) ‡	kdpD ←	sensor protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,700,783	0	G	C	5.6%	637.1 / 18.9	191	P625A (CCC→GCC) ‡	kdpD	sensor protein
Reads supporting (aligned to +/- strand): ref base G (98/82); new base C (7/4); total (105/86)
Fisher's exact test for biased strand distribution p-value = 7.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 194 total aligned reads displayed.

TGCTGCCACGCATGGCGGTCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAGGCAGGGTATCGCTGCCGT  >  NC_002947/4700636‑4700929

tgctgcCACGCATGGCGGTCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCt                                                                                                                                                  <  1:1702160/150‑1 (MQ=255)

   tgcCACGCATGGCGGTCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGc                                                                                                                                               <  1:2684860/150‑1 (MQ=255)

       aCGCATGGCGGTCAACGGCGTGCGCAAGTCATGGGAGACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAg                                                                                                                                           <  1:4582208/150‑1 (MQ=255)

           aTGGCGGTCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTg                                                                                                                                                 <  1:2519657/140‑1 (MQ=255)

             ggCGGTCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGt                                                                                                                                     <  1:182534/150‑1 (MQ=255)

              gCGGTCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCCTAGGCCAGCCCCTGGCCCAggccctg                                                                                                                                     <  2:4642636/150‑6 (MQ=255)

               cGGTCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGa                                                                                                                                   <  2:4048683/150‑1 (MQ=255)

                   cAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAg                                                                                                                               <  2:4835281/150‑1 (MQ=255)

                   cAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAg                                                                                                                               <  1:3220209/150‑1 (MQ=255)

                        gCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGcc                                                                                                                          <  2:3515613/150‑1 (MQ=255)

                             cgcAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGt                                                                                                                     <  2:3471604/150‑1 (MQ=255)

                               cAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTgg                                                                                                                   <  1:2848745/150‑1 (MQ=255)

                                   tCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCg                                                                                                               >  2:1904318/1‑150 (MQ=255)

                                      tGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTg                                                                                                            >  1:174113/1‑150 (MQ=255)

                                      tGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTg                                                                                                            >  2:1806398/1‑150 (MQ=255)

                                        ggATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGt                                                                                                                                     >  2:1878667/1‑123 (MQ=255)

                                        ggATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAg                                                                                                          >  2:1716118/1‑150 (MQ=255)

                                        ggATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAg                                                                                                          >  2:1562824/1‑150 (MQ=255)

                                        ggATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAg                                                                                                          <  1:1527299/150‑1 (MQ=255)

                                         gATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGc                                                                                                         <  1:1347076/150‑1 (MQ=255)

                                         gATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGc                                                                                                         >  1:907889/1‑150 (MQ=255)

                                           tACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCTGGCGCCGGCGTGGGTCGCTGAGCgg                                                                                                       >  1:3599194/1‑150 (MQ=255)

                                             cGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTCGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCt                                                                                                     >  2:3157156/1‑150 (MQ=255)

                                              ggAAGCCAGCAGGGCGCTGCGTAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTg                                                                                                    >  2:2191866/1‑150 (MQ=255)

                                              ggAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTg                                                                                                    >  1:2127450/1‑150 (MQ=255)

                                              ggAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTg                                                                                                    >  1:1900132/1‑150 (MQ=255)

                                                aaGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGcc                                                                                                  >  2:1211910/1‑150 (MQ=255)

                                                   ccAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGt                                                                                               >  2:610633/1‑150 (MQ=255)

                                                    cagcagGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGGCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGGCCGtt                                                                                              <  1:1165915/150‑1 (MQ=255)

                                                    cagcagGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGtt                                                                                              <  2:1487324/150‑1 (MQ=255)

                                                    cagcagGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGtt                                                                                              >  1:2741513/1‑150 (MQ=255)

                                                     agcagGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCg                                                                                                               >  1:4164353/1‑132 (MQ=255)

                                                     agcagGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCg                                                                                                               <  2:4164353/132‑1 (MQ=255)

                                                     agcagGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGTGCGGCTGCCCGTTg                                                                                             >  2:2674585/1‑150 (MQ=255)

                                                     agcagGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTg                                                                                             <  1:119597/150‑1 (MQ=255)

                                                     agcagGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTg                                                                                             >  2:3398602/1‑150 (MQ=255)

                                                     agcagGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTg                                                                                             >  1:5428758/1‑150 (MQ=255)

                                                     agcagGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTg                                                                                             >  2:549575/1‑150 (MQ=255)

                                                      gcagGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTGGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCGGCCCGTTgc                                                                                            >  1:4110973/1‑150 (MQ=255)

                                                      gcagGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTgc                                                                                            >  2:4498758/1‑150 (MQ=255)

                                                       cagGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTgcg                                                                                           >  2:51446/1‑150 (MQ=255)

                                                        agGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGGTGTTCAGCCAGGCGTGCGAGGGCCAGGG‑CCTGGGCCAGGGGCTGGGCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTTAGGGGCTGCCCGTTgcgc                                                                                          <  1:2972453/150‑1 (MQ=255)

                                                            cgcTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTGGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGgc                                                                                      <  2:2958515/150‑1 (MQ=255)

                                                              cTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGgcgc                                                                                    <  2:2922117/150‑1 (MQ=255)

                                                              cTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGgcgc                                                                                    <  1:5418179/150‑1 (MQ=255)

                                                               tGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGca                                                                                   >  2:2370162/1‑150 (MQ=255)

                                                               tGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGca                                                                                   >  2:1245590/1‑150 (MQ=255)

                                                                gcgcAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGcac                                                                                  >  1:231251/1‑150 (MQ=255)

                                                                 cgcAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGcaca                                                                                 >  1:5094706/1‑150 (MQ=255)

                                                                   cAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACAcc                                                                               <  1:1977253/150‑1 (MQ=255)

                                                                       tGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCcagc                                                                           <  1:2961771/150‑1 (MQ=255)

                                                                        gCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCcagca                                                                          >  1:4327749/1‑150 (MQ=255)

                                                                        gCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCcagca                                                                          <  2:3692254/150‑1 (MQ=255)

                                                                         cTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCcagc                                                                           >  2:4329478/1‑148 (MQ=255)

                                                                          tCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCagcagc                                                                        >  2:3756133/1‑150 (MQ=255)

                                                                          tCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCagcagc                                                                        >  1:1537240/1‑150 (MQ=255)

                                                                            ggTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAgcgc                                                                      <  1:3266074/150‑1 (MQ=255)

                                                                             gTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGcc                                                                     <  2:5314048/150‑1 (MQ=255)

                                                                               ttCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAg                                                                   >  1:582710/1‑150 (MQ=255)

                                                                                tCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGc                                                                  >  1:4763866/1‑150 (MQ=255)

                                                                                 cGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCg                                                                 <  2:4073443/150‑1 (MQ=255)

                                                                                   ccGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGc                                                               <  1:5214699/150‑1 (MQ=255)

                                                                                     gTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTg                                                             <  2:582398/150‑1 (MQ=255)

                                                                                        cAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATc                                                          >  1:3740106/1‑150 (MQ=255)

                                                                                          ggCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGt                                                        <  2:1302619/150‑1 (MQ=255)

                                                                                           gCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCt                                                                                                             <  1:682669/96‑1 (MQ=255)

                                                                                           gCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCt                                                                                                             >  2:682669/1‑96 (MQ=255)

                                                                                              tGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTAc                                                    <  1:658062/150‑1 (MQ=255)

                                                                                                  gCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGcaca                                                                                 >  2:3868778/1‑117 (MQ=255)

                                                                                                  gCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGcaca                                                                                 <  1:3868778/117‑1 (MQ=255)

                                                                                                     tCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAAt                                             <  1:3862576/150‑1 (MQ=255)

                                                                                                      ccAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATg                                            <  2:4739275/150‑1 (MQ=255)

                                                                                                        aGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGc                                                                                                                                 <  1:3433984/63‑1 (MQ=255)

                                                                                                        aGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGc                                                                                                                                 >  2:3433984/1‑63 (MQ=255)

                                                                                                        aGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCaa                                              <  2:2530088/146‑1 (MQ=255)

                                                                                                        aGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCaa                                              >  1:2530088/1‑146 (MQ=255)

                                                                                                        aGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGc                                          >  2:2358906/1‑150 (MQ=255)

                                                                                                         gCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGcc                                         >  1:4468276/1‑150 (MQ=255)

                                                                                                             ttCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGccacca                                     >  2:1033310/1‑150 (MQ=255)

                                                                                                             ttCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGccacca                                     >  1:5139983/1‑150 (MQ=255)

                                                                                                              tCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGccaccac                                    <  2:332769/150‑1 (MQ=255)

                                                                                                                  ccAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGccaccaccacc                                >  2:852541/1‑150 (MQ=255)

                                                                                                                  ccAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGccaccaccacc                                >  1:1544681/1‑150 (MQ=255)

                                                                                                                     ggCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGc                                                 <  1:1020017/130‑1 (MQ=255)

                                                                                                                     ggCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGc                             >  1:750299/1‑150 (MQ=255)

                                                                                                                             taggCCAGCC‑CCTGGCCCAGGCCCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTa                     >  1:4525891/3‑150 (MQ=255)

                                                                                                                               ggCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAgg                   <  2:5243488/150‑1 (MQ=255)

                                                                                                                               ggCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAgg                   >  1:4341483/1‑150 (MQ=255)

                                                                                                                                 tcaccgcccTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAGGc                  >  1:4351675/8‑150 (MQ=255)

                                                                                                                                 ccAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAGGCa                 >  1:2683836/1‑150 (MQ=255)

                                                                                                                                 ccAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCAGTAGGCa                 >  1:3215132/1‑150 (MQ=255)

                                                                                                                                  cAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAGGCAg                <  2:456732/150‑1 (MQ=255)

                                                                                                                                  cAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAGGCAg                <  2:695013/150‑1 (MQ=255)

                                                                                                                                   cagccccTGGCCCAGGCCCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAGGCAg                <  2:4525891/145‑1 (MQ=255)

                                                                                                                                        ccTAGGCCAGCCCCTGGCCCAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAGGCAGGGTAt           >  1:3691656/1‑150 (MQ=255)

                                                                                                                                        ccTAGGCCAGCCCCTGGCCCAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAGGCAGGGTAt           >  1:2363020/1‑150 (MQ=255)

                                                                                                                                            ggCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAGGCAGGGTATCGCt       <  1:5275704/150‑1 (MQ=255)

                                                                                                                                            ggCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAGGCAGGGTATCGCt       >  1:1853571/1‑150 (MQ=255)

                                                                                                                                               cAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCaa                                              <  1:825017/108‑1 (MQ=255)

                                                                                                                                                 ggggCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAGGCAGGGTATCGCTGCCGt  >  2:2004717/1‑150 (MQ=255)

TGCTGCCACGCATGGCGGTCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAGGCAGGGTATCGCTGCCGT  >  NC_002947/4700636‑4700929

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