BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,596,569	T→A	5.4%	L258* (TTG→TAG) ‡	PP_3169 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,596,569	0	T	A	5.4%	1372.6 / 39.4	394	L258* (TTG→TAG) ‡	PP_3169	membrane protein
Reads supporting (aligned to +/- strand): ref base T (217/155); new base A (11/11); total (228/166)
Fisher's exact test for biased strand distribution p-value = 5.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 402 total aligned reads displayed.

GGTCTGGCGGGTGCTGGCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTGGATGGTCGGTGTGC  >  NC_002947/3596420‑3596718

ggTCTGGCGGGTGCTGGCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCtt                                                                                                                                                       <  1:5472124/150‑1 (MQ=255)

      gCGGGTGCTGGCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGAACCAGCGCAAGCGCTg                                                                                                                                                 <  1:5733691/150‑1 (MQ=255)

        gggTGCTGGCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCAAGCGCTGGt                                                                                                                                               <  1:3455419/150‑1 (MQ=255)

         ggTGCTGGCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCAAGCGCTGGtg                                                                                                                                              <  2:2373462/150‑2 (MQ=255)

             cTGGCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCAcc                                                                                                                                          >  1:4850831/1‑150 (MQ=255)

              tGGCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCg                                                                                                                                         >  1:3863750/1‑150 (MQ=255)

              tGGCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCg                                                                                                                                         >  2:2454240/1‑150 (MQ=255)

               ggCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGa                                                                                                                                        <  2:2099693/150‑1 (MQ=255)

               ggCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCAACgg                                                                                                                                        <  1:5201543/150‑2 (MQ=255)

                  cAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGa                                                                                                                                     <  1:2000903/150‑1 (MQ=255)

                  cAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGa                                                                                                                                     >  2:1970706/1‑150 (MQ=255)

                  cAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGa                                                                                                                                     >  2:6431906/1‑150 (MQ=255)

                  cAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGa                                                                                                                                     >  1:6894294/1‑150 (MQ=255)

                   aaGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGAc                                                                                                                                    >  2:103869/1‑150 (MQ=255)

                    aGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACt                                                                                                                                   <  1:6640541/150‑1 (MQ=255)

                        cccGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCt                                                                                                                               >  1:7393569/1‑150 (MQ=255)

                        cccGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCt                                                                                                                               <  1:6219537/150‑1 (MQ=255)

                        cccGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCt                                                                                                                               >  2:2962224/1‑150 (MQ=255)

                        cccGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCt                                                                                                                               <  1:5414987/150‑1 (MQ=255)

                         ccGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCACGCACTGAACGATGACCCGATGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTg                                                                                                                              >  2:6967712/1‑150 (MQ=255)

                         ccGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTg                                                                                                                              >  1:503917/1‑150 (MQ=255)

                            aaTGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGt                                                                                                                           <  2:1724680/150‑1 (MQ=255)

                             aTGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCCCCCGACATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGtt                                                                                                                          <  2:3348220/150‑1 (MQ=255)

                             aTGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGtt                                                                                                                          <  1:2135358/150‑1 (MQ=255)

                              tGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTa                                                                                                                         <  1:7279051/150‑1 (MQ=255)

                               ggCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTAt                                                                                                                        <  1:1006957/150‑1 (MQ=255)

                                 cGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGa                                                                                                                      <  2:3889486/150‑1 (MQ=255)

                                 cGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGa                                                                                                                      <  2:5788444/150‑1 (MQ=255)

                                  gATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCAAGCGCTGGtgc                                                                                                                                             >  1:408841/1‑124 (MQ=255)

                                  gATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCAAGCGCTGGtgc                                                                                                                                             <  2:408841/126‑3 (MQ=255)

                                   aTTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGAcc                                                                                                                    <  1:1688445/150‑1 (MQ=255)

                                   aTTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGAcc                                                                                                                    >  2:800091/1‑150 (MQ=255)

                                        tacGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGAtt                                                                                                               >  1:4226293/1‑150 (MQ=255)

                                        tacGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGAtt                                                                                                               >  1:5299761/1‑150 (MQ=255)

                                        tacGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGAtt                                                                                                               >  1:176906/1‑150 (MQ=255)

                                        tacGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGAtt                                                                                                               >  1:6774613/1‑150 (MQ=255)

                                        tacGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGAtt                                                                                                               >  1:1561593/1‑150 (MQ=255)

                                        tacGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCAAGCGCTg                                                                                                                                                 >  1:885412/1‑116 (MQ=255)

                                        tacGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACACGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGAtt                                                                                                               >  1:4449582/1‑150 (MQ=255)

                                           gAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGc                                                                                                            <  1:5689790/150‑1 (MQ=255)

                                                 aTCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTg                                                                                                                                                 >  2:6796476/1‑107 (MQ=255)

                                                    agcagcGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACtg                                                                                                   >  2:3621878/1‑150 (MQ=255)

                                                    agcagcGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACtg                                                                                                   <  2:3060749/150‑1 (MQ=255)

                                                    agcagcGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACtg                                                                                                   <  2:2663649/150‑1 (MQ=255)

                                                     gcagcGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACtgg                                                                                                  <  2:6730728/150‑1 (MQ=255)

                                                      cagcGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACtggt                                                                                                 >  1:1839018/1‑150 (MQ=255)

                                                      cagcGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACtggt                                                                                                 <  2:1305901/150‑1 (MQ=255)

                                                       agcGTGCTCGACCAGCAACCACCCGAAAGGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACtggtg                                                                                                <  2:6843454/150‑1 (MQ=255)

                                                        gcGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACtggtgg                                                                                               <  2:3162083/150‑1 (MQ=255)

                                                               cGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGa                                                                                        >  2:4479207/1‑150 (MQ=255)

                                                                  ccAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCt                                                                                     <  1:2580688/150‑1 (MQ=255)

                                                                   cAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTg                                                                                    >  1:6328390/1‑150 (MQ=255)

                                                                    aGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTgc                                                                                   <  2:289989/150‑1 (MQ=255)

                                                                    aGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTgc                                                                                   <  2:5986149/150‑1 (MQ=255)

                                                                       aaCCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCa                                                                                >  2:6893649/1‑150 (MQ=255)

                                                                         ccaccCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATg                                                                              <  2:3933146/150‑1 (MQ=255)

                                                                          caccCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATgg                                                                             <  2:2568062/150‑1 (MQ=255)

                                                                           accCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATggg                                                                            >  1:5959044/1‑150 (MQ=255)

                                                                            cccGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGAt                                                                                                                >  2:185131/1‑113 (MQ=255)

                                                                            cccGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATgggg                                                                           >  2:6307841/1‑150 (MQ=255)

                                                                            cccGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATgggg                                                                           >  2:4782785/1‑150 (MQ=255)

                                                                             ccGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGa                                                                                                                      <  2:1520932/106‑1 (MQ=255)

                                                                             ccGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGa                                                                          <  1:6062373/150‑1 (MQ=255)

                                                                                aaaTGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATcc                                                                       >  1:2242883/1‑150 (MQ=255)

                                                                                    gCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCggg                                                                   <  2:5649761/150‑1 (MQ=255)

                                                                                     cAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCgggg                                                                  >  1:2966161/1‑150 (MQ=255)

                                                                                     cAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCgggg                                                                  >  2:7228057/1‑150 (MQ=255)

                                                                                         gCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGTCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAAtt                                                              >  1:1317762/1‑150 (MQ=255)

                                                                                            tGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTaca                                                           <  2:1412615/150‑1 (MQ=255)

                                                                                             gTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTacac                                                          <  2:6068501/150‑1 (MQ=255)

                                                                                             gTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTacac                                                          <  1:3672229/150‑1 (MQ=255)

                                                                                              tCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGt                                                                                                                           <  1:2299903/84‑1 (MQ=255)

                                                                                                 cgcAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCttt                                                      >  1:2485828/1‑150 (MQ=255)

                                                                                                      ccTGGAACTGGCGCAAGCACTGAACGATGACCCGGTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCtt                                                 <  2:3738210/150‑1 (MQ=255)

                                                                                                        tGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCa                                               <  1:7132460/150‑1 (MQ=255)

                                                                                                        tGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCa                                               >  1:7335076/1‑150 (MQ=255)

                                                                                                             cTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATg                                          >  2:156598/1‑150 (MQ=255)

                                                                                                              tGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCCTGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATgg                                         >  2:5132838/1‑150 (MQ=255)

                                                                                                                 cgcAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTc                                      <  2:4954625/150‑1 (MQ=255)

                                                                                                                        aCTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAgccaag                               >  2:7088268/1‑146 (MQ=255)

                                                                                                                        aCTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGa                               <  2:4618638/150‑1 (MQ=255)

                                                                                                                            aaCGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGatcat                           <  1:4340893/150‑1 (MQ=255)

                                                                                                                              cGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGAtcatca                         >  1:3516572/1‑150 (MQ=255)

                                                                                                                               gATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCaa                        >  1:3229253/1‑150 (MQ=255)

                                                                                                                                  gACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAggg                     <  1:5576889/150‑1 (MQ=255)

                                                                                                                                    cccGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATg                                          >  1:6434886/1‑127 (MQ=255)

                                                                                                                                    cccGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGa                               >  1:4002504/1‑138 (MQ=255)

                                                                                                                                    cccGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCa                   <  2:5411296/150‑1 (MQ=255)

                                                                                                                                       gTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGGTATGACCCGATTGGCGACGCACTGGTGGTCCCCGACCTGCGCCTGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTg                >  2:2831282/1‑150 (MQ=255)

                                                                                                                                         tGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCACAGCGCGATCATCAAGGGCAGTGGa              <  1:1451765/150‑1 (MQ=255)

                                                                                                                                             ccAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTGGATGGt          <  2:5471741/150‑1 (MQ=255)

                                                                                                                                             ccAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTGGATGGt          <  1:2758665/150‑1 (MQ=255)

                                                                                                                                              cAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCt                                                  <  2:652107/109‑1 (MQ=255)

                                                                                                                                              cAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTGGATGGTc         >  2:429953/1‑150 (MQ=255)

                                                                                                                                              cAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGAGCGATCATCAAGGGCAGTGGATGGTc         >  1:6997124/1‑150 (MQ=255)

                                                                                                                                               aGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTGGATGGTCg        >  2:2181537/1‑150 (MQ=255)

                                                                                                                                                 cgcTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTGGATGGTCGgt      <  1:4726499/150‑1 (MQ=255)

                                                                                                                                                  gcTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTGGATGGTCGgtg     <  2:7347142/150‑1 (MQ=255)

                                                                                                                                                   cTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTGGATGGTCGgtgt    >  1:774934/1‑150 (MQ=255)

                                                                                                                                                     tGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTGGATGGTCGGTGTGc  <  1:3299418/150‑1 (MQ=255)

GGTCTGGCGGGTGCTGGCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTGGATGGTCGGTGTGC  >  NC_002947/3596420‑3596718

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