BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,992,312	A→G	7.4%	I111I (ATT→ATC)	ycaC‑III ←	putative hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,992,312	0	A	G	7.4%	1101.8 / 46.5	340	I111I (ATT→ATC)	ycaC‑III	putative hydrolase
Reads supporting (aligned to +/- strand): ref base A (195/118); new base G (17/8); total (214/126)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

Only 100 of 349 total aligned reads displayed.

GCACCCGCCTGGCTCATGCGGTCATGCGCGGCATCACGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTC  >  NC_002947/5992163‑5992461

gCACCCGCCTGGCTCATGCGGTCATGCGCGGCATCACGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCa                                                                                                                                                       >  2:6173304/1‑150 (MQ=255)

atgaccGCCTGGCTCATGCGGTCATGCGCGGCATCACGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCa                                                                                                                                                       >  1:3897171/5‑150 (MQ=255)

 cACCCGCCGGGCTCATGCGGTCATGTGCGGCATCACGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCaa                                                                                                                                                      <  2:7355334/150‑1 (MQ=255)

    ccGCCTGGCTCATGCGGTCATGCGCGGCATCACGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGa                                                                                                                                                   >  2:3236643/1‑150 (MQ=255)

     cGCCTGGCTCATGCGGTCATGCGCGGCATCACGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGAt                                                                                                                                                  <  1:8656620/150‑1 (MQ=255)

      gCCTGGCTCATGCGGTCATGCGCGGCATCACGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATc                                                                                                                                                 >  2:7073403/1‑150 (MQ=255)

       ccTGGCTCATGCGGTCATGCGCGGCATCACGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCa                                                                                                                                                >  1:2630422/1‑150 (MQ=255)

            cTCATGCGGTCATGCGCGGCATCACGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAgctgc                                                                                                                                           <  1:5082389/150‑1 (MQ=255)

            cTCATGCGGTCATGCGCGGCATCACGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCAGCAATGATCAgctgc                                                                                                                                           >  1:1083690/1‑150 (MQ=255)

             tCATGCGGTCATGCGCGGCATCACGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAgctgct                                                                                                                                          >  2:9089552/1‑150 (MQ=255)

                tGCGGTCATGCGCGGCATCACGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGcttct                                                                                                                                       >  1:3684844/1‑150 (MQ=255)

                   ggTCATGCGCGGCATCACGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCACACCTCGGTCACAACACCTGCCATGATCAGCTGCTTCTTGc                                                                                                                                    >  1:9059442/1‑150 (MQ=255)

                    gTCATGCGCGGCATCACGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGcc                                                                                                                                   >  1:246915/1‑150 (MQ=255)

                     tCATGCGCGGCATCACGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCg                                                                                                                                  <  1:2650863/150‑1 (MQ=255)

                     tCATGCGCGGCATCACGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCg                                                                                                                                  >  1:8464477/1‑150 (MQ=255)

                       aTGCGCGGCATCACGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGt                                                                                                                                >  2:1135372/1‑150 (MQ=255)

                        tGCGCGGCATCACGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTc                                                                                                                               >  1:4539442/1‑150 (MQ=255)

                         gcgcGGCATCACGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCg                                                                                                                              >  1:1120508/1‑150 (MQ=255)

                          cgcgGCATCACGGGTAATTTCATTGAAGGTGCCGGAAGCATCGTTCACCACGAACACCTCAAATTCTTACTCCAAGGACGACAGCGCAGGGAACGGCACGCAAAACTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGAAGGTCGc                                                                                                                             >  2:2356479/1‑150 (MQ=255)

                           gaggCATCACGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCATACCTCGGTCACAACACCTGCAATGATCAGCTTCTTCTTGCCGGTCGcc                                                                                                                            <  1:8857295/148‑1 (MQ=255)

                                   aCGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGc                                                                                                                    >  2:3467338/1‑150 (MQ=255)

                                    cGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGcc                                                                                                                   >  2:1680270/1‑150 (MQ=255)

                                    cGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCACCTTCACCGcc                                                                                                                   >  1:4687433/1‑150 (MQ=255)

                                     gggTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCt                                                                                                                  <  2:3122809/150‑1 (MQ=255)

                                       gTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACTCCTCAAATTCTTCCTCCAAGACCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTg                                                                                                                <  2:6873848/150‑1 (MQ=255)

                                       gTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGACCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTg                                                                                                                <  2:2254427/150‑1 (MQ=255)

                                        tcatttcattGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGa                                                                                                               <  2:596763/150‑1 (MQ=255)

                                           tttcattGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGa                                                                                                            >  1:7938477/1‑150 (MQ=255)

                                            ttcattGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGaa                                                                                                           >  2:79438/1‑150 (MQ=255)

                                            ttcattGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGaa                                                                                                           <  1:984474/150‑1 (MQ=255)

                                             tcattGAAGGTGCCGGAAGCATCGGTCACCACGAACACATCAAATTCTTCCTCAAAGGCCGAAAGCGCAGGGAACGCCACGCAAACCACGGACACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTACACCGCCTTGACGaaa                                                                                                          >  2:4644218/1‑150 (MQ=255)

                                                  gAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTc                                                                                                     >  2:5110605/1‑150 (MQ=255)

                                                    aGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGt                                                                                                   <  1:1812758/150‑1 (MQ=255)

                                                     ggTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGtt                                                                                                  >  2:2696381/1‑150 (MQ=255)

                                                           ggAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCa                                                                                            <  1:4736710/150‑1 (MQ=255)

                                                           ggAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGACACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCATCGTTGTCCCa                                                                                            >  1:3010092/1‑150 (MQ=255)

                                                            gAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAg                                                                                           <  2:8904280/150‑1 (MQ=255)

                                                                cATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGt                                                                                       >  2:1590173/1‑150 (MQ=255)

                                                                 aTCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGtt                                                                                      <  2:3783030/150‑1 (MQ=255)

                                                                  tCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTg                                                                                     >  2:8316988/1‑150 (MQ=255)

                                                                   cGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGa                                                                                    <  1:1050486/150‑1 (MQ=255)

                                                                         ccacGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGc                                                                              <  1:5370302/150‑1 (MQ=255)

                                                                            cGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAg                                                                           <  1:4008846/150‑1 (MQ=255)

                                                                               acacCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAgggc                                                                        >  1:6241266/1‑150 (MQ=255)

                                                                                  ccTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAgggcggg                                                                     >  1:8186996/1‑150 (MQ=255)

                                                                                  ccTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAgggcggg                                                                     >  2:1868839/1‑150 (MQ=255)

                                                                                    tCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCa                                                                   >  1:3295636/1‑150 (MQ=255)

                                                                                    tCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCa                                                                   >  1:7399074/1‑150 (MQ=255)

                                                                                    tCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCa                                                                   >  2:5567330/1‑150 (MQ=255)

                                                                                     cAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCaa                                                                  >  1:5595762/1‑150 (MQ=255)

                                                                                      aaaTTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAAt                                                                 <  2:1986680/150‑1 (MQ=255)

                                                                                      aaaTTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTAACAACACCTGCAATGATCAGCTGATTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCATGCGATGATCTGGCCAGGGCGGGCAAt                                                                 >  1:2235041/1‑150 (MQ=255)

                                                                                         ttcttcCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTc                                                                                                     >  1:7345971/1‑111 (MQ=255)

                                                                                          tcttcCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTAc                                                             <  1:2829686/150‑1 (MQ=255)

                                                                                                tccAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCa                                                       >  2:4115336/1‑150 (MQ=255)

                                                                                                   aaGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATcc                                                    <  1:1995893/150‑1 (MQ=255)

                                                                                                        ggacaatgcCGGGAAAGCCACGCAGACTTCGGTCACAACACCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGaa                                               >  2:1033417/8‑150 (MQ=14)

                                                                                                        cGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACAGCGCATCCGGGaa                                               <  2:2482433/150‑1 (MQ=255)

                                                                                                          aCAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACa                                             <  1:6644448/150‑1 (MQ=255)

                                                                                                          aCAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACa                                             <  1:151082/150‑1 (MQ=255)

                                                                                                           cAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAg                                            <  1:6394607/150‑1 (MQ=255)

                                                                                                              tgcCGGGAAAGCCACGCAGACTTCGGTCACAACACCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGaatgtttc                                         >  2:1372450/2‑144 (MQ=14)

                                                                                                              cgcCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGc                                         <  1:755705/150‑1 (MQ=255)

                                                                                                              cgcCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGc                                         <  1:1323111/150‑1 (MQ=255)

                                                                                                                ccGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCtt                                       <  1:5889407/150‑1 (MQ=255)

                                                                                                                 cGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCttt                                      <  2:2458386/150‑1 (MQ=255)

                                                                                                                 cGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCttt                                      <  2:7229895/150‑1 (MQ=255)

                                                                                                                  gggAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTc                                     >  1:2101246/1‑150 (MQ=255)

                                                                                                                   ggAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCa                                    >  1:6540646/1‑150 (MQ=255)

                                                                                                                   ggAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCa                                    <  2:8649672/150‑1 (MQ=255)

                                                                                                                   ggAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCa                                    <  2:7759591/150‑1 (MQ=255)

                                                                                                                   ggAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCa                                    <  2:2822733/150‑1 (MQ=255)

                                                                                                                        gCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTcc                               <  1:7701029/150‑1 (MQ=255)

                                                                                                                         ccACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCg                              >  1:751069/1‑150 (MQ=255)

                                                                                                                            cGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCa                           >  2:660102/1‑150 (MQ=255)

                                                                                                                             gCAGACTTCGGTCACAACACCTGCGATGAGCAACTGCGTCTTGCCGGTGGCCTTCCCTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCCTCCGGGAATGTTTCTTTCAGCTCCGGCAc                          >  2:3345316/1‑150 (MQ=9)

                                                                                                                                  cttCGGTCACAACACCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGAATGTTTCTTTCAGCTCCGGCACCaatg                     <  2:8137616/148‑4 (MQ=9)

                                                                                                                                   cTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCgg                    <  2:6699375/150‑1 (MQ=255)

                                                                                                                                    tCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCggg                   >  2:4166556/1‑150 (MQ=255)

                                                                                                                                      ggTCACAACACCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGAATGTTTCTTTCAGCTCCGGCACCAATGGGcc                 >  2:4879863/1‑150 (MQ=9)

                                                                                                                                       gTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCACTGTACGGCGCATCCGGGAACAGGGCTGTCAGTTCCGGCACCAGCGGGCCg                >  1:3950826/1‑150 (MQ=255)

                                                                                                                                       gTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGAGCTGGCCAGTGCGGGCAATGTACGGCGCATCCGGGACCAGGGCTTTCAGTTCCGGCACCAGCGGGCCg                >  2:7510511/1‑150 (MQ=255)

                                                                                                                                       gTCACAACACCTGCAATGATCAGCTGCGTCTTGCCGGTCGCCGTCCCCGCCTTGACGAAATCCTCGGTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCg                <  2:761485/150‑1 (MQ=255)

                                                                                                                                        tCACAACACCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGAATGTTTCTTTCAGCTCCGGCACCAATGGGccat               >  1:1608516/1‑148 (MQ=9)

                                                                                                                                          acaacaCCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGAATGTTTCTTTCAGCTCCGGCACCAATGGGCCAttc             >  1:830551/1‑149 (MQ=2)

                                                                                                                                          acaacaCCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTg             <  2:5370899/150‑1 (MQ=255)

                                                                                                                                          acaacaCCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTg             >  2:6974149/1‑150 (MQ=255)

                                                                                                                                           caacaCCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGACCAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTgg            >  2:6114665/1‑150 (MQ=255)

                                                                                                                                            aacaCCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGAATGTTTCTTTCAGCTCCGGCACCAATGGGCCATTcgg           >  1:5655282/1‑150 (MQ=2)

                                                                                                                                            aacaCCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTggg           >  2:1449256/1‑150 (MQ=255)

                                                                                                                                            aacaCCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTggg           <  1:7073403/150‑1 (MQ=255)

                                                                                                                                            aacaCCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTggg           >  1:692079/1‑150 (MQ=255)

                                                                                                                                            aacaCCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTggg           >  1:5491882/1‑150 (MQ=255)

                                                                                                                                            aacaCCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTggg           >  2:4317169/1‑150 (MQ=255)

                                                                                                                                               acCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGcc        <  1:6074111/150‑1 (MQ=255)

                                                                                                                                                  tGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGAATGTTTCTTTCAGCTCCGGCACCAATGGGCCATTCGGGccgtt     >  2:6249996/1‑147 (MQ=2)

                                                                                                                                                   gCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGc    <  2:3684844/150‑1 (MQ=255)

                                                                                                                                                   gCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGc    <  1:4477822/150‑1 (MQ=255)

                                                                                                                                                    cAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCt   >  2:8788004/1‑150 (MQ=255)

                                                                                                                                                     aaTGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCACTGTACGGCGCATCCGGGAACAGGGCTGTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTc  >  2:9012743/1‑150 (MQ=255)

GCACCCGCCTGGCTCATGCGGTCATGCGCGGCATCACGGGTCATTTCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTC  >  NC_002947/5992163‑5992461

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: `‑`