BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,127,335	T→C	5.4%	L53P (CTG→CCG)	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,335	0	T	C	5.4%	580.7 / 13.7	185	L53P (CTG→CCG)	PP_4512	DMT superfamily permease
Reads supporting (aligned to +/- strand): ref base T (85/89); new base C (4/6); total (90/95)
Fisher's exact test for biased strand distribution p-value = 7.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.06e-01

Only 100 of 189 total aligned reads displayed.

CCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCAGCTTC  >  NC_002947/5127186‑5127483

ccgcaccgcCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCt                                                                                                                                                      <  2:2608806/150‑1 (MQ=255)

ccgcaccgcCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCt                                                                                                                                                      >  1:2223061/1‑150 (MQ=255)

  gcaccgcCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTgg                                                                                                                                                    <  2:4899195/150‑1 (MQ=255)

   caccgcCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGc                                                                                                                                                   >  1:4203240/1‑150 (MQ=255)

   caccgcCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGc                                                                                                                                                   >  2:2550395/1‑150 (MQ=255)

   caccgcCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGc                                                                                                                                                   <  1:3105275/150‑1 (MQ=255)

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

           tCGGCGCACTGCACCTTGGCGCGCGGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCCGCGCCAGCCCGGCGATCATCGTCGTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCgg                                                                                                                                           <  1:4147079/150‑1 (MQ=255)

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

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

            cGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCggg                                                                                                                                          <  1:5012449/150‑1 (MQ=255)

               cgcACTGCACATTGGCGCGCTGTTCTTCTGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTg                                                                                                                                       <  2:2244362/150‑1 (MQ=255)

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

                gcACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCCGGACCGGCCAGGCTgg                                                                                                                                      <  2:1609131/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:2800691/150‑1 (MQ=255)

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

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

                          ttGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGa                                                                                                                            >  1:4658752/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:3847977/1‑150 (MQ=255)

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

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

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

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

                                           ggCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCCACGCCTGAGCAGCCAGGATgcgcgc                                                                                                           >  1:4639724/1‑150 (MQ=255)

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

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

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

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

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

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

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

                                                        tCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCg                                                                                              >  2:3496842/1‑150 (MQ=255)

                                                          ttCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCACCGCGTGAGCAGCGAGGATGCGCGCCGCCTGCTGCTCggt                                                                                            >  1:2354914/1‑150 (MQ=255)

                                                          ttCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCggt                                                                                            >  1:4338619/1‑150 (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)

                                                             ggCAAACTGGCCACCATCGCCAGCCCGGCGATCATCTTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGACGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCggtggt                                                                                         >  2:2478347/1‑150 (MQ=255)

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

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

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

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

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

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

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

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

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

                                                                              cgcCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCa                                                                        >  1:2538748/1‑150 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                                                                                        cGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGc                              >  1:4323675/1‑150 (MQ=255)

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

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

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

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

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

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

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

                                                                                                                                 gcTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCAc                     <  2:5087108/150‑1 (MQ=255)

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

                                                                                                                                         tGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCt             >  1:4825852/1‑150 (MQ=255)

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

                                                                                                                                           ttCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTc           <  1:416685/150‑1 (MQ=255)

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

                                                                                                                                            ttgCCAGCCCGGACCGGCCAGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCg          >  2:4466984/3‑150 (MQ=255)

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

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

                                                                                                                                              gCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGcc        >  2:1925965/1‑150 (MQ=255)

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

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

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

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

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

CCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCAGCTTC  >  NC_002947/5127186‑5127483

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