BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,127,345	G→A	5.6%	P56P (CCG→CCA)	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,345	0	G	A	5.6%	555.9 / 14.7	179	P56P (CCG→CCA)	PP_4512	DMT superfamily permease
Reads supporting (aligned to +/- strand): ref base G (81/87); new base A (4/6); total (86/93)
Fisher's exact test for biased strand distribution p-value = 7.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

Only 100 of 186 total aligned reads displayed.

CTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCAGCTTCCCGGCCTTTA  >  NC_002947/5127196‑5127493

 tCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCgg                                                                                                                                                     >  2:2103229/1‑150 (MQ=255)

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

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

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

      gcACTGCACATTGGCGCGCGGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTgg                                                                                                                                                <  1:530992/150‑1 (MQ=255)

       cACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCCGGACCGGCCAGGCTGGc                                                                                                                                               <  2:2059860/150‑1 (MQ=255)

            cacaTTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTg                                                                                                                                                 >  2:573615/1‑143 (MQ=255)

            cacaTTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTg                                                                                                                                                 <  1:573615/143‑1 (MQ=255)

              caTTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTTGCCAGCCCGGACCGGCCAGGCTGGcgaacagg                                                                                                                                        <  2:1263565/150‑8 (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)

                  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)

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

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

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

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

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

                                   ccTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGACCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGcgccg                                                                                                                   >  2:4662992/1‑150 (MQ=255)

                                     ttACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTgg                                                                                                                                                >  1:2718516/1‑119 (MQ=255)

                                     ttACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTgg                                                                                                                                                <  2:2718516/119‑1 (MQ=255)

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

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

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

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

                                         cGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCctgct                                                                                                             >  1:2609671/1‑150 (MQ=255)

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

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

                                               cTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCgg                                                                                                       <  2:483792/150‑1 (MQ=255)

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

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

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

                                                   ggCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCggtggt                                                                                                   <  2:5336155/150‑1 (MQ=255)

                                                    gCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGgtggtg                                                                                                  <  1:3188431/150‑1 (MQ=255)

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

                                                             ccaccaGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTgg                                                                                         <  1:3805749/150‑1 (MQ=255)

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

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

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

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

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

                                                               accaGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTggcc                                                                                       <  1:3702286/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)

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

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

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

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

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

                                                                                    atcGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCttcttc                                                                  >  1:3365788/1‑150 (MQ=255)

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

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

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

                                                                                        tCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCg                                                              >  1:3241148/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)

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

                                                                                                   gcCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTg                                                   >  1:4448735/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                                          <  1:4148724/150‑1 (MQ=255)

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

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

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

                                                                                                                tACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCa                                      <  2:1434750/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:2166275/1‑150 (MQ=255)

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

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

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

                                                                                                                     gcgcTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGcc                                 >  1:1747110/1‑150 (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)

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

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

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

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

                                                                                                                                     ccAGCCCGGACCGGCCAGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCa                 >  2:4787729/1‑150 (MQ=255)

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

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

                                                                                                                                             gCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCAGCTTCCCg         >  2:1230823/1‑150 (MQ=255)

                                                                                                                                                    cGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCAGCTTCCCGGCCTTTa  <  1:3376628/150‑1 (MQ=255)

CTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCAGCTTCCCGGCCTTTA  >  NC_002947/5127196‑5127493

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