BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,803,269	G→C	5.3%	R871R (CGC→CGG)	PP_3361 ←	aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,803,269	0	G	C	5.3%	1348.5 / 45.2	396	R871R (CGC→CGG)	PP_3361	aminotransferase
Reads supporting (aligned to +/- strand): ref base G (190/185); new base C (11/10); total (201/195)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

Only 100 of 413 total aligned reads displayed.

AGCAGGTAATGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAG  >  NC_002947/3803120‑3803418

aGCAGGTAATGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTcg                                                                                                                                                       >  1:6062576/1‑150 (MQ=255)

 gCAGGTAATGGCCGACACTCAAGGCGTTTTCTTTCAACTTGTCGCGCTCGTTTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGATAGC‑TTCGGCAACCTcgc                                                                                                                                                      <  2:3777376/150‑1 (MQ=255)

 gCAGGTAATGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTcgc                                                                                                                                                      >  2:7242642/1‑150 (MQ=255)

     gTAATGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGt                                                                                                                                                  >  1:1903074/1‑150 (MQ=255)

     gTAATGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGt                                                                                                                                                  <  2:6351469/150‑1 (MQ=255)

      tAATGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTg                                                                                                                                                 >  2:3555773/1‑150 (MQ=255)

        aTGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGt                                                                                                                                               >  2:7391081/1‑150 (MQ=255)

         tGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTc                                                                                                                                              >  2:4280781/1‑150 (MQ=255)

                 acTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTAccgc                                                                                                                                      >  1:1466595/1‑150 (MQ=255)

                   tCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCCCTAccgccg                                                                                                                                    >  1:8720526/1‑150 (MQ=255)

                   tCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTAccgccg                                                                                                                                    >  1:655884/1‑150 (MQ=255)

                   tCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTAccgccg                                                                                                                                    >  2:4630371/1‑150 (MQ=255)

                       ggCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACAt                                                                                                                                >  1:6284047/1‑150 (MQ=255)

                          gTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCgg                                                                                                                             <  1:678251/150‑1 (MQ=255)

                          gTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCgg                                                                                                                             <  2:2185678/150‑1 (MQ=255)

                          gTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCgg                                                                                                                             <  1:5394117/150‑1 (MQ=255)

                           ttttCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCggg                                                                                                                            >  2:8716696/1‑150 (MQ=255)

                            tttctttcAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGt                                                                                                                           >  1:6413564/1‑150 (MQ=255)

                               ctttcAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGc                                                                                                                        <  2:8363376/150‑1 (MQ=255)

                                tttcAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGcc                                                                                                                       >  2:2681284/1‑150 (MQ=255)

                                 ttcAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCg                                                                                                                      >  2:4878148/1‑150 (MQ=255)

                                     aCTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTAc                                                                                                                  >  1:4414077/1‑150 (MQ=255)

                                     aCTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTAc                                                                                                                  <  1:7215238/150‑1 (MQ=255)

                                     aCTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTAc                                                                                                                  <  2:7341622/150‑1 (MQ=255)

                                        ggTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCaa                                                                                                               <  2:1013796/150‑1 (MQ=255)

                                              gcTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGtt                                                                                                         >  2:1150504/1‑150 (MQ=255)

                                                 cGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGc                                                                                                      <  1:1580233/150‑1 (MQ=255)

                                                 cGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATCCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGc                                                                                                      <  1:2864857/150‑1 (MQ=255)

                                                  gatgcaaTCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGcc                                                                                                     >  2:2813634/7‑150 (MQ=255)

                                                  gATTGCATCGAGTACCGCCAGGCCTACCGCGCGTGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGcc                                                                                                     >  1:6441347/1‑150 (MQ=255)

                                                     tGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCAt                                                                                                  <  1:3057219/150‑1 (MQ=255)

                                                      gCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATg                                                                                                 <  2:3979013/150‑1 (MQ=255)

                                                      gCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATg                                                                                                 <  2:3375513/150‑1 (MQ=255)

                                                      gCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATg                                                                                                 <  2:2926901/150‑1 (MQ=255)

                                                       cATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGc                                                                                                <  1:897680/150‑1 (MQ=255)

                                                        aTCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCt                                                                                               >  1:8353514/1‑150 (MQ=255)

                                                         tCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTg                                                                                              >  1:5793156/1‑150 (MQ=255)

                                                         tCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTg                                                                                              >  2:320423/1‑150 (MQ=255)

                                                         tCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTg                                                                                              >  2:8476877/1‑150 (MQ=255)

                                                           gAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGAc                                                                                            <  1:1255677/150‑1 (MQ=255)

                                                           gAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGAc                                                                                            <  2:5336902/150‑1 (MQ=255)

                                                           gAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGAc                                                                                            <  2:6831933/150‑1 (MQ=255)

                                                           gAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGAc                                                                                            <  1:3774651/150‑1 (MQ=255)

                                                                ccGCCAGGCCGACCGCGCATGACACCGGGTGACCGGCGACGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATACTGACGGCAt                                                                                       <  2:6500830/150‑1 (MQ=255)

                                                                 cGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATc                                                                                      >  2:8215116/1‑150 (MQ=255)

                                                                 cGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATc                                                                                      >  1:8600200/1‑150 (MQ=255)

                                                                   ccAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCg                                                                                    >  1:5144622/1‑150 (MQ=255)

                                                                   ccAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGCATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCg                                                                                    >  2:8655649/1‑150 (MQ=255)

                                                                    cAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGt                                                                                                                     <  1:2737038/116‑1 (MQ=255)

                                                                      ggCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGa                                                                                 <  1:5688352/150‑1 (MQ=255)

                                                                      ggCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGa                                                                                 <  1:5338418/150‑1 (MQ=255)

                                                                       gCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGAc                                                                                >  2:7808253/1‑150 (MQ=255)

                                                                       gCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGAc                                                                                <  1:8546073/150‑1 (MQ=255)

                                                                       gCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGAc                                                                                <  2:1627356/150‑1 (MQ=255)

                                                                       gCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGAc                                                                                <  1:6939382/150‑1 (MQ=255)

                                                                          tACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCAc                                                                             >  1:263249/1‑150 (MQ=255)

                                                                          tACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCAc                                                                             >  1:1644290/1‑150 (MQ=255)

                                                                               cgcATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCt                                                                        <  2:6039404/150‑1 (MQ=255)

                                                                               cgcATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCt                                                                        <  2:5528162/150‑1 (MQ=255)

                                                                                 cATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTgg                                                                      <  1:4589202/150‑1 (MQ=255)

                                                                                        ccGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGa                                                               <  1:1172893/150‑1 (MQ=255)

                                                                                          gggTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAAt                                                             <  1:1021436/150‑1 (MQ=255)

                                                                                          gggTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAAt                                                             <  1:4878148/150‑1 (MQ=255)

                                                                                           ggTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATg                                                            <  2:704610/150‑1 (MQ=255)

                                                                                                ccaccGAAGGTGTTGAAATATTCCATGCCATTGTT‑AAAGCTGTCTGCCACTTCCCGCGTTGTGACAACCGCAGACATTGGATGACCGTTACCAATGGGTTTGCCCATAGTGACGATGTCAGGGACAACATCCTGCAACTCGAAAGCccac                                                       >  1:3524834/5‑150 (MQ=255)

                                                                                                     gAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTg                                                  >  2:2521071/1‑150 (MQ=255)

                                                                                                     gAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTg                                                  <  1:8853095/150‑1 (MQ=255)

                                                                                                        ggTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCt                                               <  2:8134667/150‑1 (MQ=255)

                                                                                                           gttGAAATATTCCATGCCATTGTT‑AAAGCTGTCTGCCACTTCCCGCGTTGTGACAACCGCAGACATTGGATGACCGTTACCAATGGGTTTGCCCATAGTGACGATGTCAGGGACAACATCCTGCAACTCGAAAGCCCACCAGTGGCTGcc                                            >  2:9158772/2‑150 (MQ=255)

                                                                                                           gttGAAATATTCCATGCCATTGTT‑AAAGCTGTCTGCCACTTCCCGCGTTGTGACAACCGCAGACATTGGATGACCGTTACCAATGGGTTTGCCCATAGTGACGATGTCAGGGACAACATCCTGCAACTCGAAAGCCCACCAGTGGCTGcc                                            >  2:7957202/2‑150 (MQ=255)

                                                                                                            ttGAAATATTCCATGCCATTGTT‑AAAGCTGTCTGCCACTTCCCGCGTTGTGACAACCGCAGACATTGGATGACCGTTACCAATGGGTTTGCCCATAGTGACGATGTCAGGGACAACATCCTGCAACTCGAAAGCCCACCAGTGGCTGcct                                           >  1:7742547/1‑149 (MQ=255)

                                                                                                            ttGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCg                                           >  1:8239666/1‑150 (MQ=255)

                                                                                                            ttGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCAGGGTCTCGAATGCCCACCAGAGGCTGccc                                           >  2:885458/1‑149 (MQ=255)

                                                                                                             tGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGa                                          <  1:2411287/150‑1 (MQ=255)

                                                                                                              gAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGAc                                         >  1:8013415/1‑150 (MQ=255)

                                                                                                              gAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGAc                                         <  2:5089644/150‑1 (MQ=255)

                                                                                                               aaaTACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGAcc                                        >  1:5065788/1‑150 (MQ=255)

                                                                                                                aaTATTCCATGCCATTGTT‑AAAGCTGTCTGCCACTTCCCGCGTTGTGACAACCGCAGACATTGGATGACCGTTACCAATGGGTTTGCCCATAGTGACGATGTCAGGGACAACATCCTGCAACTCGAAAGCCCACCAGTGGCTGCCTAcgc                                       <  1:3967686/150‑3 (MQ=255)

                                                                                                                aaTACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGAccc                                       >  1:4647939/1‑150 (MQ=255)

                                                                                                                    cTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGTCGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGcc                                   >  2:1508085/1‑150 (MQ=255)

                                                                                                                     tCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCg                                  >  1:2943289/1‑150 (MQ=255)

                                                                                                                     tCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCg                                  >  2:4353165/1‑150 (MQ=255)

                                                                                                                        aTGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAAc                               <  1:2150716/150‑1 (MQ=255)

                                                                                                                          gCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAAccc                             >  1:4267922/1‑150 (MQ=255)

                                                                                                                            cGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGa                           <  2:1903074/150‑1 (MQ=255)

                                                                                                                                 ttgAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGa                      >  1:412387/1‑150 (MQ=255)

                                                                                                                                  tgAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGaa                     <  2:387300/150‑1 (MQ=255)

                                                                                                                                    aaaGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAAcc                   >  2:4805264/1‑150 (MQ=255)

                                                                                                                                     aaGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCt                  >  1:7607815/1‑150 (MQ=255)

                                                                                                                                      aGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTc                 <  1:4195060/150‑1 (MQ=255)

                                                                                                                                      aGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTc                 >  2:8329027/1‑150 (MQ=255)

                                                                                                                                       gC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCa                <  2:4565078/150‑1 (MQ=255)

                                                                                                                                        c‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCAt               <  2:5195120/150‑1 (MQ=255)

                                                                                                                                               cAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCa         <  1:7041551/150‑1 (MQ=255)

                                                                                                                                                   cTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGca     <  2:7654473/150‑1 (MQ=255)

                                                                                                                                                   cTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGca     >  1:322553/1‑150 (MQ=255)

                                                                                                                                                    tCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGcac    >  2:1306465/1‑150 (MQ=255)

                                                                                                                                                     cgcgcgTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGcaca   <  2:655884/150‑1 (MQ=255)

                                                                                                                                                     cgcgcgTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGcaca   >  2:6248180/1‑150 (MQ=255)

                                                                                                                                                     cgcgcgTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGcaca   <  2:5865009/150‑1 (MQ=255)

AGCAGGTAATGGCCGACACTCAAGGCGTTTTCTTTCAACTGGTCGCGCTCGATTGCATCGAGTACCGCCAGGCCTACCGCGCATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGC‑ATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAG  >  NC_002947/3803120‑3803418

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