BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,127,340	G→C	5.5%	G55R (GGT→CGT) ‡	PP_4512 →	DMT superfamily permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,127,340	0	G	C	5.5%	606.8 / 17.6	182	G55R (GGT→CGT) ‡	PP_4512	DMT superfamily permease
Reads supporting (aligned to +/- strand): ref base G (84/88); new base C (4/6); total (88/94)
Fisher's exact test for biased strand distribution p-value = 7.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 183 total aligned reads displayed.

CCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCAGCTTCCCGGC  >  NC_002947/5127191‑5127488

  gcCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGt                                                                                                                                                    >  2:4561886/1‑150 (MQ=255)

     cTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCg                                                                                                                                                 >  2:610746/1‑150 (MQ=255)

      tCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCgg                                                                                                                                                <  1:1627191/150‑1 (MQ=255)

      tCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCgg                                                                                                                                                <  1:1627214/150‑1 (MQ=255)

           gcACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCCGGACCGGCCAGGCTgg                                                                                                                                           <  1:5154194/150‑1 (MQ=255)

           gcACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCCGGACCGGCCAGGCTgg                                                                                                                                           <  2:1609131/150‑1 (MQ=255)

                   caTTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCt                                                                                                                                   <  2:623701/150‑1 (MQ=255)

                   caTTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCt                                                                                                                                   <  2:2800691/150‑1 (MQ=255)

                   caTTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCCGGACCGGCCAGGCTGGcgaacagg                                                                                                                                   <  1:3445294/150‑8 (MQ=255)

                     ttGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGa                                                                                                                                 >  1:4658744/1‑150 (MQ=255)

                     ttGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGa                                                                                                                                 >  1:5442035/1‑150 (MQ=255)

                       ggCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGACAACGCCTGagc                                                                                                                               >  2:3161571/1‑150 (MQ=255)

                        gcgcgcTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGagca                                                                                                                              <  2:1747379/150‑1 (MQ=255)

                        gcgcgcGGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGagca                                                                                                                              <  2:824907/150‑1 (MQ=255)

                          gcgcTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGagcagc                                                                                                                            >  2:5008272/1‑150 (MQ=255)

                             cgggtcgtcGGCCTGACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAg                                                                                                                         <  2:1241513/143‑1 (MQ=255)

                              tGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAgg                                                                                                                        >  2:1586338/1‑150 (MQ=255)

                              tGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAgg                                                                                                                        >  1:746508/1‑150 (MQ=255)

                              tGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAgg                                                                                                                        >  1:3566079/1‑150 (MQ=255)

                                ttcttgGGCCTTACCGGCTTCTTGGGCAAACGGGCCACCAGCGCCAGCCCGGCGATCATCTTCTTGGGTCGCGCCACCTTCGCCTTACTCGCGCTGGCCCTTTTCGCAAGCCTGGCGGGTCCGGGGTGGCAACGCCTGAGCAGCCAGGAt                                                                                                                      <  2:2123236/150‑1 (MQ=255)

                                  cttcGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATgc                                                                                                                    <  1:784774/150‑1 (MQ=255)

                                   ttcGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATgcg                                                                                                                   >  1:3251117/1‑150 (MQ=255)

                                    tcGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATgcgc                                                                                                                  <  1:4691651/150‑1 (MQ=255)

                                          ttACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCgccgcc                                                                                                            >  2:4573087/1‑150 (MQ=255)

                                             ccGGCGTCTTCGGCAAACTGGCCACCCGCGCCAGCCCGGCGATCATCGTCGTCGGTCGCGCCACCTTCGCCGTACTCGCGCGGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCctgc                                                                                                         <  2:2505187/150‑1 (MQ=255)

                                             ccGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCctgc                                                                                                         <  2:2475528/150‑1 (MQ=255)

                                                gCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCtgctgc                                                                                                      >  2:3943004/1‑150 (MQ=255)

                                                gCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCtgctgc                                                                                                      <  1:1222360/150‑1 (MQ=255)

                                                 cGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTgctgct                                                                                                     <  1:1662553/150‑1 (MQ=255)

                                                     ttCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCggt                                                                                                 <  1:400457/150‑1 (MQ=255)

                                                     ttCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCggt                                                                                                 <  2:3621239/150‑1 (MQ=255)

                                                      tCGGCAAACTGGCCGCCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCggtg                                                                                                <  1:589925/150‑1 (MQ=255)

                                                       cGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCggtgg                                                                                               >  1:1545967/1‑150 (MQ=255)

                                                        ggCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCggtggt                                                                                              >  1:673608/1‑150 (MQ=255)

                                                         gCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGgtggtg                                                                                             >  2:5354674/1‑150 (MQ=255)

                                                         gCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGgtggtg                                                                                             >  1:4801057/1‑150 (MQ=255)

                                                          cAAACGGGCCACCAGCGCCAGCCCGGCGAGCATCGTCGTCGGGCGCGCCACCTTCGCCGTACACGCGCAGGCCCCGTGCGCCAGCCGGGCCGGGCCGGGCTGGCCACGCCTGAGCAGCCCGGAGGCGCGCCGCCTGCTCCTCGgtggtgg                                                                                            >  2:2771562/1‑149 (MQ=255)

                                                                 gCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGctgctg                                                                                     <  1:2796961/150‑1 (MQ=255)

                                                                   caccaGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTggc                                                                                   >  1:4326706/1‑150 (MQ=255)

                                                                   caccaGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTggc                                                                                   >  1:3008050/1‑150 (MQ=255)

                                                                    accaGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTggcc                                                                                  >  2:1311666/1‑150 (MQ=255)

                                                                    accaGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTggcc                                                                                  <  2:673608/150‑1 (MQ=255)

                                                                     ccaGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTggccg                                                                                 <  1:4388886/150‑1 (MQ=255)

                                                                       aGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTggccggc                                                                               <  2:2943734/150‑1 (MQ=255)

                                                                       aGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTggccggc                                                                               >  2:1896478/1‑150 (MQ=255)

                                                                        gcgcCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTggccggcc                                                                              >  1:3824617/1‑150 (MQ=255)

                                                                             aGCCCGGCGATCATCGTCTTCGGTCGCGCCAGCTTCGGCGTACTCGCGCTGGCCGTGTTCGCCAGCCTGGCCGGTCCGGGCTTGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTgg                                                                         <  1:3933261/150‑1 (MQ=255)

                                                                             aGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTgg                                                                         <  1:2728879/150‑1 (MQ=255)

                                                                              gCCCGGCGATCATTGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGACTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTggg                                                                        >  2:3825934/1‑150 (MQ=255)

                                                                                ccGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTc                                                                      <  1:3627329/150‑1 (MQ=255)

                                                                                ccGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTc                                                                      >  2:2064514/1‑150 (MQ=255)

                                                                                     gATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTg                                                                          <  2:2236855/141‑1 (MQ=255)

                                                                                      atcatcGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGcttc                                                                <  1:319135/150‑1 (MQ=255)

                                                                                      atcatcGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGcttc                                                                <  1:3796732/150‑1 (MQ=255)

                                                                                      atcatcGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGACAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGcttc                                                                <  1:3161571/150‑1 (MQ=255)

                                                                                        catcGTCTTCGGTCGCGCCACCGTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGcttctt                                                              <  1:459194/150‑1 (MQ=255)

                                                                                          tcGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCa                                                            >  1:1394532/1‑150 (MQ=255)

                                                                                            gTCTTCGGTCGCGCCACGTTCGCCGTACTCGCGCTGGCCGTGTTCGCCAGCCTGGCCGGTCCGGGCTGGGAACGCCTGTGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTTCTGCTGGCAGGCCACTGTGTCAGCTTCTTGATc                                                          <  2:5267655/150‑1 (MQ=255)

                                                                                            gTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATc                                                          <  1:127439/150‑1 (MQ=255)

                                                                                             tCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCg                                                         >  2:4968330/1‑150 (MQ=255)

                                                                                              cTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGc                                                        <  2:3288191/150‑1 (MQ=255)

                                                                                               ttCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGcc                                                       >  1:96001/1‑150 (MQ=255)

                                                                                                tCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCg                                                      >  2:2337335/1‑150 (MQ=255)

                                                                                                 cGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTggccggc                                                                               <  2:4625078/124‑1 (MQ=255)

                                                                                                 cGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTggccggc                                                                               >  1:4625078/1‑124 (MQ=255)

                                                                                                  ggTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTc                                                    >  2:3902464/1‑150 (MQ=255)

                                                                                                    tCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCaa                                                  >  1:2349700/1‑150 (MQ=255)

                                                                                                     cgcgCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCaaa                                                 <  2:3268256/150‑1 (MQ=255)

                                                                                                     cgcgCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCaaa                                                 <  1:2326185/150‑1 (MQ=255)

                                                                                                       cgcCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGt                                               >  2:4825891/1‑150 (MQ=255)

                                                                                                        gcCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTg                                              >  2:2981006/1‑150 (MQ=255)

                                                                                                         ccACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTgg                                             <  2:3566079/150‑1 (MQ=255)

                                                                                                            ccTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGggcg                                          <  1:2550395/150‑1 (MQ=255)

                                                                                                            ccTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGggcg                                          <  1:1784350/150‑1 (MQ=255)

                                                                                                                 gCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGgcggcgat                                     <  2:746508/150‑3 (MQ=255)

                                                                                                                 gCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTa                                     <  2:5011985/150‑1 (MQ=255)

                                                                                                                 gCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTa                                     <  2:5011897/150‑1 (MQ=255)

                                                                                                                 gCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTa                                     <  1:4148724/150‑1 (MQ=255)

                                                                                                                     tACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCa                                 <  1:4109547/150‑1 (MQ=255)

                                                                                                                       cTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATc                               <  1:5008272/150‑1 (MQ=255)

                                                                                                                        tCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGtgc                                                                                          <  1:2447534/90‑1 (MQ=255)

                                                                                                                        tCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCg                              <  2:938603/150‑1 (MQ=255)

                                                                                                                        tCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCg                              >  2:2166275/1‑150 (MQ=255)

                                                                                                                          gcgcTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTgctgct                                                                                      <  1:1997075/92‑1 (MQ=255)

                                                                                                                          gcgcTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTgctgct                                                                                      >  2:1997075/1‑92 (MQ=255)

                                                                                                                          gcgcTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGcc                            >  2:4283864/1‑150 (MQ=255)

                                                                                                                              tGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCAccc                        >  1:4182898/1‑150 (MQ=255)

                                                                                                                                  ccTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCgg                    >  2:4505206/1‑150 (MQ=255)

                                                                                                                                   cTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGc                   >  2:3525147/1‑150 (MQ=255)

                                                                                                                                     gTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCtt                 <  1:2064514/150‑1 (MQ=255)

                                                                                                                                      gttgCCAGCCCGGACCGGCCAGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTc                >  1:4081977/4‑150 (MQ=255)

                                                                                                                                        cGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGc              <  2:4322956/150‑1 (MQ=255)

                                                                                                                                         gCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCATCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGcc             <  2:4207446/150‑1 (MQ=255)

                                                                                                                                          ccAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCa            <  2:4338619/150‑1 (MQ=255)

                                                                                                                                          ccAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCa            >  1:100327/1‑150 (MQ=255)

                                                                                                                                          ccAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCa            >  1:2700138/1‑150 (MQ=255)

                                                                                                                                           cAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCAg           >  2:3126022/1‑150 (MQ=255)

                                                                                                                                            aGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCAGc          >  1:5426485/1‑150 (MQ=255)

                                                                                                                                               cTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCAGCTTc       >  1:5187446/1‑150 (MQ=255)

                                                                                                                                                tGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCg                              <  2:4730680/126‑1 (MQ=255)

CCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCAGCTTCCCGGC  >  NC_002947/5127191‑5127488

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