BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,700,782	G→C	5.7%	P625R (CCC→CGC) ‡	kdpD ←	sensor protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,700,782	0	G	C	5.7%	639.3 / 19.6	191	P625R (CCC→CGC) ‡	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 196 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)

    gcCACGCATGGCGGTCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGcc                                                                                                                                              >  1:2531145/1‑150 (MQ=255)

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

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

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

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

             ggCGGTCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCCGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGt                                                                                                                                     <  1:3562301/150‑1 (MQ=255)

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

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

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

                              gcAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGc                                                                                                                                               >  1:3230379/1‑123 (MQ=255)

                              gcAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTg                                                                                                                    >  2:3957126/1‑150 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                         gATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGc                                                                                                         <  1:1347076/150‑1 (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)

                                               gAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGc                                                                                                   >  1:3733037/1‑150 (MQ=255)

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

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

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

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

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

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

                                                     agcagGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCg                                                                                                               >  1:4164353/1‑132 (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                                                                                             >  2:549575/1‑150 (MQ=255)

                                                     agcagGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTg                                                                                             >  1:5428758/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)

                                                          ggCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCgg                                                                                        <  1:3620062/150‑1 (MQ=255)

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

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

                                                              cTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGgcgc                                                                                    <  2:2922117/150‑1 (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)

                                                                     gTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCca                                                                             <  2:3475085/150‑1 (MQ=255)

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

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

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

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

                                                                          tCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCagcagc                                                                        >  1:1537240/1‑150 (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)

                                                                                  gCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCgg                                                                >  1:4075956/1‑150 (MQ=255)

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

                                                                                    cGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCt                                                              >  2:2589186/1‑150 (MQ=255)

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

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

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

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

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

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

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

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

                                                                                                        aGCTGTTCAGCCAGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGc                                          <  1:3867692/150‑1 (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)

                                                                                                                    aGGCGTGCGCGGGCCAGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAg                              >  2:3806329/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                >  1:384675/1‑150 (MQ=255)

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

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

                                                                                                                                   aGGG‑CCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAGGCAgg               <  2:3673226/150‑1 (MQ=255)

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

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

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

                                                                                                                                            ggCCAGGGGCTGGCCGAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAGGCAGGGTATCGCt       <  1:3286805/150‑1 (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: `‑`