BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,298,647	G→A	5.2%	D506D (GAC→GAT)	sbcC ←	exonuclease subunit SbcC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,298,647	0	G	A	5.2%	349.5 / 6.2	115	D506D (GAC→GAT)	sbcC	exonuclease subunit SbcC
Reads supporting (aligned to +/- strand): ref base G (55/54); new base A (3/3); total (58/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.33e-01

Only 100 of 117 total aligned reads displayed.

GTTGCTGGCGTTGCAGCATGGCCTGCTGCTGCCGTTCGCGCAGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAGCTCTTCTCGGCCCTTGGTCAGCCGGCC  >  NC_002947/2298499‑2298796

gTTGCTGGCGTTGCCGCCTGGCCTGCTGCTTCGGTTTGCGCTGAGCGTCGGGCTGCTGCCGCAGCTCGTGTTCACGGCGGTGCAGGCGTGGCATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCCGGt                                                                                                                                                      <  2:1121124/150‑1 (MQ=255)

    cTGGCGTTGCAGCATGGCCTGCTGCTGCCGTTCGCGCAGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGCCATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGTCGAGCAGGTCGAt                                                                                                                                                  >  2:342792/1‑150 (MQ=255)

     tGGCGTTGCAGCATGGCCTGCTGCTGCCGTTCGCGCAGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGATCGacc                                                                                                                                                 >  1:615137/1‑148 (MQ=255)

      ggCGTTGCAGCATGGCCTGCTGCTGCCGTTCGCGCAGGGCATCCAGCTTCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCt                                                                                                                                                >  2:701436/1‑150 (MQ=255)

       gCGTTGCAGCATGGCCTGCTTCTGCCGTTCGCGCAGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTg                                                                                                                                               >  1:80775/1‑150 (MQ=255)

         gTTGCAGCATGGCCTGCTGCTGCCGTTCGCGCAGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCt                                                                                                                                             <  2:1628452/150‑1 (MQ=255)

          ttGCAGCATGGCCTGCTGCTGCCGTTCGCGCAGGGCATCCAGCTACTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGc                                                                                                                                              >  1:2315058/1‑148 (MQ=255)

           tGCAGCATGGCCTGCTGCTGCCGTTCGCGCAGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGATCGACCTGCTCg                                                                                                                                           <  1:1141439/150‑1 (MQ=255)

             cagcaTGGCCTGCTGCTGCCGTTCGCGCAGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGATCGACCTGCTCGGc                                                                                                                                         <  2:615137/150‑1 (MQ=255)

                  tGGCCTGCTGCTGCCGTTCGCGCAGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAggg                                                                                                                                    >  1:1540812/1‑150 (MQ=255)

                    gCCTGCTGCTGCCGTTCGCGCAGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGcc                                                                                                                                  >  2:1768706/1‑150 (MQ=255)

                     ccTGCTGCTGCCGTTCGCGCAGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCt                                                                                                                                 >  2:1265010/1‑150 (MQ=255)

                     ccTGCTGCTGCCGTTCGCGCAGGGCATCCAGATGCTGGCGACGCTCCTATACTCGGCGGTCCGGGCGGGCGCTACCCGCTTCGGCTCTCTGCTGGTTGCGGGTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGccc                                                                                                                                 >  2:2102533/1‑149 (MQ=255)

                      ctgctgctgCCGTTCGCGCAGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTg                                                                                                                                <  1:1068460/150‑1 (MQ=255)

                        gctgctgcCGTTCGCGCAGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGATCGATCTGCTCGGCCAGGGCCTggg                                                                                                                              >  1:798428/1‑150 (MQ=255)

                        gctgctgcCGTTCGCGCAGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCATGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAAGGccgaac                                                                                                                              >  1:459053/1‑146 (MQ=255)

                            ctgcCGTTCGCGCAGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTc                                                                                                                          >  1:1194255/1‑150 (MQ=255)

                            ctgcCGTTCGCGCAGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGAGCTGCTGTTTGCGGTTTTCCTGCAACAGGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTc                                                                                                                          <  1:1616347/150‑1 (MQ=255)

                             tgcCGTTCGCGCAGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTgg                                                                                                                               <  1:2295437/144‑1 (MQ=255)

                             tgcCGTTCGCGCAGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTgg                                                                                                                               >  2:2295437/1‑144 (MQ=255)

                               ccGTTCGCGCAGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGc                                                                                                                       >  2:321233/1‑150 (MQ=255)

                                   tCGCGCAGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTg                                                                                                                   <  1:1815274/150‑1 (MQ=255)

                                   agcagcAGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTg                                                                                                                   <  2:808495/146‑1 (MQ=255)

                                     gcgcAGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGc                                                                                                                 <  2:2396857/150‑1 (MQ=255)

                                       gcAGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCtt                                                                                                               >  2:242919/1‑150 (MQ=255)

                                         aGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCa                                                                                                             <  2:213993/150‑1 (MQ=255)

                                         aGGGCATCCAGCTGCTGGAGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCa                                                                                                             >  1:149024/1‑150 (MQ=255)

                                          gggCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCAc                                                                                                            <  1:321233/150‑1 (MQ=255)

                                          gggCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGATCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCAc                                                                                                            <  2:798428/150‑1 (MQ=255)

                                              aTCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGaa                                                                                                        <  2:317313/150‑1 (MQ=255)

                                              aTCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGaa                                                                                                        >  1:214045/1‑150 (MQ=255)

                                                       ctgGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTc                                                                                               >  1:1586026/1‑150 (MQ=255)

                                                            gcAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAgg                                                                                          >  1:86220/1‑150 (MQ=255)

                                                            gcAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAgg                                                                                          >  1:1269078/1‑150 (MQ=255)

                                                                  ccTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTc                                                                                               <  1:2108569/139‑1 (MQ=255)

                                                                  ccTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTc                                                                                               >  2:2108569/1‑139 (MQ=255)

                                                                    tGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTcgagg                                                                                  <  1:1265010/150‑4 (MQ=255)

                                                                    tGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCCCCGTGCAGGTCGATCTGCTACGACAACCCCTAGTGATCTTCCTTGGCTTCACGGCACAACAGGTCGTAGGCGTCgcgc                                                                                  >  1:535854/1‑150 (MQ=255)

                                                                      ttCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTc                                                                                <  2:1684633/150‑1 (MQ=255)

                                                                       tCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCg                                                                               <  2:149024/150‑1 (MQ=255)

                                                                         aCGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGACCTTgg                                                                                                                  >  2:26824/1‑113 (MQ=255)

                                                                         aCGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGACCTTgg                                                                                                                  <  1:26824/113‑1 (MQ=255)

                                                                           ggCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGaacaa                                                                                                     >  1:1400620/1‑124 (MQ=255)

                                                                           ggCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGaacaa                                                                                                     <  2:1400620/124‑1 (MQ=255)

                                                                           ggCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCt                                                                           >  2:2267431/1‑150 (MQ=255)

                                                                           ggCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCt                                                                           >  2:2268173/1‑150 (MQ=255)

                                                                            gCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTg                                                                          >  2:1761458/1‑150 (MQ=255)

                                                                             ccGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGc                                                                         <  1:2101042/150‑1 (MQ=255)

                                                                             ccGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGc                                                                         <  2:1586026/150‑1 (MQ=255)

                                                                                 ggaGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGGGCAGGTCGATCTGCTCGGCCAGGGCCCGGGGTTCGGCCTTGGCTTCACGGACCAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAgc                                                                     <  2:66721/148‑1 (MQ=255)

                                                                                  cAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAgcg                                                                    >  1:2400324/1‑150 (MQ=255)

                                                                                      cGTGACATTTCCTCTATGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCa                                                                <  1:1100103/150‑1 (MQ=255)

                                                                                       gTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCgg                                                                                                                        >  2:1409318/1‑93 (MQ=255)

                                                                                       gTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCAt                                                               >  1:2237944/1‑150 (MQ=255)

                                                                                       gTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCAt                                                               >  1:2405296/1‑150 (MQ=255)

                                                                                           cATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGc                                                           >  1:440603/1‑150 (MQ=255)

                                                                                           cATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGc                                                           <  1:537675/150‑1 (MQ=255)

                                                                                            aTTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGcc                                                          >  2:2132584/1‑150 (MQ=255)

                                                                                             tttCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCct                                                         >  1:2297255/1‑150 (MQ=255)

                                                                                               tCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCctgg                                                       <  1:2132584/150‑1 (MQ=255)

                                                                                                ccTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCctggc                                                      >  2:2392590/1‑150 (MQ=255)

                                                                                                 ctctACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCAc                                                                                                            >  2:1775867/1‑95 (MQ=255)

                                                                                                 ctctACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCAc                                                                                                            <  1:1775867/95‑1 (MQ=255)

                                                                                                    tACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTgg                                                                                                                  >  2:178906/1‑86 (MQ=255)

                                                                                                    tACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTgg                                                                                                                  <  1:178906/86‑1 (MQ=255)

                                                                                                        gcgcgcTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGc                                                                                                                       <  2:1003460/77‑1 (MQ=255)

                                                                                                        gcgcgcTGCTGTTTGCGGTTGTCCTGCAACATGCAGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGc                                                                                                                       >  1:1003460/1‑77 (MQ=255)

                                                                                                            gctgctGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTggcctgggcct                                          >  2:1428157/1‑150 (MQ=255)

                                                                                                            gctgctGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTggcctgggcct                                          >  1:2228807/1‑150 (MQ=255)

                                                                                                                   tttGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGcc                                   >  1:269918/1‑150 (MQ=255)

                                                                                                                     tGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGTCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGACCAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCgg                                 <  1:1094962/150‑1 (MQ=255)

                                                                                                                      gCGGTTGTTCGGCAATAGGCCGCCGTTCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCggg                                <  2:1957880/150‑1 (MQ=255)

                                                                                                                      gCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTACAGGCCggg                                <  1:242919/150‑1 (MQ=255)

                                                                                                                       cGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGt                               <  2:80775/150‑1 (MQ=255)

                                                                                                                        ggTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTg                                                                                                                                <  1:954661/52‑1 (MQ=255)

                                                                                                                        ggTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAg                                      >  2:2235721/1‑142 (MQ=255)

                                                                                                                        ggTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAg                                      <  1:2235721/142‑1 (MQ=255)

                                                                                                                        ggTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTa                              >  1:1127508/1‑150 (MQ=255)

                                                                                                                         gTTGTCCTGCAACTTGTCGGCGTGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAg                             <  2:2274961/150‑1 (MQ=255)

                                                                                                                              ccTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAGCtctt                        <  1:2185221/150‑1 (MQ=255)

                                                                                                                                tGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAGCtcttct                      >  2:843598/1‑150 (MQ=255)

                                                                                                                                tGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAGCtcttct                      >  2:2371878/1‑150 (MQ=255)

                                                                                                                                 gCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAGCTCTtctc                     <  2:1835709/150‑1 (MQ=255)

                                                                                                                                  cAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAGCTCTTCTCg                    >  1:578042/1‑150 (MQ=255)

                                                                                                                                   aaCATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAg                                      >  1:1864819/1‑131 (MQ=255)

                                                                                                                                   aaCATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAGCTCTTCTCgg                   >  1:1579243/1‑150 (MQ=255)

                                                                                                                                    aCATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAGCTCTTCTCGGc                  >  1:1574535/1‑150 (MQ=255)

                                                                                                                                    aCATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGAGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAGCTCTTCTCGGc                  <  1:1841684/150‑1 (MQ=255)

                                                                                                                                     cATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAGCTCTTCTCGGcc                 <  2:2400324/150‑1 (MQ=255)

                                                                                                                                      aTGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAGCTCTTCTCGGccc                <  1:1398543/150‑1 (MQ=255)

                                                                                                                                        gccgccGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAGCTCTTCTCGGCCCtt              <  2:1269078/150‑1 (MQ=255)

                                                                                                                                        gccgccGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAGCTCTTCTCGGCCCtt              <  2:86220/150‑1 (MQ=255)

                                                                                                                                         tggccgAGCAGATCGACCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAGCTCTTCTCGGCCCTTg             >  2:1649342/3‑150 (MQ=255)

                                                                                                                                         ccgccgAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAGCTCTTCTCGGCCCTTg             <  2:214045/150‑1 (MQ=255)

                                                                                                                                         ccgccgAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAGCTCTTCTCGGCCCTTg             <  2:1540812/150‑1 (MQ=255)

                                                                                                                                          cgccgAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTggcctgggc                                            <  2:186209/118‑1 (MQ=255)

                                                                                                                                          cgccgAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTggcctgggc                                            >  1:186209/1‑118 (MQ=255)

                                                                                                                                                 cAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAGCTCTTCTCGGCCCTTGGTCAGCCg     <  2:2228807/150‑1 (MQ=255)

                                                                                                                                                    gTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAGCTCTTCTCGGCCCTTGGTCAGCCGGcc  >  1:1397165/1‑150 (MQ=255)

                                                                                                                                                    gTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAGCTCTTCTCGGCCCTTGGTCAGCCGGcc  >  1:2076571/1‑150 (MQ=255)

GTTGCTGGCGTTGCAGCATGGCCTGCTGCTGCCGTTCGCGCAGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAGCTCTTCTCGGCCCTTGGTCAGCCGGCC  >  NC_002947/2298499‑2298796

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