BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,836,768	G→A	9.6%	R219Q (CGG→CAG)	yfiP →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,836,768	0	G	A	9.6%	755.2 / 49.7	260	R219Q (CGG→CAG)	yfiP	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (97/138); new base A (8/17); total (105/155)
Fisher's exact test for biased strand distribution p-value = 4.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 278 total aligned reads displayed.

CCGCCTGCGGCGCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTCCCTTGA  >  NC_002947/1836619‑1836916

ccGCCTGCGGCGCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAgcggcg                                                                                                                                                      <  2:3756766/150‑1 (MQ=255)

 cGCCTGCGGCGCTCCCAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGcggcgg                                                                                                                                                     <  2:1565202/150‑1 (MQ=255)

 cGCCTGCGGCGCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGcggcgg                                                                                                                                                     >  2:4676998/1‑150 (MQ=255)

     tGCGGCGCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGCTCGTCAAGCGGCAGCCGc                                                                                                                                                 <  1:6312611/150‑1 (MQ=255)

        ggcgCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGc                                                                                                                                              >  2:448109/1‑150 (MQ=255)

          cgcTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGTCAAGCGGCAGCCGCCGCtt                                                                                                                                            <  1:6083019/150‑1 (MQ=255)

           gcTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCAGCCGCCGCTTg                                                                                                                                           <  1:1766864/150‑1 (MQ=255)

           gcTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCAGCCGCCGCTTg                                                                                                                                           <  2:3050808/150‑1 (MQ=255)

            cTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCAGCCGCCGCTTgc                                                                                                                                          <  1:2453230/150‑2 (MQ=255)

             tCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGAc                                                                                                                                         <  1:3085666/150‑1 (MQ=255)

             tCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGAc                                                                                                                                         <  1:5273205/150‑1 (MQ=255)

                aaGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAg                                                                                                                                      <  2:2875295/150‑1 (MQ=255)

                  gcgcGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGCTCGTCAAGCGGCAGCCGCCGCTTGCCgctca                                                                                                                                    <  2:1313562/150‑5 (MQ=39)

                  gcgcGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCAGCCGCCGCTTGCCgctca                                                                                                                                    <  1:2025143/150‑5 (MQ=255)

                   cgcgATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGCTCGTCAAGCGGCAGCCGCCGCTTGCCgctcag                                                                                                                                   <  2:2093265/150‑6 (MQ=38)

                    gcgTTGATCATTTCTGTACGGCAGAAGTTGCGGCGAGGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCTCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGATCAgg                                                                                                                                  <  1:3591504/150‑1 (MQ=255)

                    gcgATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAgg                                                                                                                                  >  1:2647100/1‑150 (MQ=255)

                    gcgATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAgg                                                                                                                                  <  1:3262259/150‑1 (MQ=255)

                     cgATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGga                                                                                                                                 <  1:1865593/150‑1 (MQ=255)

                     cgATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGga                                                                                                                                 <  1:3535706/150‑1 (MQ=255)

                      gatgatCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGgac                                                                                                                                >  1:2494411/1‑150 (MQ=255)

                      gatgatCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGgac                                                                                                                                >  1:898911/1‑150 (MQ=255)

                         gatCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGgacgac                                                                                                                             <  2:6271635/150‑1 (MQ=255)

                         gatCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGgacgac                                                                                                                             >  2:5423890/1‑150 (MQ=255)

                              tttCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGAcaccca                                                                                                                        <  2:553007/150‑1 (MQ=255)

                               ttCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGAcacccac                                                                                                                       <  2:644610/150‑1 (MQ=255)

                               ttCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGAcacccac                                                                                                                       <  2:2735006/150‑1 (MQ=255)

                                  tGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAg                                                                                                                    <  2:2362557/150‑1 (MQ=255)

                                       ggCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGc                                                                                                                                                 <  1:2359627/116‑1 (MQ=255)

                                       ggCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGttt                                                                                                               <  1:6449167/150‑1 (MQ=255)

                                       ggCAGAAGGTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGCTGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGttt                                                                                                               <  1:738242/150‑1 (MQ=255)

                                        gCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTg                                                                                                              <  2:2231778/150‑1 (MQ=255)

                                        gCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTg                                                                                                              <  1:5637966/150‑1 (MQ=255)

                                        gCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTg                                                                                                              <  2:4249034/150‑1 (MQ=255)

                                        gCAGAAGTCGCGGCGCCGTGCCTGGACCGTGCTGGCGAGGGCCAGGCTTCGCAGGCGCTGGATGCCGATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTg                                                                                                              <  2:5744779/150‑1 (MQ=255)

                                          aGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCa                                                                                                            >  2:3377296/1‑150 (MQ=255)

                                            aaGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATa                                                                                                          <  1:5456508/150‑1 (MQ=255)

                                             aGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATAc                                                                                                         <  2:41818/150‑1 (MQ=255)

                                                 gcggcgATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCacctggtc                                                                                                     >  1:1211823/1‑143 (MQ=255)

                                                 gcggcgATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTc                                                                                                     <  2:6166807/150‑1 (MQ=255)

                                                      gATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCTGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCtatc                                                                                                <  1:3699893/150‑2 (MQ=255)

                                                      gATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCtata                                                                                                <  1:6759105/150‑1 (MQ=255)

                                                                gACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAgggg                                                                                      <  2:2558580/150‑1 (MQ=255)

                                                                    ttGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGcaca                                                                                  >  2:5022129/1‑150 (MQ=255)

                                                                        tGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCg                                                                                                                                               <  2:6388251/85‑1 (MQ=255)

                                                                        tGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAAccc                                                                              >  1:4152580/1‑150 (MQ=255)

                                                                        tGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAAccc                                                                              >  2:5196505/1‑150 (MQ=255)

                                                                         ggCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAAcccc                                                                             >  2:476260/1‑150 (MQ=255)

                                                                          gCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCt                                                                            <  2:4139109/150‑1 (MQ=255)

                                                                                gCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACgctca                                                                                                                                    >  1:2674824/1‑84 (MQ=255)

                                                                                 ccAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTgg                                                                     <  2:6639091/150‑1 (MQ=255)

                                                                                  cAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGt                                                                    <  2:2514708/150‑1 (MQ=255)

                                                                                     gCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCa                                                                 <  2:6737196/150‑1 (MQ=255)

                                                                                       ttCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATa                                                               >  2:527712/1‑150 (MQ=255)

                                                                                       ttCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATa                                                               >  1:3173405/1‑150 (MQ=255)

                                                                                        tCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATaa                                                              <  1:430635/150‑1 (MQ=255)

                                                                                         cGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAAt                                                             <  2:4376055/150‑1 (MQ=255)

                                                                                          gCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCGGCTTTGGTTCATAATg                                                            >  2:1511316/1‑150 (MQ=255)

                                                                                           cAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGa                                                           >  1:1899958/1‑150 (MQ=255)

                                                                                               cgcTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGa                                                       <  2:5380557/150‑1 (MQ=255)

                                                                                                   ggATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTg                                                   >  2:1525794/1‑150 (MQ=255)

                                                                                                    gATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTgg                                                  >  1:458754/1‑150 (MQ=255)

                                                                                                           aTCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAg                                           <  2:6154182/150‑1 (MQ=255)

                                                                                                            tCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACACACCACCGTTTGCATACCTACCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAgg                                          >  1:4003412/1‑150 (MQ=255)

                                                                                                             cTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCCAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAggg                                         <  1:5413580/150‑1 (MQ=255)

                                                                                                             cTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAggg                                         <  1:4676998/150‑1 (MQ=255)

                                                                                                                 aTGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAgcaggat                                     <  2:4847573/150‑7 (MQ=255)

                                                                                                                  tgtgTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGc                                    <  1:4387750/150‑1 (MQ=255)

                                                                                                                    tgtTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCga                                  <  2:1864553/150‑1 (MQ=255)

                                                                                                                      ttCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCCATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCgaga                                >  2:3441481/1‑150 (MQ=255)

                                                                                                                         aGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACgg                             <  1:1551614/150‑1 (MQ=255)

                                                                                                                           ccTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCa                           >  1:1291601/1‑150 (MQ=255)

                                                                                                                             tGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAgg                         <  2:1099025/150‑1 (MQ=255)

                                                                                                                             tGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAgg                         >  1:1151901/1‑150 (MQ=255)

                                                                                                                              gCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAggg                        >  1:2302910/1‑150 (MQ=255)

                                                                                                                               cATGACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTAGAGTCCCAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGc                       <  2:4545257/150‑1 (MQ=255)

                                                                                                                                aTTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCg                      <  2:2207335/150‑1 (MQ=255)

                                                                                                                                 ttACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGa                     >  1:4624758/1‑150 (MQ=255)

                                                                                                                                   aCCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGAtt                   <  1:1788425/150‑1 (MQ=255)

                                                                                                                                   aCCCGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGAtt                   <  2:5661056/150‑1 (MQ=255)

                                                                                                                                      tgctcgtcAAGCGGCAGCCGCCGCTTGCCGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCGGCTGTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGAGTGGc                >  2:4042901/8‑150 (MQ=255)

                                                                                                                                      tGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGc                >  1:6660154/1‑150 (MQ=255)

                                                                                                                                      tGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGc                >  1:3608827/1‑150 (MQ=255)

                                                                                                                                       gAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATa                                                               <  1:5921667/102‑1 (MQ=255)

                                                                                                                                       gAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATa                                                               >  2:5921667/1‑102 (MQ=255)

                                                                                                                                       gAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCt               >  1:3552961/1‑150 (MQ=255)

                                                                                                                                        aGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCtt              <  1:4939310/150‑1 (MQ=255)

                                                                                                                                         tcgtcAAGCGGCAGCCGCCGCTTGCCGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTg             >  1:2548906/5‑150 (MQ=255)

                                                                                                                                          cGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGa            <  1:965740/150‑1 (MQ=255)

                                                                                                                                          cGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGa            <  1:3383808/150‑1 (MQ=255)

                                                                                                                                           gtcAAGCGGCAGCCGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGAt           >  2:1943496/3‑150 (MQ=255)

                                                                                                                                           gtcAAGCGGCAGCCGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGAt           >  1:4462304/3‑150 (MQ=255)

                                                                                                                                              aaGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGAtctc        >  2:481761/1‑150 (MQ=255)

                                                                                                                                              aaGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGAtctc        <  1:3801160/150‑1 (MQ=255)

                                                                                                                                                gcggcgGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTccc      <  1:666810/150‑1 (MQ=255)

                                                                                                                                                gcggcgGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTccc      <  1:481761/150‑1 (MQ=255)

                                                                                                                                                gcggcgGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTACATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTccc      <  2:4115663/150‑1 (MQ=255)

                                                                                                                                                 cggcggCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTCCCt     >  2:1194045/1‑150 (MQ=255)

                                                                                                                                                  ggcggcTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCt                                                    <  2:688018/102‑1 (MQ=255)

                                                                                                                                                    cggcTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTCCCTTGa  >  1:5246058/1‑150 (MQ=255)

CCGCCTGCGGCGCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTCCCTTGA  >  NC_002947/1836619‑1836916

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