BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,127,338	C→A	5.8%	A54D (GCC→GAC)	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,338	0	C	A	5.8%	597.0 / 17.9	181	A54D (GCC→GAC)	PP_4512	DMT superfamily permease
Reads supporting (aligned to +/- strand): ref base C (84/86); new base A (5/6); total (89/92)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 185 total aligned reads displayed.

CACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCAGCTTCCCG  >  NC_002947/5127189‑5127486

caccgcCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGc                                                                                                                                                      >  1:4203240/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                                                                                                                                              >  2:2103229/1‑150 (MQ=255)

        tCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCgg                                                                                                                                              <  1:1627191/150‑1 (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                                                                                                                                          <  1:573615/143‑1 (MQ=255)

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

                     caTTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTTGCCAGCCCGGACCGGCCAGGCTGGcgaacagg                                                                                                                                 <  2:1263565/150‑8 (MQ=255)

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

                     caTTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCt                                                                                                                                 <  2:2959829/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)

                          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)

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

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

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

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

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

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

                                            ttACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTgg                                                                                                                                         >  1:2718516/1‑119 (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                                                                                                    >  2:3943004/1‑150 (MQ=255)

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

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

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

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

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

                                                       ttCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCggt                                                                                               >  1:4338619/1‑150 (MQ=255)

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

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

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

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

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

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

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

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

                                                                     caccaGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTggc                                                                                 >  1:278345/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)

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

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

                                                                                     gCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGc                                                                 <  2:4326706/150‑1 (MQ=255)

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

                                                                                            tcGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCa                                                          >  2:3610465/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)

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

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

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

                                                                                                    ggTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTc                                                  <  2:3129937/150‑1 (MQ=255)

                                                                                                     gTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCa                                                 >  2:3024823/1‑150 (MQ=255)

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

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

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

                                                                                                                       tACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCa                               <  2:1434750/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)

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

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

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

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

                                                                                                                                        gttgCCAGCCCGGACCGGCCAGGCTGGCGACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTc              >  1:2025542/4‑150 (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          <  2:1685357/150‑1 (MQ=255)

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

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

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

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

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

CACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTCTTCGGTCGCGCCACCTTCGCCGTACTCGCGCTGGCCCTGTTCGCCAGCCTGGCCGGTCCGGGCTGGCAACGCCTGAGCAGCCAGGATGCGCGCCGCCTGCTGCTCGGTGGTGTGCTGCTGGCCGGCCACTGGGTCAGCTTCTTCATCGCCGTCAAAGTGGGCGGCGTAGCCATCGCCACCCTCGGCTTCGCCAGCTTCCCG  >  NC_002947/5127189‑5127486

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