BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,264,638	C→G	5.6%	G263A (GGC→GCC)	PP_1106 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,264,638	0	C	G	5.6%	679.6 / 22.4	197	G263A (GGC→GCC)	PP_1106	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (110/76); new base G (9/2); total (119/78)
Fisher's exact test for biased strand distribution p-value = 2.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 207 total aligned reads displayed.

CGGTGAGGTAGCGCACCAGTACCTCGACCGAGCCGTGAGTGACCATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCGCAGCAGTGGATCGTTGCGCGGAATATCGCCTGCGTAGCGGG  >  NC_002947/1264489‑1264785

  gTGAGGTAGCGCACCAGTACCTCGACCGAGCCGTGAGTGACCATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGTGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGcgccg                                                                                                                                                   >  1:128328/1‑150 (MQ=255)

   tGAGGTAGCGCACCAGTACCTCGACCGAGCCGTGAGTGACCATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGcgccgc                                                                                                                                                  >  2:3296722/1‑150 (MQ=255)

     aGGTAGCGCACCAGTACCTCGACCGAGCCGTGAGTGACCATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGt                                                                                                                                                <  1:284520/150‑1 (MQ=255)

      ggTAGCGCACCAGTACCTCGACCGAGCCGTGAGTGACCATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCGGCGCCCGCCGGC‑GGCGCGGCGtg                                                                                                                                               <  1:2266159/150‑2 (MQ=255)

          gcgcACCAGTACCTCGACCGAGCCGTGAGTGACCATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTagcag                                                                                                                                           <  1:3417466/150‑1 (MQ=255)

           cgcACCAGTACCTCGACCGAGCCGTGAGTGACCATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAACCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTaccacc                                                                                                                                          <  1:1177988/150‑3 (MQ=255)

                cAGTACCTCGACCGAGCCGTGAGTGACCATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCa                                                                                                                                     >  2:3355087/1‑150 (MQ=255)

                cAGTACCTCGACCGAGCCGTGAGTGACCATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCa                                                                                                                                     >  2:3219041/1‑150 (MQ=255)

                cAGTACCTCGACCGAGCCGTGAGTGACCATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCa                                                                                                                                     >  2:1365613/1‑150 (MQ=255)

                  gTACCTCGACCGAGCCGTGAGTGACCATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGc                                                                                                                                   >  1:1416205/1‑150 (MQ=255)

                   tACCTCGACCGAGCCGTGAGTGACCATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGcc                                                                                                                                  >  1:2213592/1‑150 (MQ=255)

                      cTCGACCGAGCCGTGAGTGACCATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCt                                                                                                                               >  1:1225773/1‑150 (MQ=255)

                      cTCGACCGAGCCGTGAGTGACCATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCt                                                                                                                               >  1:3172279/1‑150 (MQ=255)

                       tCGACCGAGCCGTGAGTGACCATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTg                                                                                                                              >  2:2418979/1‑150 (MQ=255)

                       tCGACCGAGCCGTGAGTGACCATCACCCGTTGGGCACCGGACTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTg                                                                                                                              <  1:1880567/150‑1 (MQ=255)

                        cGACCGAGCCGTGAGTGACCATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTgg                                                                                                                             <  2:4172356/150‑1 (MQ=255)

                         gACCGAGCCGTGAGTGACCATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGc                                                                                                                            <  1:2380172/150‑1 (MQ=255)

                             gAGCCGTGAGTGACCATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAg                                                                                                                        <  1:249444/150‑1 (MQ=255)

                                 cGTGAGTGACCATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCAt                                                                                                                    >  2:1166617/1‑150 (MQ=255)

                                 cGTGAGTGACCATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCAt                                                                                                                    >  1:414314/1‑150 (MQ=255)

                                 cGTGAGTGACCATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCAt                                                                                                                    >  1:1776949/1‑150 (MQ=255)

                                   tgagtgaCCATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCa                                                                                                                  <  2:1636107/150‑1 (MQ=255)

                                     agtgaCCATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACt                                                                                                                <  2:2555646/150‑1 (MQ=255)

                                     agtgaCCATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACt                                                                                                                <  1:433322/150‑1 (MQ=255)

                                     agtgaCCATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACt                                                                                                                <  2:4361719/150‑1 (MQ=255)

                                      gtgaCCATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTg                                                                                                               >  1:216100/1‑150 (MQ=255)

                                       tgaCCATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGt                                                                                                              <  1:1144694/150‑1 (MQ=255)

                                         aCCATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTaa                                                                                                            >  1:597699/1‑150 (MQ=255)

                                         aCCATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTaa                                                                                                            >  2:970307/1‑150 (MQ=255)

                                           cATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATc                                                                                                          >  2:877964/1‑150 (MQ=255)

                                             tCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGc                                                                                                        <  2:442540/150‑1 (MQ=255)

                                                cccGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCaa                                                                                                     >  2:1207606/1‑150 (MQ=255)

                                                cccGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCaa                                                                                                     >  2:1338982/1‑150 (MQ=255)

                                                  cGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCaaaa                                                                                                   >  2:95409/1‑150 (MQ=255)

                                                         cACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGGCAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCa                                                                                            >  2:2824835/1‑150 (MQ=255)

                                                          aCCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCAt                                                                                           >  1:3081494/1‑150 (MQ=255)

                                                          aCCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCAt                                                                                           >  2:2659893/1‑150 (MQ=255)

                                                           ccGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATc                                                                                          <  1:1267849/150‑1 (MQ=255)

                                                            cGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATcc                                                                                         >  1:1007706/1‑150 (MQ=255)

                                                                 tGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTa                                                                                    >  1:3758544/1‑150 (MQ=255)

                                                                 tGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTa                                                                                    >  1:4489214/1‑150 (MQ=255)

                                                                     cAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTCCTGc                                                                                >  2:3230663/1‑150 (MQ=255)

                                                                     cAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGc                                                                                >  1:1433299/1‑150 (MQ=255)

                                                                      aaTGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGcc                                                                               >  2:2150310/1‑150 (MQ=255)

                                                                         ggCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCACGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCAcc                                                                            >  1:2969253/1‑150 (MQ=255)

                                                                                agCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGAc                                                                     <  1:3634474/150‑1 (MQ=255)

                                                                                  caggccaggcCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACgg                                                                   >  1:3965725/1‑150 (MQ=255)

                                                                                     gccaggccagTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACggggg                                                                >  2:3674544/1‑150 (MQ=255)

                                                                                     gccaggccagTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACggggg                                                                >  2:3077750/1‑150 (MQ=255)

                                                                                       caggccagTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGgcg                                                              >  1:3483787/1‑150 (MQ=255)

                                                                                       caggccagTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGgcg                                                              >  1:3307892/1‑150 (MQ=255)

                                                                                       caggccagTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGgcg                                                              >  2:2027225/1‑150 (MQ=255)

                                                                                         ggccagTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGcc                                                            <  2:4053071/150‑1 (MQ=255)

                                                                                         ggccagTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGcc                                                            <  2:1749467/150‑1 (MQ=255)

                                                                                         ggccagTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGcc                                                            >  1:656726/1‑150 (MQ=255)

                                                                                          gccagTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCa                                                           <  2:3096938/150‑1 (MQ=255)

                                                                                          gccagTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCa                                                           >  1:4244529/1‑150 (MQ=255)

                                                                                           ccagTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGACAg                                                          <  1:2143964/150‑1 (MQ=255)

                                                                                             agTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGa                                                        <  2:4317739/150‑1 (MQ=255)

                                                                                                       gTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTg                                              <  1:257668/150‑1 (MQ=255)

                                                                                                       gTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTg                                              >  1:3531992/1‑150 (MQ=255)

                                                                                                        tCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTgg                                             >  2:41395/1‑150 (MQ=255)

                                                                                                        tCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTgg                                             <  1:1610755/150‑1 (MQ=255)

                                                                                                        tCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTgg                                             <  1:1299156/150‑1 (MQ=255)

                                                                                                         cGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGgc                                            <  2:1512696/150‑1 (MQ=255)

                                                                                                         cGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGgc                                            <  2:3229693/150‑1 (MQ=255)

                                                                                                          ggACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGgcg                                           <  2:1002107/150‑1 (MQ=255)

                                                                                                            aCAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCgca                                         <  1:3783960/150‑1 (MQ=255)

                                                                                                               gCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCgcagca                                      <  2:607221/150‑1 (MQ=255)

                                                                                                               gCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCgcagca                                      <  2:607273/150‑1 (MQ=255)

                                                                                                                 acaAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCGCAGCAGt                                    >  1:1045807/1‑150 (MQ=255)

                                                                                                                 acaAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCGCAGCAGt                                    >  1:406969/1‑150 (MQ=255)

                                                                                                                 acaAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCGCAGCAGt                                    >  1:845443/1‑150 (MQ=255)

                                                                                                                    aaGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCGCAGCAGTGGa                                 <  2:2433722/150‑1 (MQ=255)

                                                                                                                        cccGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCGCAGCAGTGGATCGt                             <  1:1120226/150‑1 (MQ=255)

                                                                                                                         ccGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCGCAGCAGTGGATCGtt                            >  2:301890/1‑150 (MQ=255)

                                                                                                                         ccGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCGCAGCAGTGGATCGtt                            >  2:3763896/1‑150 (MQ=255)

                                                                                                                         ccGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCGCAGCAGTGGATCGtt                            >  2:3728593/1‑150 (MQ=255)

                                                                                                                          cGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCGCAGCAGTGGATCGTTg                           >  1:2792793/1‑150 (MQ=255)

                                                                                                                             tCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCCTCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCGCAGCAGTGGATCGTTgcgc                        >  2:3490790/1‑150 (MQ=255)

                                                                                                                             tCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCGCAGCAGTGGATCGTTgcgc                        >  1:96675/1‑150 (MQ=255)

                                                                                                                             tCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCGCAGCAGTGGATCGTTgcgc                        <  2:4454329/150‑1 (MQ=255)

                                                                                                                               tgctacgcg‑gcgcCCGCCGGCGGCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCGCAGCAGTGGATCGTTGCGCgg                      >  2:588333/6‑150 (MQ=255)

                                                                                                                               tgctacgcg‑gcgcCCGCCGGCGGCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCGCAGCAGTGGATCGTTGCGCgg                      >  1:3027704/6‑150 (MQ=255)

                                                                                                                                gctacgcg‑gcgcCCGCCGGCGGCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCGCAGCAGTGGATCGTTGCGCGGa                     >  2:3534670/5‑150 (MQ=255)

                                                                                                                                aCGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCGCAGCAGTGGATCGTTGCGCGGa                     >  2:3943713/1‑150 (MQ=255)

                                                                                                                                 cgccgcg‑ccgCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCGCAGCAGTGGATCGTTGCGCGGaa                    <  1:1171041/150‑1 (MQ=255)

                                                                                                                                   ccgcg‑ccgcCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCGCAGCAGTGGATCGTTGCGCGGAata                  >  1:3713418/1‑150 (MQ=255)

                                                                                                                                   acgcg‑gcgcCCGCCGGCGGCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCGCAGCAGTGGATCGTTGCGCGGAata                  >  2:2494328/2‑150 (MQ=255)

                                                                                                                                     gcg‑gcgcCCGCCGGCGGCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCGCAGCAGTGGATCGTTGCGCGGAATATc                >  2:1646/1‑150 (MQ=255)

                                                                                                                                     gcg‑gcgcCCGCCGGCGGCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCGCAGCAGTGGATCGTTGCGCGGAATATc                >  1:743541/1‑150 (MQ=255)

                                                                                                                                     gcg‑ccgcCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCGCAGCAGTGGATCGTTGCGCGGAATATc                <  1:2944761/150‑1 (MQ=255)

                                                                                                                                       g‑ccgccGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCGCAGCAGTGGATCGTTGCGCGGAATATCGc              >  1:1864051/1‑150 (MQ=255)

                                                                                                                                          cgccgGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCGCAGCAGTGGATCGTTGCGCGGAATATCGCCt            >  2:1737324/1‑150 (MQ=255)

                                                                                                                                          cgccgGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCGCAGCAGTGGATCGTTGCGCGGAATATCGCCt            >  2:3548252/1‑150 (MQ=255)

                                                                                                                                          cgccgGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCGCAGCAGTGGATCGTTGCGCGGAATATCGCCt            >  2:1270636/1‑150 (MQ=255)

                                                                                                                                           gccgGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCGCAGCAGTGGATCGTTGCGCGGAATATCGCCTg           <  1:105722/150‑1 (MQ=255)

                                                                                                                                           gccgGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACGGTAATCGCCAAAACGGCGCATCCAGCTACTGCCCCCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCGCAGCAGTGGATCGTTGCGCGGACGATCGCCTg           >  2:3439073/1‑150 (MQ=255)

                                                                                                                                                  ggcgCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCGCAGCAGTGGATCGTTGCGCGGAATATCGCCTGCGTAGCg    >  1:2686824/1‑150 (MQ=255)

                                                                                                                                                   gcgcCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGAC‑GGGGCGCCAGGACCAGTGCCTGGCGCAGCAGTGGATCGTTGCGCGGAATATCGCCTGCGTAGCggg  >  2:2754668/1‑150 (MQ=255)

CGGTGAGGTAGCGCACCAGTACCTCGACCGAGCCGTGAGTGACCATCACCCGTTGGGCACCGGTCTGGCCAATGGCCCAGAGCAGGCCAGGCCAGTCGGCGTGGTCGGACAGCACAAAGCCCCGGTCCACGCCGCG‑CCGCCGGCGGGCGCCGCGTAGCAGCATCCAGCCGCTGGCGAAGGCATCACTGTAATCGCCAAAACGGCGCATCCAGCTACTGCCACCGGCTGACGGGGGCGCCAGGACCAGTGCCTGGCGCAGCAGTGGATCGTTGCGCGGAATATCGCCTGCGTAGCGGG  >  NC_002947/1264489‑1264785

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