BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,984,862	G→A	22.5%	P86P (CCC→CCT)	ycaC‑II ←	putative hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,984,862	0	G	A	22.5%	200.5 / 51.1	102	P86P (CCC→CCT)	ycaC‑II	putative hydrolase
Reads supporting (aligned to +/- strand): ref base G (40/39); new base A (11/12); total (51/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.57e-01

Only 100 of 105 total aligned reads displayed.

CGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGA  >  NC_002947/5984714‑5985010

cGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCg                                                                                                                                                           <  8:4218/149‑1 (MQ=255)

       ttCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACggg                                                                                                                                                    >  8:38458/1‑149 (MQ=255)

       ttCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTAGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACggg                                                                                                                                                    >  7:9352/1‑149 (MQ=255)

           aagccTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAAt                                                                                                                                                >  7:206223/4‑149 (MQ=255)

           aaaCCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACgg                                                                                                                                                     <  5:380333/144‑1 (MQ=255)

           aaaCCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACgg                                                                                                                                                     >  6:380333/1‑144 (MQ=255)

           aaaCCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAAt                                                                                                                                                >  2:355845/1‑149 (MQ=255)

               cTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTAc                                                                                                                                            <  7:263503/149‑1 (MQ=255)

               cTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTAc                                                                                                                                            <  7:126159/149‑1 (MQ=255)

                  ccGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACgg                                                                                                                                          <  7:336195/148‑1 (MQ=255)

                      cAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCAt                                                                                                                                     <  1:355845/149‑1 (MQ=255)

                       aGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATc                                                                                                                                    >  6:251158/1‑149 (MQ=255)

                         cgcCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCg                                                                                                                                  <  4:123113/149‑1 (MQ=255)

                          gcCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCgg                                                                                                                                 <  8:206223/149‑1 (MQ=255)

                             gACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGaa                                                                                                                              >  1:311527/1‑149 (MQ=255)

                              aCAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAAc                                                                                                                             <  3:224220/149‑1 (MQ=255)

                                      gggAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCttt                                                                                                                     <  8:272183/149‑1 (MQ=255)

                                        gAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTCGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCa                                                                                                                   <  8:9352/149‑1 (MQ=255)

                                         aaGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAg                                                                                                                  >  1:39401/1‑149 (MQ=255)

                                            gCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTc                                                                                                               >  5:125406/1‑149 (MQ=255)

                                            gCCACGCACACATCGGTCACCACACCGGCGAAGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTc                                                                                                               >  1:42708/1‑149 (MQ=255)

                                             ccACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCg                                                                                                              <  5:251158/149‑1 (MQ=255)

                                              cACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCgg                                                                                                             >  1:258316/1‑149 (MQ=255)

                                               aCGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCATGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCggg                                                                                                            <  5:31321/149‑1 (MQ=255)

                                                  cacacaTCGGTCACCACTCCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CAc                                                                                                      <  8:81594/149‑1 (MQ=255)

                                                    cacaTCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTCACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCa                                                                                                    <  4:373132/149‑1 (MQ=255)

                                                     acaTCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCa                                                                                                                   <  6:333138/136‑1 (MQ=255)

                                                     acaTCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATCTACGGCGCATCCGGGAACATCTCTTTCa                                                                                                                   >  5:333138/1‑136 (MQ=255)

                                                         cGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCgg‑g                                                                                             >  6:305306/1‑149 (MQ=255)

                                                           gTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTc                                                                                                                    >  5:364263/1‑129 (MQ=255)

                                                           gTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTc                                                                                                                    <  6:364263/129‑1 (MQ=255)

                                                              accacACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCAtt                                                                                        >  3:381649/1‑149 (MQ=255)

                                                                cacaCCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCg                                                                                      >  3:52780/1‑149 (MQ=255)

                                                                 acacCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCgg                                                                                     <  2:42708/149‑1 (MQ=255)

                                                                     cGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCt                                                                                 >  3:205370/1‑149 (MQ=255)

                                                                     cGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCt                                                                                 >  5:140999/1‑149 (MQ=255)

                                                                     cGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCt                                                                                 >  4:35237/1‑149 (MQ=255)

                                                                        cGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGc                                                                              <  7:161544/149‑1 (MQ=255)

                                                                        cGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGc                                                                              <  8:27213/149‑1 (MQ=255)

                                                                         gatgatCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCt                                                                             >  6:260696/1‑149 (MQ=255)

                                                                             atCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCg                                                                           >  5:182692/1‑147 (MQ=255)

                                                                              tCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAg                                                                        >  4:261754/1‑149 (MQ=255)

                                                                              tCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAg                                                                        >  4:266994/1‑149 (MQ=255)

                                                                                aGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCt                                                                      >  6:167555/1‑149 (MQ=255)

                                                                                 gTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCtg                                                                     >  3:10022/1‑149 (MQ=255)

                                                                                 gTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCtg                                                                     >  2:368668/1‑149 (MQ=255)

                                                                                 gTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCtg                                                                     >  4:107213/1‑149 (MQ=255)

                                                                                     cttgcGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCtggtgg                                                                 <  5:305306/149‑1 (MQ=255)

                                                                                      ttgcGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTggaggt                                                                >  7:208326/1‑149 (MQ=255)

                                                                                         cGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGa                                                              <  2:2631/148‑1 (MQ=255)

                                                                                         cGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGa                                                              <  4:205370/148‑1 (MQ=255)

                                                                                          ggCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAgg                                                            <  4:381649/149‑1 (MQ=255)

                                                                                          ggCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAg                                                             >  1:334892/1‑148 (MQ=255)

                                                                                           gCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGa                                                           <  1:368668/149‑1 (MQ=255)

                                                                                            cccGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGAt                                                          <  5:260696/149‑1 (MQ=255)

                                                                                              cGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATgg                                                        >  7:315786/1‑149 (MQ=255)

                                                                                                      tcaccgCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCtgg                                                                    >  1:79639/6‑129 (MQ=11)

                                                                                                      tcaccgCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCtgg                                                                    <  2:79639/124‑1 (MQ=11)

                                                                                                       gATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAg                                                 <  7:38458/147‑1 (MQ=255)

                                                                                                       gATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAg                                                 <  2:39401/147‑1 (MQ=255)

                                                                                                        accgCCTAGACGAAATCCTCCTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGcc                                              <  8:261008/146‑1 (MQ=11)

                                                                                                        aTTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCACGCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTc                                                      <  1:81514/142‑1 (MQ=255)

                                                                                                        aTTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTc                                                      >  2:81514/1‑141 (MQ=255)

                                                                                                            ccTTGACGAAATCCTCGTTGTCCCCGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAg                                          <  6:214536/149‑1 (MQ=14)

                                                                                                             cTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAg                                          >  3:270087/1‑148 (MQ=255)

                                                                                                                gACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGc                                                                                   >  4:53795/1‑104 (MQ=14)

                                                                                                                gACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGc                                                                                   <  3:53795/104‑1 (MQ=14)

                                                                                                                 aCGAACGCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACtt                                     <  2:152754/149‑1 (MQ=255)

                                                                                                                  cGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCGTGc                                                                              >  4:330258/1‑107 (MQ=255)

                                                                                                                  cGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTg                                    <  2:52431/149‑1 (MQ=255)

                                                                                                                   gAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTGGTGTTAGAGGA‑CCGGAGGCAATAGGGCCTTGATCGAAGCTGGTGGTGAGGATGTTCGGCAGGCCAAAGAAc                                       >  1:324330/1‑149 (MQ=18)

                                                                                                                   gAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGt                                                       >  5:100416/1‑129 (MQ=255)

                                                                                                                   gAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGt                                                       <  6:100416/129‑1 (MQ=255)

                                                                                                                    aaGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGc                                  <  6:125406/149‑1 (MQ=255)

                                                                                                                     aGTCTTCGTTGTCACAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGc                                  >  5:129721/1‑148 (MQ=255)

                                                                                                                      atcctCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGc                                  >  5:120010/5‑147 (MQ=255)

                                                                                                                      atcctCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCa                                <  7:69621/145‑1 (MQ=255)

                                                                                                                        cctCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAgg                                                <  4:180980/129‑1 (MQ=255)

                                                                                                                        cctCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAgg                                                >  3:180980/3‑131 (MQ=255)

                                                                                                                            gTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTc                            <  5:324430/147‑1 (MQ=255)

                                                                                                                             ttGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGcc                         >  1:149483/1‑149 (MQ=255)

                                                                                                                             ttGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGcc                         >  3:20122/1‑149 (MQ=255)

                                                                                                                               gTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGccc                        <  2:109564/148‑1 (MQ=255)

                                                                                                                               gTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGccc                        >  1:7861/1‑148 (MQ=255)

                                                                                                                                  ccAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAgg                    >  3:48654/1‑149 (MQ=255)

                                                                                                                                  ccAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAgg                    <  2:311527/149‑1 (MQ=255)

                                                                                                                                  ccAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAgg                    >  6:291486/1‑149 (MQ=255)

                                                                                                                                  ccAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAgg                    >  7:194435/1‑149 (MQ=255)

                                                                                                                                  ccAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAg                     <  4:152974/148‑1 (MQ=255)

                                                                                                                                  ccAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAg                     <  2:351151/148‑1 (MQ=255)

                                                                                                                                   cAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGc                   >  6:372938/1‑149 (MQ=255)

                                                                                                                                   cAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGc                   >  3:75535/1‑149 (MQ=255)

                                                                                                                                    aGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGcc                  <  6:140999/149‑1 (MQ=255)

                                                                                                                                       cGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGc               <  6:182692/149‑1 (MQ=255)

                                                                                                                                         ttGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCAc             <  5:366760/149‑1 (MQ=255)

                                                                                                                                         ttGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCAc             <  5:372920/149‑1 (MQ=255)

                                                                                                                                           gATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCATCACgt           >  5:211444/1‑149 (MQ=255)

                                                                                                                                            aTCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACgt           <  7:320704/148‑1 (MQ=255)

                                                                                                                                                 gCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACgttg         <  6:129721/145‑1 (MQ=255)

                                                                                                                                                    gggACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCTGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGTACGTTGTTCTTGa  <  1:351950/149‑1 (MQ=255)

CGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGA  >  NC_002947/5984714‑5985010

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